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Computer  aided  manufacturing  plays  a significant  role  in 
the  assembly  of  printed  circuit  (PC)  boards.  In  this  work  I 
develop  a model  of  a PC  board  assembly  problem  for  a class 
of  automated  component  placement  machines.  These  machines 
consist  of  a moving  table,  a moving  carrier,  and  a robotic 
arm  or  a turret.  A supply  of  parts  is  arranged  and  delivered 
by  the  carrier  to  the  turret.  The  table  positions  the  PC 
board  at  a prespecified  location  for  the  part  placement.  The 
PC  board  assembly  problem  is  divided  into  two  subproblems: 
the  placement  sequence  problem  and  the  sequence  mending 
problem.  The  placement  sequence  problem  determines  the  order 
in  which  a single  part  type  should  be  placed  on  the  PC 
board.  It  is  formulated  as  a Hamiltonian  cycle  problem  in  a 
Tchebyshev  grid  graph  (in  this  context  also  called  the 
occupancy  graph).  An  optimal  polynomial  time  algorithm  is 
presented  for  a special  case  of  Tchebyshev  (4,  n)-grid 


graphs.  Another  algorithm,  also  optimal  and  polynomial  in 
performance,  is  presented  for  the  rectilinear  (4,  n)-grid 
graph.  The  solution  of  the  placement  sequence  problem  is 
used  in  a greedy  heuristic  to  solve  the  sequence  mending 
problem.  The  solution  of  the  sequence  mending  problem  is  a 
grand  assembly  sequence  that  determines  both  the  order  of 
placement  of  part  types  on  the  PC  board  and  the  order  in 
which  the  parts  should  by  arranged  on  the  carrier. 
Computational  experience  with  the  algorithms  is  presented. 

The  machinery  being  modelled  has  a high  degree  of 
concurrency.  This  makes  the  PC  board  assembly  problem  quite 
difficult.  My  approach  to  the  solution  of  the  problem  takes 
advantage  of  the  special  characteristics  of  the  machines, 
and,  using  a 'divide  and  conquer'  strategy,  provides  simple 
heuristic  solutions. 


IX 


CHAPTER  1 
INTRODUCTION 


Computer  aided  manufacturing  plays  a significant  role 
in  the  assembly  of  printed  circuit  (PC)  boards.  The 
development  of  flexible  manufacturing  systems  with  a high 
degree  of  concurrency  in  their  functional  operations  enables 
manufacturers  to  increase  dramatically  the  rate  of 
production.  However,  due  to  the  high  concurrency  of  the 
machinery  an  efficient  use  of  such  manufacturing  systems  can 
benefit  from  the  formulation  and  solution  of  complex 
optimization  problems.  The  concurrency  of  the  manufacturing 
system  must  be  properly  understood,  otherwise  the  high 
"burst  rate"  claimed  by  the  manufacturing  systems  can  be 
severely  degraded  in  actual  production  (Ahmadi,  Grotzinger, 
and  Johnson  [1988]). 

In  this  work  I consider  a problem  arising  in  the 
manufacturing  of  PC  boards.  During  assembly  N parts  are  to 
be  placed  on  a rectangular  board  by  an  automatic  robot. 

Since  such  an  assembly  is  repeated  a large  number  of  times, 
it  is  desirable  to  minimize  the  time  required  to  populate  a 
single  board.  Two  factors  influence  the  total  assembly  time: 
the  part  retrieval  from  the  supply  source  and  the  placement 
sequence  of  the  parts.  Thus,  the  PC  board  assembly  problem 
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consists  of  two  related  subproblems.  The  parts  have  to  made 
available  to  the  assembly  device  in  an  efficient  fashion, 
and  an  efficient  placement  sequence  must  be  known. 

The  component  placement  machine  generally  consists  of 
three  subsystems:  a part  feeder,  carrier  or  magazine,  a 
rotating  robotic  turret  or  an  arm,  and  a table.  The  parts 
are  arranged  in  some  type  of  a feeder,  carrier  or  magazine 
according  to  the  part  type.  The  feeder  aligns  the  needed 
part  with  a pick-up  position,  where  a head  on  a rotating 
robotic  turret  picks  up  the  part  and  moves  with  it  to  the 
placement  position.  Meanwhile,  the  table  aligns  the  part 
location  on  the  PC  board  with  the  placement  position  of  the 
head  in  order  for  the  head  to  place  the  part.  In  some 


z-direction 

carrier 


f 


pick-up  position 


turret 


placement  position 


y-direction 


x-direction 


table 


Figure  1.1.  Schematic  top  view  of  the  Panasonic  Mkl  machine 
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machines  the  table  may  be  stationary,  and  the  robotic  arm 
moves  to  the  proper  placement  location  on  the  PC  board. 

Since  neither  the  part  in  the  magazine  nor  the  part  location 
on  the  PC  board  are  instantaneously  reachable,  the 
arrangement  of  parts  in  the  magazine  influences  the 
placement  sequence  and  vice  versa.  Often  the  robotic  turret 
or  arm  moves  through  several  positions  between  the  pick-up 
and  placement  positions,  which  further  compounds  the 
problem. 

Heuristic  solutions  may  be  devised  by  considering  the 
characteristics  of  a particular  type  of  component  placement 
machine,  such  as  the  Panasonic  Mkl,  and  using  additional 
information  in  modelling  the  PC  board  assembly  problem.  In 
the  case  of  the  Panasonic  Mkl  component  placement  machine,  I 
obtain  an  interesting  combinatorial  optimization  problem.  I 
decompose  the  PC  board  assembly  problem  into  assembly  by 
part  type.  The  placement  sequence  problem  for  each  part  type 
can  be  formulated  as  a Hamiltonian  cycle  problem  on  a 
"Tchebyshev"  grid  graph.  Given  solutions  for  each  part  type, 
the  overall  PC  assembly  problem  can  then  be  solved  using  a 
traveling  salesman  problem  "mending"  technique,  where  cycles 
corresponding  to  different  part  types  need  to  be  mended 
together  into  a single  grand  assembly  sequence. 

I characterize  some  of  the  combinatorial  properties  of 
"Tchebyshev"  grid  graphs,  and  develop  an  algorithm  to  solve 
the  placement  sequence  problem  under  the  Mkl  conditions . I 


then  present  a mending  heuristic  to  solve  the  PC  board 
assembly  problem. 
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Grid  graphs  closely  related  to  the  "Tchebyshev"  grid 
graphs  appear  in  the  literature,  and  the  Hamiltonian  cycle 
problem  in  the  general  case  is  found  to  be  NP-complete.  I 
show  that  2-connectedness  is  a necesssary  and  sufficient 
condition  for  a special  case  of  a Tchebyshev  grid  graphs  to 
be  Hamiltonian.  This  new  characteristic  leads  to  a 
polynomial-time  Hamiltonian  cycle  algorithm  for  this  special 
case  of  grid  graphs. 

Chapter  by  chapter  my  work  proceeds  as  follows.  In 
Chapter  2,  I present  a survey  of  literature  related  to  PC 
board  assembly  and  Hamiltonian  cycles,  and  formulate  the  PC 
board  assembly  problem  for  the  Panasonic  Mkl,  and,  in 
Chapter  3,  I show  that  the  Hamiltonian  cycle  formulation  of 
the  placement  sequence  subproblem  of  the  PC  board  assembly 
problem  is  NP-complete. 

In  Chapter  4,  I show  that  a Tchebyshev  (4,  n)-grid 
graph  is  Hamiltonian  if  and  only  if  it  is  2-connected  and 
present  a polynomial-time  algorithm  based  on  this  result.  I 
also  present  some  related  Hamiltonian  cycle  results  for 
special  case  of  rectilinear  grid  graphs  in  Chapter  4.  In  the 
case  of  the  Hamiltonian  cycle  problem  has  no  solution,  a 
shortest  Tchebyshev  traveling  salesman  tour  is  a solution  to 
the  placement  sequence  problem.  In  Chapter  5,  I develop  an 
'efficient  set'  heuristic  solution  to  the  Tchebyshev 
traveling  salesman  problem  (TSP).  This  heuristic  is  a 
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modification  of  the  convex  hull  heuristic  for  the  TSP  where 
the  convex  hull  is  replaced  by  an  efficient  set  (Francis, 
McGinnis,  and  White  [1992]). 

In  Chapter  6,  I present  a mending  sequence  heuristic  to 
find  the  grand  assembly  sequence.  I also  describe  a computer 
implementation  of  the  heuristic  approach  developed  in  this 
dissertation  and  examine  its  computational  performance. 
Finally,  I propose  future  research. 


CHAPTER  2 

PC  BOARD  ASSEMBLY  PROBLEM 


In  this  chapter  I review  the  relevant  literature  on  the 
PC  board  assembly  problem  and  Hamiltonian  cycles.  The 
proposed  method  of  solution  exploits  the  characteristics  of 
the  specific  component  assembly  machine,  a Panasonic  Mkl, 
which  is  also  described  in  this  chapter.  The  PC  board 
assembly  problem  on  the  Panasonic  Mkl  is  formulated  in  the 
last  section  of  this  chapter.  The  approach  is,  however,  not 
limited  to  a specific  machine,  as  other  machines  operate  in 
a similar  manner. 


2 . 1 Literature  Review 

In  this  section  I present  a review  of  literature 
related  to  PC  board  assembly  and  Hamiltonian  cycles. 

2.1.1  PC  Board  Assembly  Literature 

The  increasing  automatization  in  the  production  of  PC 
boards  presents  several  optimization  problems  to  be  solved 
simultaneously  to  achieve  productivity  gains.  The  literature 
also  manifests  the  use  of  many  different  machines  used  in 
the  electronics  industry.  Even  though  the  general  analysis 
is  similar  for  all  of  the  machines,  the  diverse 
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characteristics  of  the  individual  machines  are  studied  and 
exploited  by  some  authors . 

Bedini,  Lisini  and  Sterpos  [1979]  consider  a specific 
industrial  robot  equipped  with  two  independent  arms . The 
robot  can  be  used  for  drilling  or  automatic  assembly.  The 
work  cycle  consists  of  a sequence  of  operations  defined  by 
location,  tool  and  time.  The  problem  is  to  minimize  total 
production  time  including  setup,  tool  retrieval,  and  delays 
due  to  arm  interference.  The  constraints  are  related  to  the 
characteristics  of  the  robot  and  layout  of  the  tools.  An 
assignment  problem  is  solved  as  the  key  part  of  the 
heuristic  solution  procedure. 

Drezner  and  Nof  [1984]  analyze  the  pick-move-insert 
sequence  in  robotic  assembly.  Their  objective  is  to 
determine  the  physical  arrangement  of  assembly  cell 
components  with  the  objective  to  optimize  pick,  place,  and 
insert  motions . They  decompose  the  problem  into  an 
assignment  problem  to  determine  location  of  parts,  and  a 
traveling  salesman  problem  (TSP).  After  the  assignment 
problem  is  solved,  the  TSP  solution  determines  the  insertion 
sequence.  In  a similar  application,  Walas  and  Askin  [1984] 
formulate  the  problem  as  a combination  of  a TSP  and  a 
quadratic  assignment  problem.  Since  both  of  these  problems 
are  solved  independently,  the  authors  suggest  an  iterative 
approach. 

In  two  related  papers  Sarin  and  Wilhelm  [1984]  and 
Wilhelm  and  Sarin  [1985]  provide  a structure  for  analysis  of 
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problems  involved  in  sequencing  operations  and  design  of 
layouts  for  a robot  tending  to  machines  in  a robotic  cell. 
Sarin  [1987]  models  the  sequencing  problem  with  specific 
parameters  as  a linear  0-1  integer  program. 

Lofgren  and  McGinnis  [1986a]  and  Lofgren  and  McGinnis 
[1986b]  consider  sequencing  of  batches  of  jobs  that  have  to 
pass  through  automated  assembly.  Each  type  of  job  requires 
different  parts  or  different  assembly  tools.  The  authors 
consider  assignment  of  tools  to  assembly  machines,  or 
sequencing  batches  through  one  assembly  machine  to  minimize 
the  set-up  times. 

Tang  [1986]  studies  the  job  scheduling  problem 
considered  by  Lofgren  and  McGinnis.  A given  number  of  job 
batches  are  assembled  on  one  machine  with  a limited  tool 
magazine.  Tang  presents  two  integer  programming  formulations 
with  different  criteria  to  be  optimized.  He  gives  a near- 
optimal  heuristic  to  solve  one  problem  and  an  optimal 
algorithm  to  solve  the  other. 

Chauny,  Haurie,  Wagneur,  and  Loulou  [1987]  use 
Bartholdi's  and  Platzman's  [1982,  1988]  idea  of  a space- 
filling curve  to  develop  a heuristic  to  solve  the  sequencing 
of  operations  to  populate  a product.  The  assignment  of  parts 
in  a carrier  is  not  considered. 

Ball  and  Magazine  [1988]  consider  assembly  with  parts 
dispensed  from  bins  on  one  side  of  the  PC  board.  They 
formulate  the  sequencing  problem  as  a rural  postman  problem. 
They  present  a linear  order  algorithm,  which  is  optimal  with 


rectilinear  distances,  and  has  a constant  performance  bound 
with  Tchebyshev  distances. 
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Ahmadi,  Grotzinger,  and  Johnson  [1988]  model  a highly 
concurrent  dual  delivery  placement  machine  used  by  IBM.  The 
authors  identify  three  optimization  problems  related  to  the 
productivity  increase  on  this  machine.  They  address  the 
allocation  of  components  to  the  carriers  and  the  reel 
sequencing  problem.  The  authors  emphasize  accuracy  of  the 
developed  model,  and  formulate  a mixed  integer  program 
solved  by  standard  OR  methods.  The  third  problem  models  the 
placement  of  fixtures  on  the  component  carrier.  It  is 
formulated  as  a nonlinear  integer  program. 

Younis  and  Cavalier  [1989]  address  the  problem  of 
locating  bins  containing  parts  in  a constrained  area.  The 
authors  use  Lagrangian  relaxation  to  solve  a 0-1  mixed 
integer  program. 

Bard,  Clayton,  and  Feo  [1989],  Leipala  and  Nevalainen 
[1989],  and  Hu  and  Carter  [1989]  consider  PC  board  assembly 
in  an  environment  where  parts  located  in  a moving  carrier  or 
carousel  are  delivered  by  a robot  to  a proper  position  on 
the  PC  board.  A table  moving  in  two  directions 
simultaneously  ensures  proper  location.  Because  of  the  part 
feeder  movement,  besides  the  insertion  sequence  (TSP)  a 
quadratic  integer  problem  must  be  solved  to  determine  the 
assignment  and  order  of  retrieval  of  parts  in  the  feeding 
mechanism.  Given  an  insertion  sequence,  Bard,  Clayton,  and 
Feo  [1989]  use  a Lagrangian  relaxation  scheme  to  jointly 
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solve  the  assignment  and  retrieval  problems.  Based  on  the 
solution,  they  suggest  resequencing  the  table  tour. 

Leipala  and  Nevalainen  [1989]  model  the  PC  board 
assembly  problem  as  a generalization  of  a traveling  salesman 
problem  and  a quadratic  assignment  problem.  They  suggest  a 
method  of  successively  fixing  the  insertion  sequence  and  the 
assignment  of  parts  to  the  feeders  to  find  an  approximate 
solution  of  the  PC  assembly  problem. 

Hu  and  Carter  [1989]  suggest  a similar  approach  as 
Leipala  and  Nevalainen  [1989]  for  another  type  of  machine 
with  a rotating  carousel  as  a part  feeder  device.  The 
circular  nature  of  the  feeder  allows  them  to  model  the 
assignment  problem  as  a "linear  arrangement  problem"  which 
permits  faster  solution  procedures  than  the  ordinary 
quadratic  assignment  problem.  They  include  a statistical 
evaluation  of  their  heuristic  as  well  as  comparison  of  their 
solutions  with  the  real  production. 

Francis,  Lee,  Hamacher,  and  Yerelan  [1994]  consider  the 
repetitive  placement  of  n parts  on  a rectangular  workpiece. 
They  formulate  a problem  to  find  locations  for  bins 
containing  the  part  supply  and  the  placement  sequence  to 
minimize  the  total  assembly  time.  The  authors  model  this 
problem  as  a traveling  salesman  problem  with  special 
structure  which  enables  the  computation  of  the  lower  bound 
on  the  minimum  total  assembly  time.  The  authors  describe  an 
asymptotically  optimal  algorithm  of  the  order  0(n  log  n); 
they  also  give  closed-form  expressions  for  the  expected 


value  of  the  lower  bound  in  the  case  when  the  parts  are 
uniformly  distributed  on  the  workpiece. 
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2.1.2 Hamiltonian  Cycle  and  Geometric  Traveling  Salesman 

Literature 

I will  formulate  the  assembly  sequence  problem  as  a 
Hamiltonian  cycle  problem  in  a special  class  of  graphs.  In 
this  section  I present  a review  of  literature  related  to 
Hamiltonian  cycles. 

Since  necessary  and  sufficient  conditions  for  a general 
graph  to  be  Hamiltonian  have  not  been  found,  Erdos  and  Renyi 
[1960]  posed  the  following  question.  If  n and  m are  the 
number  of  vertices  and  edges,  respectively,  how  fast  must 
m(n)  grow  in  order  to  make  "almost  all"  graphs  with  n 
vertices  and  m(n)  edges  Hamiltonian?  "Almost  all"  in  this 
context  means  that  the  existence  of  a Hamiltonian  cycle  is 
ensured  with  probability  approaching  1.  In  the  search  for  an 
answer  to  this  question,  Komlos  and  Szemeredi  [1973] 
developed  a so-called  extension-rotation  operation.  Korsunov 
[1976]  designed  an  efficient  algorithm  using  the  extension- 
rotation  operation  to  find  Hamiltonian  cycles  in  "almost 
all"  graphs  with  n vertices  and  3 n log(n)  edges.  Angluin 
and  Valiant  [1979]  presented  a probabilistic  algorithm  based 
on  the  extension-rotation  operation. 

Chvatal  [1973]  gives  three  necessary  conditions  for  a 
graph  to  be  Hamiltonian.  If  the  graph  G is  not  "1-tough" 
then  G is  not  Hamiltonian;  if  the  graph  G has  no  "2-factor" 
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then  G is  not  Hamiltonian,  and  if  G is  not  "3-cyclable"  then 
G is  not  Hamiltonian.  For  a more  detailed  treatise  of 
necessary  conditions  see  Lawler,  Lenstra,  Rinnooy  Kan  and 
Shmoys  [1985].  This  reference  also  contains  a treatise  on 
sufficient  conditions. 

Karp  [1972]  shows  that  the  general  Hamiltonian  cycle 
problem  is  NP-complete.  Itai,  Papadimitriou  and  Szwarcfiter 
[1982]  show  that  the  Hamiltonian  cycle  problem  remains  NP 
complete  in  (Euclidean)  grid  graphs.  Their  work  also  appears 
in  Lawler,  Lenstra,  Rinnooy  Kan  and  Shmoys  [1985].  Garey, 
Graham  and  Johnson  [1976]  show  that  the  related  geometric 
traveling  salesman  problem  with  cities  located  on  the  points 
of  the  plane  with  integer  coordinates  using  either  the 
discretized  Euclidean,  rectilinear  or  Tchebyshev  metric  is 
NP-complete. 

A natural  approach  to  solve  the  placement  problem  is  a 
geometric  type  of  a TSP  heuristic.  Bozer,  Schorn,  and  Sharp 
[1990]  study  the  computational  performance  of  several  such 
heuristics  for  the  Tchebyshev  TSP  in  a material  handling 
context.  They  include  the  convex  hull  heuristic,  the  band 
heuristic,  the  band  insertion  heuristic,  the  sweep 
heuristic,  and  the  spacefilling  heuristic.  The  authors  also 
study  the  performance  of  these  heuristics  with  a tour 
improvement  routine  known  as  2 -opt  introduced  by  Croes 
[1958]  and  generalized  by  Lin  [1965].  The  authors  conclude 
that  the  hull  heuristic  has  the  best  overall  performance  in 
tour  length  although  it  is  the  most  complicated  heuristic  on 
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a relative  scale.  The  improvement  routine  is  not  effective 
with  the  hull  heuristic,  which  generates  optimal  or  near- 
optimal  solutions.  A modification  of  the  band  insertion 
heuristic,  the  1/2  band  insertion  heuristic,  has  a 
comparable  performance.  The  1/2  band  insertion  heuristic 
coupled  with  a 2-opt  routine  yields  solutions  within  one 
percent  of  the  hull  heuristic  solutions  and  runs  1.34  to 
1.93  times  faster.  Without  improvement  the  solutions  are 
within  six  percent  and  the  1/2  band  insertion  heuristic  runs 
2.5  to  2.8  times  faster.  The  remaining  heuristics  coupled 
with  an  improvement  routine  yield  tours  comparable  in  length 
at  the  expense  of  longer  running  time. 

Francis,  Hamacher,  Lee,  and  Yerelan  [1994]  describe  a 
robot  tour  problem  and  propose  a traveling  salesman  type 
asymptotically  optimal  heuristic.  The  robot  tour  problem  is 
another  type  of  the  PC  board  assembly  problem  where  the 
robotic  arm  has  to  travel  to  bins  to  pick  up  parts  to  be 
placed  next.  The  bins  are  located  outside  of  the  rectangular 
area  occupied  by  the  PC  board.  The  clock  heuristic  consists 
of  two  phases.  During  the  first  phase,  the  part  locations 
are  projected  onto  the  closest  side  of  the  rectangle 
occupied  by  the  PC  board.  The  clock  sequence  is  determined 
in  the  second  phase  by  following  the  projected  points  in  a 
clockwise  fashion  around  the  PC  board.  This  heuristic  is 
similar  in  nature  to  the  band  heuristic  discussed  above. 

Bartholdi,  Platzman,  Collins,  and  Warden  [1983] 
describe  a problem  related  to  the  placement  sequence 
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problem.  The  authors  successfully  implemented  a new 
traveling  salesman  heuristic  based  on  the  idea  of 
spacefilling  curves  (e.g.  Bartholdi  and  Platzman  [1982]).  A 
spacefilling  curve  can  be  imagined  as  a tour  that  visits 
every  point  of  a unit  square.  A placement  sequence  (or 
traveling  salesman  tour)  is  determined  by  the  order  in  which 
locations  appear  on  the  spacefilling  curve.  The  great 
advantage  of  the  spacefilling  curve  approach  is  its 
flexibility.  If  a particular  spacefilling  curve  is 
determined  for  the  rectangular  grid  graph,  the  placement 
sequence  for  each  part  type  can  be  easily  determined.  With  a 
large  number  of  different  PC  board  designs  this  can  be  a 
very  practical  characteristic.  This  method,  however,  does 
not  take  advantage  of  the  particular  structure  of  a grid 
graph  - this  structure  should  be  exploited  in  the 
construction  of  the  spacefilling  curve. 

If  the  underlying  grid  graph  structure  is  disregarded, 
the  placement  sequence  problem  is  equivalent  to  the 
rectilinear  or  Tchebyshev  TSP.  Allison  and  Noga  [1984]  note 
that  the  rectilinear  TSP  has  not  received  the  same  amount  of 
attention  as  the  related  Euclidean  version  of  the  problem. 
The  authors  note  that  the  convex  hull  heuristic  performs 
consistently  well  for  the  Euclidean  TSP  and  present  an 
adaptation  of  this  heuristic  for  the  rectilinear  TSP.  From 
their  empirical  analysis,  the  expected  tour  lengths  obtained 
by  the  rectilinear  hull  heuristic  are  better  than  those  of 
the  nearest  and  farthest  insertion  heuristics. 
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2.2  Problem  Description 

In  this  section  I describe  the  characteristics  of  the 
Panasonic  Mkl  component  placement  machine  that  can  be 
exploited  in  modelling  the  PC  board  assembly  problem.  The  PC 
board  assembly  problem  on  the  Panasonic  Mkl  is  formulated  in 
the  last  subsection. 

2.2.1  Machine  Description 

The  Panasonic  Mkl  machine  consists  of  three  subsystems: 
a table,  a translating  carrier  with  a supply  of  parts,  and  a 
rotating  robotic  turret  with  10  heads.  All  three  parts  of 
the  machine  operate  concurrently. 

The  turret  heads  pick  specified  parts  at  a designated 
pick  up  position,  where  the  parts  are  delivered  by  the 
carrier.  They  then  move  the  parts  to  a placement  position. 
The  PC  board  is  positioned  on  a table  that  can  move  in  two 
orthogonal  directions  at  the  same  time.  The  table  moves  so 
that  the  part  location  for  the  currently  placed  part  is 
aligned  with  the  placement  position  of  the  head  carrying  the 
part.  Thus  the  appropriate  metric  to  compute  distances 

acceleration/decceleration  effects)  is  the  Tchebyshev  or 
loo  -metric,  denoted  by  d^.  The  carrier  consists  of  slots 

containing  reels  with  individual  parts.  Each  reel  contains 
between  part  locations  on  the  PC  board  ( ignoring  one  type  of 
part,  and  several  slots  may  be  loaded  with  the  same  type  of 
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z-direction 

carrier 


I 


pick-up  position 


turret 


placement  position 


y-direction 


x-direction 


table 


Figure  2.1.  Schematic  top  view  of  the  Panasonic  Mkl  machine 

reel.  The  carrier  moves  back  and  forth  to  align  parts  to  be 
placed  with  the  pick  up  position  of  the  robot  head.  The 
movements  of  the  three  subsystems  are  synchronized  by  a cam 
shaft.  The  cam  shaft  has  to  rotate  once  between  the 
placement  of  two  successive  parts.  The  cam  profiles 
determine  a "time  window"  during  which  the  table  must  reach 
the  next  placement  location.  If  the  table  cannot  reach  its 
destination  during  this  time  window,  a brake  is  applied  to 
the  cam  shaft,  stopping  it  just  before  the  Mkl  locks  the 
table  into  position.  This  is  referred  to  as  a late  table. 
The  carrier  operates  similarly.  If  the  carrier  cannot  reach 
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its  destination  during  its  time  window,  a brake  is  again 
applied,  resulting  in  a late  carrier. 

2.2.2  Statement  of  the  Problem 

Given  an  a x b PC  board  with  N part  locations  P]^,  ..., 
Pn,  I wish  to  determine  a placement  sequence  o(l,...,N), 
that  minimizes  the  total  assembly  time  (TAT).  Assume  the 
times  are  scaled  so  that  the  table  moves  with  a speed  of  one 
unit  of  length  per  unit  of  time  in  each  direction.  Further 
assume  that  no  delays  are  experienced  due  to  the  late 
carrier.  For  any  given  sequence  the  TAT  can  be  computed  as 
the  sum  of  the  Tchebyshev  distances  between  consecutive  part 
locations . Let  r be  the  nominal  rotation  time  of  the  cam 
shaft,  or  cycle  time,  and  let  w denote  the  length  of  the 
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in  place  locked 


Figure  2.2.  Time  interval  for  placement  of  single  part 
when  table  is  on  time 

time  window  allowed  by  the  cam  shaft  for  the  table  to  reach 
its  next  location.  If  the  table  aligns  the  part  location 
with  the  placement  position  during  this  time  window  as 
illustrated  in  Figure  2.2.,  the  time  necessary  to  place  a 
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part  is  r.  This  means  that  a lower  bound  on  the  TAT  is  N r. 
If  the  table  is  late,  the  time  necessary  to  place  a part  i+1 
after  part  i is 


r + (dT(Pa( , Pa(i+l))  - w) . 

I refer  to  the  quantity  dT(Pa^j,  Pa(i+i))  - w as  the 
delay  time.  I can  express  the  delay  time  in  general  as 

t i - max { 0 , dij(Pa^j,  Pa(i+1))  - w}, 

and  the  time  to  place  part  i+1  after  part  i as  r + tj_.  The 
TAT  then  can  be  described  as 

TAT  - 2a (i/...fN)  tr  + t-i ] 

- N r + 2a(l,...,N)  max{0,  dip(Pa^-Lj,  Pa  (i+1))  - w}  • 

The  time  interval  needed  for  placement  of  single  type  in  the 
event  that  the  table  is  late  is  illustrated  in  Figure  2.3. 
Since  N r is  constant,  the  problem  to  minimize  TAT  can  be 
reformulated  as 

minimize  ^ # ^ ^ max{0,  drji(Pa^j,  Pa(i+l))  - w}* 

Note  that  a lower  bound  on  the  last  formulation  is 
zero.  This  lower  bound  can  be  achieved  if  every  part  in  the 
placement  sequence  is  placed  within  the  time  window  w. 
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Figure  2.3.  Time  interval  for  placement  of  single  part  when 
table  is  late 


From  the  preceding  analysis  it  is  obvious  that  an 
economical  placement  sequence  and  allocation  of  parts  to  the 
parts  carrier  can  decrease  the  TAT  necessary  to  populate  a 
PC  board,  and  thus  increase  productivity.  A 'good' 
allocation  of  parts  to  the  carrier  will  reduce  the  delay 
caused  by  a late  carrier.  Similarly,  a 'good'  placement 
sequence  will  reduce  the  delay  caused  by  a late  table. 
However,  there  is  a significant  interaction  between  these 
problems.  Given  an  allocation  of  parts  on  the  parts  carrier, 
a placement  sequence  that  does  not  delay  the  carrier  may 
cause  the  table  to  be  late.  If  a placement  sequence  is 
specified,  an  allocation  that  does  not  delay  the  table  may 
delay  the  carrier.  Thus  to  optimize  the  assembly  process  one 
has  to  solve  a complex  problem  comprised  of  two  notorious 
combinatorial  optimization  subproblems,  the  traveling 
salesman  problem  and  a nonlinear  assignment  problem.  For 
possible  mathematical  formulations  of  such  problems  see  Hu 
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and  Carter  [1989]  or  Leipala  and  Nevalainen  [1989].  To 
manage  this  difficult  problem  I consider  a "divide  and 
conquer"  strategy. 

To  reduce  the  dependence  of  the  two  subproblems,  let 
the  parts  be  allocated  to  the  reels  on  the  parts  carrier  so 
that  if  one  part  occupies  more  than  one  reel,  then  these 
reels  are  in  adjacent  slots  on  the  carrier.  Thus  the  carrier 
slots  are  assigned  part  types,  say  part  type  I,  II,  etc.,  in 
a consecutive  fashion.  This  amounts  to  determining  the 
assignment  of  parts  to  the  carrier,  and  keeps  carrier 
movement  small. 

The  parts  are  then  placed  on  the  PC  board  according  to 
part  types.  Starting  with  the  part  type  I occupying  a slot 
at  one  end  of  the  carrier,  all  the  parts  of  this  type  are 
placed  first.  The  carrier  then  aligns  the  slot  containing 
parts  of  type  II  (this  slot  is  adjacent  to  the  slot 
containing  parts  of  type  I)  with  the  pickup  position,  and 
all  of  the  parts  of  type  II  are  placed  on  the  board.  This 
process  continues  one  part  type  at  a time  until  all  parts 
required  are  placed  on  the  board.  This  way,  the  carrier 
moves  only  from  one  reel  slot  to  another  adjacent  slot 
either  when  a reel  is  depleted  or  when  all  parts  of  one  type 
are  placed.  The  overall  PC  board  assembly  is  then  determined 
by 

the  order  in  which  parts  of  each  type  are  placed, 
the  order  in  which  the  different  types  of  parts  are 
placed. 


An  added  benefit  of  this  strategy  is  minimum  carrier 
movement  and  simplicity. 


21 


The  PC  board  assembly  problem  can  now  be  modelled  as  a 
"clustered"  traveling  salesman  problem  (TSP)  described  in 
Chisman  [1975].  This  problem  can  be  stated  as  follows.  In  an 
N city  traveling  salesman  problem,  the  salesman  must  not 
only  visit  each  city  once  and  only  once,  but  a subset 
(cluster)  of  these  cities  must  be  visited  contiguously. 

There  may  be  more  than  one  cluster  within  a given  problem. 

In  the  PC  board  assembly  context  the  part  types  correspond 
to  the  clusters . 

The  clustered  TSP  formulation  of  the  PC  board  assembly 
problem  can  be  solved  suboptimally  with  the  following 
procedure: 

Step  1 : Find  the  optimal  placement  sequence  for  each  part 
type. 

Step  2 : Given  the  optimal  placement  sequence  for  each  part 
type  determine  the  sequence  in  which  the  part  types 
should  be  placed  on  the  board. 

At  Step  1,  the  placement  sequence  problem  for  single 
part  type  can  be  formulated  as  a Hamiltonian  cycle  problem 
in  a "Tchebyshev"  grid  graph.  This  problem  is  shown  to  be  NP 
complete  in  Chapter  3.  However,  in  a special  case  this 
problem  can  be  solved  in  polynomial  time.  This  result, 
together  with  a polynomial  time  algorithm  that  finds  the 
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Hamiltonian  cycle,  if  it  exists,  is  presented  in  Chapter  4. 
This  solution  seems  to  be  adequate  in  the  practical  sense. 
Information  obtained  in  Step  1 is  used  in  Step  2 to 
determine  the  sequence  in  which  the  part  types  should  be 
placed  on  the  board.  A sequence  mending  heuristic  is 
presented  in  Chapter  5 for  this  purpose. 


CHAPTER  3 

HAMILTONIAN  PROBLEMS  IN  GRID  GRAPHS 


In  Chapter  2,  I decomposed  the  PC  board  assembly 
problem  into  two  subproblems . One  of  these  was  a placement 
sequence  problem  for  individual  part  types.  In  this  chapter, 
I further  analyze  the  placement  sequence  problem  for  a 
single  part  type.  I formulate  the  problem  as  a Hamiltonian 
cycle  problem  in  "Tchebyshev"  grid  graphs  and  show  that  this 
problem  is  NP-complete. 


3 . 1 Terminology 

In  the  following  I introduce  some  graph  theory  used  to 
further  describe  the  problem.  I follow  closely  the 
definitions  in  Harary  [1969]  and  Itai,  Papadimitriou  and 
Szwarcfiter  [1982]. 

Let  G(V,  E)  denote  an  arbitrary  graph  with  vertex  set  V 
and  edge  set  E.  A subgraph  of  G(V,  E)  is  a graph  having  all 
vertices  and  edges  in  G(V,  E).  A spanning  subgraph  is  a 
subgraph  containing  all  vertices  of  G(V,  E).  For  any  set  S 
of  vertices  in  V,  the  induced  subgraph  <S>  is  the  edge- 
maximal  subgraph  of  G(V,  E)  with  vertex  set  S.  Thus  any  two 
vertices  are  adjacent  in  <S>  if  and  only  if  they  are 
adjacent  in  G(V,  E).  A walk  of  a graph  G(V,  E)  of  length  n 
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is  an  alternating  sequence  of  vertices  and  edges  vq,  ej_,  V]_, 
e2,...,  vn_^,  en/  vn,  beginning  and  ending  with  vertices,  in 
which  each  edge  is  incident  with  the  two  vertices 
immediately  preceding  and  following  it.  A walk  may  also  be 
denoted  simply  by  the  sequence  of  vertices  vq,...,  vn  (the 
edges  being  evident  by  the  context).  A walk  is  closed  if  vq 
= vn.  It  is  a trail  if  edges  are  distinct,  and  a path  if  the 
vertices  are  distinct.  A closed  path  with  more  than  two 
vertices  is  a cycle . A Hamiltonian  cycle  is  a spanning 
cycle.  A graph  is  said  to  be  Hamiltonian  if  it  contains  a 
Hamiltonian  cycle. 

A graph  G is  connected  if  for  every  pair  {u,  v}  of 
distinct  vertices  there  is  a path  between  u and  v.  A 
(connected)  component  of  a graph  is  a maximal  connected 
subgraph.  A cutvertex  is  a vertex  whose  deletion  increases 
the  number  of  components.  G is  2-connected  if  G is  connected 
and  has  no  cutvertex,  or,  equivalently,  G is  2-connected  if 
for  every  pair  {u,  v}  of  distinct  vertices  there  are  two 
internally-dis joint  paths  between  u and  v. 

Let  G00  be  the  infinite  graph  whose  vertex  set  consist 
of  all  the  points  of  the  plane  with  integer  coordinates.  Two 
vertices  of  G00  are  adjacent  if  and  only  if  the  Euclidean 
distance  between  them  is  1 . A rectilinear  grid  graph  is  a 
finite  induced  subgraph  of  G00.  Let  vx  and  vy  be  the 
coordinates  of  the  vertex  v.  Let  R(m,  n)  be  the  grid  graph 
whose  vertex  set  is  V = {v:  1 <;  vx  <;  m,  1 s vy  ^ n} . A 
rectangular  rectilinear  graph  with  dimensions  m x n is  a 
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rectilinear  grid  graph  isomorphic  to  R(m,  n)  for  some  m and 
n-  A rectilinear  (m,  n)-grid  graph,  shown  in  Figure  3.1,  is 
a rectilinear  grid  graph  induced  by  a rectangular 
rectilinear  graph  with  dimensions  m x n. 


Figure  3.1.  A rectilinear  (4,  5) -grid  graph 

The  Tchebyshev  versions  of  these  graphs  are  defined  in 
the  same  way  using  a Tchebyshev  distance.  Two  vertices  of 
G00  are  adiacent  if  and  only  if  the  Tchebyshev  distance 
between  them  is  equal  to  1 . For  vertices  u and  v with 
coordinates  (ux,  Uy)  and  (vx,  Vy),  respectively,  the 
Tchebyshev  distance  between  u and  v is  give  by 

max  { | ux  - vx | , | uy  - vy | } . 

Thus  a Tchebyshev  grid  graph  is  defined  as  a finite 
induced  subgraph  of  G00.  Let  vx  and  Vy  be  the  coordinates 
of  the  vertex  v.  Let  T(m,  n)  be  the  Tchebyshev  grid  graph 
whose  vertex  set  is  V = {v:  1 <;  vx  <;  m,  1 <;  vy  s n}.  A 

rectangular  Tchebyshev  graph  with  dimensions  m x n is  a 
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Tchebyshev  grid  graph  isomorphic  to  T(m,  n)  for  some  m and 
n.  A Tchebyshev  (m,  n)-arid  graph . illustrated  in  Figure  3.2, 
is  a Tchebyshev  grid  graph  induced  by  a rectangular 
Tchebyshev  graph  with  dimensions  m x n. 


Figure  3.2.  A Tchebyshev  (4,  5) -grid  graph 

3.2  Placement  Sequence  Problem  for  Single  Part  Type 

Consider  the  placement  of  N parts  of  a single  type  at 
locations  P^,  ...,  PN  on  an  a x b PC  board.  Assuming  the 
times  are  scaled  so  that  the  table  moves  with  a speed  of  one 
unit  of  length  per  unit  of  time,  subdivide  the  PC  board  into 
w/2  x w/2  cells  (like  a checker  board),  where  w is  the  time 
window  available  for  placing  a part.  Assume  2 a and  2 b 
(twice  the  board  dimensions)  are  divisible  by  w to  obtain  an 
even  subdivision  of  the  PC  board.  Figure  3.3  shows  an 
example  of  such  a subdivision. 
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a)  b) 


c) 


Figure  3 . 3 . Modelling  of  the  PC  board  assembly  problem  on 

Panasonic  Mkl 

(a)  PC  board  with  part  locations 

(b)  partitioned  PC  board 

(c)  occupancy  graph  (Tchebyshev  grid  graph) 

Let  q be  the  number  of  cells  obtained  in  this  fashion. 
The  cells  are  said  to  be  neighbors  if  they  share  a common 
side  or  common  corner.  A cell  is  said  to  be  unoccupied  if  it 
does  not  contain  any  placement  location  Pj_  and  occupied  if 
it  does.  It  follows  that  the  time  required  to  move  between 
any  two  points  in  neighboring  cells  is  at  most  w.  Thus  if 
consecutive  parts  in  the  placement  sequence  are  located  in 
neighboring  cells  then  TAT  = N r.  Therefore,  a sufficient 
condition  for  an  optimal  solution  is  that  there  exists  a 
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sequence  with  the  property  that  consecutively  placed  parts 
lie  in  the  same  or  neighboring  cells. 

We  can  construct  a rectangular  graph  G associated  with 
the  partitioned  PC  board  as  follows.  With  every  cell 
associate  the  corresponding  vertex  of  G.  Two  vertices  are 
adjacent  in  G if  and  only  if  the  represented  cells  are 
neighbors.  Thus  each  edge  in  G represents  a possible 
movement  of  the  table  shorter  than  the  time  window  w.  Let  T 
be  a subgraph  of  the  rectangular  graph  G induced  by  the 
vertices  associated  with  nonempty  cells  as  shown  in  Figure 
3.3.  The  graph  r represents  the  nonempty  cells  of  the 
partitioned  PC  board  with  edges  indicating  possible 
movements  of  the  table  that  can  be  executed  in  time  less 
than  w.  I will  refer  to  such  a graph  as  an  occupancy  graph. 
Note  that  G is  a Tchebyshev  rectangular  graph,  and  the 
occupancy  graph  T is  a Tchebyshev  grid  graph.  The  placement 
sequence  problem  can  now  be  restated  in  graph  theoretical 
terms.  In  particular,  is  the  occupancy  graph  T Hamiltonian 
or,  equivalently,  is  there  a placement  sequence  with 
consecutive  parts  located  in  neighboring  cells? 

3.3  NP-Completness  of  the  Placement  Sequence  Problem 

In  this  section  I will  show  that  the  Hamiltonian  path 
and  Hamiltonian  cycle  problems  for  general  Tchebyshev  grid 
graphs  are  NP-complete. 
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Theorem  3.1;  The  Hamiltonian  cycle  problem  for  general 
Tchebyshev  grid  graphs  is  NP-complete. 

Proof ; I will  show  that  the  Hamilton  cycle  problem  for  grid 
graphs  polynomially  transforms  to  the  Hamiltonian  cycle 
problem  for  Tchebyshev  grid  graphs. 

Since  a nondeterministic  algorithm  needs  to  only  guess 
an  ordering  of  vertices  and  check  in  polynomial  time  that 
all  the  required  edges  belong  to  the  edge  set  of  the  given 
graph,  HC  e NP. 

Suppose  G = (V,  E)  is  a grid  graph.  The  vertex  set  of  a 
corresponding  Tchebyshev  grid  graph  T = (V',  E)  is 
determined  by 


f;  (vx,  vy) 


“Mvx,  vy) 


1 -1 

1 1 


and  two  vertices  in  r are  adjacent  if  and  only  if  the 
corresponding  vertices  are  adjacent  in  G.  The  transformation 
f can  be  computed  by  a polynomial  time  algorithm.  An  example 
of  the  transformation  f is  shown  in  Figure  3.4. 

It  is  also  easy  to  see  that  G contains  a Hamiltonian 
cycle  if  and  only  if  r does.  The  Tchebyshev  grid  graph  T is 
obtained  from  the  grid  graph  G by  rotating  G by  45  degrees 
and  scaling  it  by  a factor  of  square  root  of  2.  This  means 
that  G and  r are  isomorphic.  Thus  HC  on  a grid  graph 
polynomially  transforms  to  HC  on  a Tchebyshev  grid  graph.  ■ 
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Figure  3.4.  Transformation  f 

Corollary  3.1;  The  Hamiltonian  path  problem  for  a general 
Tchebyshev  grid  graphs  is  NP-complete. 

Proof:  The  Hamiltonian  path  problem  for  a general  grid  graph 
is  NP-complete  (Itai,  Papadimitriou,  and  Szwarcfiter 
[1982]).  ■ 

In  this  chapter,  I have  modelled  the  placement  sequence 
problem  as  a Hamiltonian  cycle  problem  in  a graph  with  grid 
structure,  or  grid  graphs.  Also,  I have  shown  this  problem 
belongs  to  the  class  of  NP-complete  problems.  In  the  next 
chapter,  I will  study  a special  case  of  grid  graphs, 
Tchebyshev  (4,  n)-grid  graphs,  and  describe  the  necessary 
and  sufficient  conditions  for  the  existence  of  the 
Hamiltonian  cycle  in  such  graphs.  Based  on  this  result  I 
will  formulate  a heuristic  algorithm  to  find  a placement 
sequence  in  a Tchebyshev  (4,  n)-grid  graph. 


CHAPTER  4 

HAMILTONIAN  GRID  GRAPHS 


The  literature  on  Hamiltonian  cycles  in  ordinary  grid 
graphs  is  very  sparse;  I have  not  found  any  literature  on 
the  Hamiltonian  cycle  problem  in  a Tchebyshev  grid  graph. 
Luccio  and  Mugnai  [1978]  provide  conditions  for  the 
existence  of  Hamiltonian  cycles  and  paths  in  a rectangular 
graph.  Itai,  Papadimitriou,  and  Szwarcfiter  [1982]  show  that 
the  Hamiltonian  cycle  problem  in  a rectilinear  grid  graph  is 
NP-complete.  I have  shown  in  Chapter  3 that  the  Hamiltonian 
cycle  problem  in  a general  Tchebyshev  grid  graph  is  NP- 
complete.  However,  in  a special  case,  when  the  Tchebyshev 
grid  graph  is  induced  by  a 4 x n Tchebyshev  rectangular 
graph,  this  problem  can  be  solved  in  polynomial  time.  For 
this  special  case,  I will  present  necessary  and  sufficient 
conditions  for  the  existence  of  a Hamiltonian  cycle  in  a 
Tchebyshev  grid  graph,  together  with  a polynomial-time 
algorithm.  I will  also  present  a polynomial-time  algorithm 
to  find  a Hamiltonian  cycle  in  a grid  graph  induced  by  a 
rectangular  rectilinear  graph  of  dimensions  4 x n. 


31 


32 


4 . 1 Tchebvshev  Grid  Graphs  and  Hamiltonian  Cycles 

In  Chapter  3,  I modelled  the  placement  sequence  problem 
for  a single  part  type  as  a Hamiltonian  cycle  problem  on  an 
occupancy  graph  (Tchebyshev  grid  graph).  The  PC  board  can  be 
divided  into  cells.  If  a cell  is  occupied,  it  is  represented 
by  a vertex  of  the  occupancy  graph.  Two  vertices  of  the 
occupancy  graph  are  adjacent  if  the  corresponding  occupied 
cells  are  neighbors.  A Hamiltonian  cycle  in  the  occupancy 
graph  provides  a placement  sequence  with  minimal  assembly 
time . 

In  this  section  I study  the  properties  of  Tchebyshev 
(4,  n)-grid  graphs.  I develop  necessary  and  sufficient 
conditions  for  such  graphs  to  be  Hamiltonian  and  present  a 
polynomial-time  algorithm  to  find  a Hamiltonian  cycle  if  it 
exists . 

4.1.1  Hamiltonian  Cycles  in  Tchebyshev  ( 4 . n)-Grid  Graphs 

The  rectilinear  grid  graphs  and  Tchebyshev  grid  graphs 
differ  only  in  their  respective  edge  sets.  The  Tchebyshev 
grid  graph,  as  opposed  to  the  rectilinear  grid  graph, 
contains  diagonal  edges  between  vertices . These  extra  edges 
make  a Tchebyshev  grid  graph  "more  likely"  to  contain  a 
Hamiltonian  cycle. 

Luccio  and  Mugnai  [1978]  show  that  a Hamiltonian  cycle 
in  a rectangular  rectilinear  graph  exists  if  and  only  if  the 
number  of  vertices  in  the  graph  is  even.  Thus  a rectangular 
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rectilinear  graph  with  dimensions  m x n is  Hamiltonian  as 
long  as  at  least  one  dimension  is  even.  A rectangular 
Tchebyshev  graph,  on  the  other  hand,  is  always  Hamiltonian. 

Proposition  4.1;  A Tchebyshev  rectangular  graph  with 
dimensions  m x n is  Hamiltonian. 

Proof ; Let  G be  a rectangular  Tchebyshev  graph  with 
dimensions  m x n.  If  either  m or  n is  even,  a Hamiltonian 
cycle  in  G can  be  constructed  the  same  way  as  described  by 
Luccio  and  Mugnai  [1978]  for  rectangular  rectilinear  graphs. 
An  example  of  a Hamiltonian  cycle  in  this  case  is  shown  in 
Figure  4.1.  below. 


Figure  4 . 1 . Hamiltonian  cycle  in  a m x n rectangular 
Tchebyshev  graph  with  either  m or  n even 

If  both  dimensions  of  G are  odd,  a path  can  be 
constructed  which  is  similar  to  the  path  constructed  in  the 
case  when  m or  n is  even.  Initially,  the  path  consists  of 
the  upper  left  corner  vertex  and  extends  in  the  "down-right- 
up  to  the  second  row-right  again"  pattern  until  it  reaches 
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column  (m  - 2).  Since  m is  odd,  the  path  is  extending 
downward  in  this  column.  Extend  the  path  through  the  bottom 
row  vertices  of  the  last  two  columns  and  continue  in  the 
following  fashion:  one  vertex  up  - one  vertex  left  - one 

vertex  up  - one  vertex  right  - (repeat)  until  the  path 

reaches  the  second  row  from  the  top.  Since  n is  odd,  the 
last  vertex  of  the  path  constructed  thus  far  lies  in  column 
(m  - 1).  Extend  the  path  to  the  upper  right  corner  vertex 
using  a diagonal  edge,  and  complete  the  cycle  by  extending 
the  path  through  the  first  row  of  vertices.  The  constructed 
Hamiltonian  cycle  is  shown  in  Figure  4.2.  below.  ■ 


Figure  4 . 2 . Hamiltonian  cycle  in  a m x n rectangular 
Tchebyshev  graph  with  both  m and  n odd 


As  I noted  earlier,  the  equivalent  proposition  is  not 
true  for  rectangular  rectilinear  graphs.  For  example,  a 
3x3  rectangular  rectilinear  graph  is  not  Hamiltonian. 
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4 » 1 » 2 Necessary  and  Sufficient  Conditions  for  a Hamiltonian 

Cycle  in  Tchebvshev  (4.  nl-Grid  Graphs 

In  this  section  I will  present  the  main  result,  namely 
necessary  and  sufficient  conditions  for  a Tchebyshev  (4,  n)- 
grid  graph  to  be  Hamiltonian.  As  I discussed  in  Chapter  3,  a 
Tchebyshev  (4,  n)-grid  graph  represents  the  nonempty  cells 
of  a PC  board  that  has  been  partitioned  into  four  rows  with 
n cells  in  each  row.  The  minimum  and  maximum  dimensions  of 
the  PC  board  mounting  area  allowed  by  the  Panasonic  Mkl  are 
90  x 52  and  508  x 373  millimeters,  respectively.  The  table 
can  move  at  three  different  speeds:  100,  200,  and  400 
millimeters  per  second.  Recall,  that  the  time  window  w is  a 
time  interval  during  which  the  table  can  move  to  the  next 
part  location  without  being  "late."  The  Tchebyshev  (4,  n)- 
grid  graph  is,  therefore,  large  enough  to  model  the  PC  board 
assembly  problem  if  the  length  of  the  time  window  w is  at 
least  one  second. 

I first  introduce  some  terminology.  A graph  is  said  to 
be  embedded  in  a surface  S when  it  is  drawn  on  S so  that  no 
two  edges  intersect.  A graph  is  planar  if  it  can  be  embedded 
in  the  plane;  a plane  graph  is  a drawing  of  a planar  graph 
in  the  plane. 

If  the  curves  and  points  corresponding,  respectively,  to 
the  edges  and  vertices  of  a plane  graph  are  deleted,  then 
the  remainder  of  the  plane  falls  into  connected  components, 
called  faces . The  unbounded  region  is  called  the  exterior 
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Figure  4.3.  Planar  graph  and  plane  graph 

face . A boundary  of  a face  is  the  set  of  edges  in  its 
closure  (see  Figure  4.4.). 


exterior  face 


Figure  4.4.  Plane  graph 

I have  shown  that  the  Hamiltonian  cycle  problem  in  a 
general  Tchebyshev  grid  graph  is  an  NP-complete  problem.  It 
is  common  for  NP-complete  problems  to  have  special  instances 
when  the  problem  becomes  polynomially  solvable.  The 
Hamiltonian  cycle  problem  is  no  exception.  I will  show  that 
2-connectedness  is  a sufficient  condition  for  a Tchebyshev 
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Hamiltonian.  Since  every  Hamiltonian  graph  is  2-connected, 
it  follows  that  2-connectedness  is  also  a necessary 
condition. 

The  following  Lemmas  will  be  used  to  show  this  result. 

Lemma  4.1:  The  boundary  of  any  face  of  a 2-connected  plane 
graph  is  a cycle. 

Proof:  Let  G be  a 2-connected  plane  graph  with  N vertices. 
Let  s be  the  boundary  of  a face  f of  G.  The  boundary  s 
defines  a walk  w = xj.,  X2,...,  x^  = that  can  be 
constructed  as  follows.  Let  x^  be  any  vertex  in  s and  e = 
(xl'  x2 ) be  an  edge  in  s incident  to  x^ . Assume  f lies  to 
the  "right"  of  e,  when  e is  considered  as  a directed  edge 
from  x^  to  X2 . The  walk  w is  then  extended  to  the  first  edge 
in  G incident  to  X2  found  in  the  "counterclockwise" 
direction  from  e.  Since  G is  2-connected  and,  thus,  the 
degree  of  every  vertex  in  G is  at  least  two,  the  walk  w can 
be  extended  until  it  reaches  a vertex  already  included  in 
the  walk  w,  say  x-[  (Figure  4.5.). 


Figure  4.5.  Walk  w 
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Suppose  x-l  is  the  first  repeated  vertex  in  the  walk  w, 

i.e.  w = X]_,  — , xi, , xj  = xi  with  w'  = x^, , xj  = x-^. 

Then  w'  is  closed  and  contains  no  repeated  vertices,  hence 
w'  is  a cycle  (Figure  4.6.).  If  x\  = xj^,  then  s = w = w'  is 
a cycle,  and  there  is  nothing  to  prove. 

Suppose  x^  * x^.  Since  G can  be  embedded  on  the  sphere, 
a North  Pole  of  the  sphere  can  be  chosen  in  such  a way  that 
G is  stereographically  projected  onto  a plane  with  the 
vertices  of  w'  visited  in  the  counterclockwise  direction  and 
the  face  f on  the  outside  of  w ' (Figure  4.6.).  From  the 
construction  of  w it  follows  that  there  are  no  edges 
incident  to  w'  - x^  on  the  outside  of  w'  . This  implies  that 
all  paths  in  G joining  and  x-^+i  have  to  include  vertex 

xi • Thus  x-l  is  a cutvertex,  which  contradicts  the  assumption 
that  G is  2-connected.  ■ 


Figure  4.6.  Cycle  w' 

It  is  well-known  that  every  2-connected  planar  graph 
can  be  embedded  in  the  plane  so  that  any  specified  face  is 


the  exterior.  Thus  the  Lemma  4.1  implies  the  following 
corollary. 
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Corollary  4.1;  The  boundary  of  the  exterior  face  of  a 2- 
connected  plane  graph  is  a cycle. 


Lemma  4.2;  Every  vertex  in  the  exterior  face  of  a 2- 
connected  plane  graph  is  in  the  boundary  of  the  exterior 
face. 

Prooi:  Let  G be  a 2 -connected  plane  graph.  Let  fe  be  the 
exterior  face.  Let  v be  a vertex  in  fe  and  let  u be  an 
arbitrary  vertex  of  G.  Since  G is  2-connected,  there  are  two 
paths  between  v and  u with  no  vertices  in  common  except  v 
and  u.  Let  s be  the  cycle  formed  by  these  two  paths.  Since 
the  interior  of  the  cycle  s is  bounded  by  the  cycle,  s is 
the  boundary  of  a union  of  interior  faces  of  G.  Since  the 
vertex  v is  contained  both  in  fe  and  in  the  boundary  of  an 
interior  face,  it  must  lie  in  the  boundary  of  fe.  Since  v 
was  chosen  arbitrarily  in  the  exterior  face  fe,  every  vertex 
in  fe  must  lie  in  the  boundary  of  fe.  Thus  every  vertex  in 
the  exterior  face  of  a 2-connected  plane  graph  is  in  the 
boundary  of  the  exterior  face.  ■ 

As  a direct  consequence  of  Corollary  4.1  and  Lemma  4.2, 

I have  the  following  lemma. 


Lemma  4.3;  Every  2-connected  plane  graph  contains  a cycle 
with  no  vertex  interior  to  the  cycle  exterior. 
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A Tchebyshev  grid  graph  G = (V,  E)  is  not  planar,  in 
general.  To  be  able  to  use  the  preceding  result,  I will 
introduce  a subgraph  of  G,  a plane  graph  G~,  that  is  closely 
related  to  G and  retains  some  of  the  properties  of  G.  In 
particular,  G~  is  2-connected  if  and  only  if  G is  2- 
connected. 

The  nonplanarity  of  a Tchebyshev  grid  graph  is  caused 
by  the  diagonal  edges  between  any  four  mutually  adjacent 
vertices  arranged  in  a square.  Such  a set  of  four  mutually 
adjacent  vertices  arranged  in  a square  is  referred  to  as  a 
4-cluster . In  other  words,  a 4 -cluster  is  a clique  (maximal 
complete  subgraph)  with  four  vertices.  The  plane  graph 
G~  = (V,  E ’ ) is  obtained  by  removing  one  diagonal  edge  in 
every  4-cluster  contained  in  G.  The  next  proposition  shows 
that  G~  is  2-connected  if  and  only  if  G is. 

Proposition  4.2:  The  graph  G~  obtained  from  a Tchebyshev 
grid  graph  G by  removing  one  diagonal  edge  in  every  4- 
cluster  is  2-connected  if  and  only  if  G is. 

Proof ; Assume  G is  2-connected.  Let  x and  y be  two  vertices 
in  different  connected  components  of  G~\{v}.  Note  x and  y 
are  vertices  of  G.  Since  G is  2-connected,  there  exist  two 
internally-dis joint  paths  P and  Q in  G joining  the  vertices 
x and  y.  At  least  one  of  these  paths,  say  P,  contains  at 
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least  one  diagonal  edge  in  some  4-cluster  {a,  b,  c,  d},  say 
e = (b,  c),  that  is  not  contained  in  the  edge  set  of  G-. 
Otherwise,  P and  Q are  disjoint  paths  in  both  G and  G-,  and 
there  is  nothing  to  prove.  A pair  of  paths  P'  and  Q'  joining 
x to  y that  do  not  contain  e can  be  constructed  from  P and  Q 
as  follows.  If  the  vertex  a (or  d)  does  not  belong  to  the 
path  Q,  then  a path  P'  from  x to  y can  be  constructed  from  P 
by  replacing  the  diagonal  edge  e by  a pair  of  edges  (b,  a) 
and  (a,  c)  (or  (b,  d)  and  (d,  c)).  In  this  case  path  Q'  is 
the  same  as  Q.  if  both  vertices  a and  d belong  to  Q,  let  P^ 
and  P2  be  the  subpaths  of  P from  x to  b and  from  c to  y, 
respectively,  and  let  Qi  and  Q2  be  the  subpaths  of  Q from  x 
to  a and  from  d to  y,  respectively.  Then  P'  = Pi,  b,  d,  Q2 
and  Q'  = Qi,  a,  c,  P2  are  two  internally-dis joint  paths  from 
x to  y that  do  not  contain  the  edge  e.  These  constructions 
are  illustrated  in  Figure  4.7. 

If  either  P'  or  Q ' contains  other  edges  that  are  not 
included  in  the  edge  set  of  G~,  the  paths  can  be  adjusted  in 
the  same  fashion  until  all  such  edges  are  eliminated  from 
the  two  paths.  Let  P'  and  Q'  be  the  resulting  internally- 
dis  joint  paths  between  x and  y.  Since  P'  and  Q'  contain  no 
diagonal  edges  between  vertices  of  any  4-cluster  contained 
in  the  vertex  set  V,  p*  and  Q'  are  in  G~.  Thus  G-  is  2- 
connected. 

Conversely,  suppose  G~  is  2-connected.  Note  that  G-  is 
a subgraph  of  G with  the  same  vertex  set.  Let  u and  v be 
vertices  in  G.  Since  u and  v are  vertices  in  G~,  there  exist 
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Figure  4.7.  Alternate  paths  P'  and  Q'  in  G~ 

two  internally-dis joint  paths  P and  Q joining  u and  v.  Since 
the  edge  set  of  G contains  the  edge  set  of  G~,  P and  Q are 
internally-dis joint  paths  in  G joining  u and  v.  Thus,  if  G~ 
is  2 -connected,  then  G is  2 -connected.  ■ 


I am  now  ready  to  state  and  prove  the  main  result. 

Proposition  4.3:  Every  2-connected  Tchebyshev  (4,  n)-grid 
graph  is  Hamiltonian. 

Proof ; Let  G(V,  E)  be  a 2-connected  Tchebyshev  (4,  n)-grid 
graph  and  let  G-  be  the  plane  graph  associated  with  G.  Let  S 
be  a set  of  cycles  in  G with  the  property  that  every  vertex 
in  each  such  cycle  is  on  the  cycle  or  in  its  interior.  By 


43 


the  preceding  Lemma  4.3  at  least  one  such  cycle  exists  in 
G~,  and,  since  G-  is  a subgraph  of  G with  the  same  vertex 
set,  the  set  S is  not  empty. 

Consider  the  Tchebyshev  (4,  n)-grid  graph  G.  I claim 
that  the  largest  (the  one  with  most  vertices)  cycle  C in  the 
set  S contains  all  vertices  in  the  vertex  set  V.  Suppose 
there  exists  a vertex  v that  is  not  contained  in  C.  Since 
there  are  no  vertices  of  the  vertex  set  V in  the  exterior  of 
C,  v is  in  the  interior  of  C.  The  vertex  v cannot  lie  in  the 
first  or  the  fourth  row  of  G,  since  these  vertices  belong  to 
every  cycle  in  S.  Assume  the  vertex  v lies  in  the  second 
row.  The  case  when  v lies  in  the  third  row  is  symmetric. 
Since  v is  in  the  cycle  interior,  vertex  u directly  above  v 
must  be  in  the  cycle  C.  Since  all  the  vertices  adjacent  to  u 
in  C are  also  adjacent  to  v in  G,  the  cycle  C = (...,  x,  u, 
y,...)  can  be  extended  to  a larger  cycle  C'  = (...,  x,  v,  u, 
y,...)  with  all  the  vertices  of  V on  it  or  in  its  interior. 
But  this  is  a contradiction,  since  C is  the  largest  cycle 
with  all  the  vertices  on  it  or  its  interior.  Hence  C 
contains  all  vertices  in  the  vertex  set  V,  and  the 
Tchebyshev  grid  graph  G is  Hamiltonian.  ■ 

The  proof  of  Proposition  4.3  implies  the  following 
Corollary  4.2. 

Corollary  4.2;  Let  G be  a Tchebyshev  (4,  n)-grid  graph,  and 
let  C be  a cycle  in  G such  that  every  vertex  of  G is 
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contained  either  in  the  cycle  C or  in  its  interior.  Then 
every  vertex  in  the  interior  of  C is  adjacent  to  two 
vertices  that  are  adjacent  in  C. 

Proof : Let  v be  a vertex  in  the  interior  of  C.  Then  v lies 
in  the  second  or  third  row  of  G.  Assume  the  v lies  in  the 
second  row.  The  case  when  v lies  in  th  third  row  is 
symmetric.  Vertex  u directly  above  v in  the  first  row  must 
be  in  C.  Vertex  u is  adjacent  to  two  vertices  in  C.  Since 
all  the  vertices  adjacent  to  u in  C are  also  adjacent  to  v 
in  G,  it  follows  that  v is  is  adjacent  to  two  vertices  that 
are  adjacent  in  C,  namely  u and  another  vertex  adjacent  to  u 
in  C.  ■ 


Figure  4.8  A nonhamiltonian  Tchebyshev  (5,  5) -grid  graph 

The  Tchebyshev  (5,  5) -grid  graph  in  Figure  4.8  is  a 2- 
connected  nonhamiltonian  graph.  This  graph  illustrates  that 
2-connectedness  fails  as  a sufficient  condition  for 
Tchebyshev  (m,  n)-grid  graphs  with  m ss  5 to  be  Hamiltonian. 
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4.1.3  Hamiltonian  Cycle  Algorithm 

The  proof  of  Proposition  4.3  constitutes  an  algorithm. 
This  is  essentially  the  "convex  hull  insertion  procedure" 
due  to  Stewart  [1977],  and  described  in  Lawler,  Lenstra, 
Rinnooy  Kan,  and  Shmoys  [1985],  where  the  convex  hull  is 
replaced  by  an  envelope.  An  envelope  in  this  context  means  a 
cycle  in  G,  such  that  all  vertices  of  G are  either  in  the 
cycle  or  in  its  interior. 

Let  G ( V,  E)  be  the  Tchebyshev  (4,  n)-grid  graph.  The 
Hamiltonian  cycle  algorithm  for  G can  be  summarized  as 
follows : 

Step  1.  Construct  an  envelope  C of  G. 

Step  2.  If  C = V,  Stop.  Otherwise  go  to  Step  3. 

Step  3.  Let  x be  a vertex  in  V \ C.  Corollary  4.2  implies 
that  x must  be  adjacent  to  two  vertices,  say  u and 
v,  that  are  adjacent  in  C.  Thus  x can  be  inserted 
in  C between  the  vertices  u and  v.  Go  to  Step  2. 

For  the  implementation  of  the  algorithm  it  is 
convenient  to  adopt  the  following  convention.  Let  T be  the 
smallest  4 x n rectangular  Tchebyshev  graph  containing  G. 
Note  that  n is  the  number  of  columns  of  G.  Label  the 
vertices  of  T from  1 to  4 n starting  in  the  "northwest" 
corner  of  T from  left  to  right  one  row  at  a time.  The 
vertices  of  G are  then  indexed  using  the  vertex  labels  of  T. 
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A straightforward  implementation  of  this  algorithm 
starts  with  the  construction  of  a list  C that  contains  the 
vertex  indices.  Initially,  x is  the  lowest  indexed  vertex  of 
G and  the  adjacent  vertices  are  scanned  starting  from  the 
"north"  of  x.  If  x is  the  last  vertex  of  C,  the  vertices  of 
T adjacent  to  x are  scanned  in  a clockwise  direction 
starting  at  the  arc  joining  x and  the  predecessor  of  x in 
the  list  until  a vertex  in  G is  found.  This  vertex  is  added 
to  the  end  of  the  list  C.  Repeat  this  procedure  until  the 
vertex  at  the  beginning  of  the  list  is  added  to  the  list  for 
the  second  time.  At  the  end  of  this  "loop"  the  list  C 
contains  the  indices  of  the  envelope  vertices.  During  the 
execution  of  this  loop,  the  list  C is  updated  by  adding  a 
vertex  of  G to  the  end  of  the  list.  For  each  update  at  most 

seven  vertices  are  scanned  if  they  belong  to  G.  Since  the 

envelope  can  contain  up  to  four  vertices  in  some  columns,  C 
can  be  updated  at  most  4 n times.  Thus,  Step  1 can  be 
executed  in  0(n)  running  time. 

Now,  for  each  index  i of  a vertex  of  G in  the  second 
row  of  T not  included  in  C,  insert  i in  the  list  C between 
vertex  (i  - n)  and  its  predecessor.  For  each  index  i of  a 

vertex  of  G in  the  third  row  of  T not  included  in  C,  insert 

i into  the  list  C between  vertex  (i  + n)  and  its 
predecessor.  Recall  that  the  first  vertex  in  the  list 
appears  in  the  list  twice,  hence  every  vertex  in  the  list 
has  a predecessor.  Since  there  are  at  most  2 (n  - 2) 
vertices  of  G in  the  second  and  third  row  of  T that  are  not 
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in  C at  the  end  of  Step  1,  and  each  insertion  into  the  list 
can  be  done  in  constant-time  if  an  array  is  kept  indicating 
the  location  of  each  vertex  of  G in  C,  Step  3 can  be 
executed  in  0(n)  running  time. 

Since  each  sequential  step  of  the  algorithm  requires 
0(n)  running  time,  the  complexity  of  the  algorithm  is  0(n), 
where  n is  the  number  of  columns  in  G. 

4.2  Hamiltonian  Cycles  in  Rectilinear  (4,  n)-Grid  Graphs 

In  the  Panasonic  Mkl  machine  application  the  distances 
between  part  locations  on  the  PC  board  are  computed  using 
the  Tchebyshev  metric.  I shall  now  consider  the  same 
application  where  the  distances  are  computed  using  the 
rectilinear  or  L]^  metric.  The  rectilinear  metric  is 
applicable  to  those  machines  which  place  components  onto  the 
PC  board  in  a rectilinear  motion. 

The  PC  board  assembly  problem  analysis  in  Chapters  2 
and  3 is  independent  of  the  metric  used.  Hence  the  placement 
sequence  problem  for  a single  part  type  using  a rectilinear 
metric  can  be  formulated  as  a Hamiltonian  cycle  problem  in  a 
rectilinear  grid  graph.  In  this  section  I develop  a 
characterization  of  nonrectangular  rectilinear  (4,  n)-grid 
graphs  in  terms  of  the  "segments"  they  consist  of.  I will 
then  present  necessary  and  sufficient  conditions  for  a 
Hamiltonian  cycle  to  exist  in  this  special  case  of 
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rectilinear  grid  graphs.  The  results  presented  in  this 
section  are  proved  in  the  Appendix  A. 

Let  G be  a nonrectangular  rectilinear  (4,  n)-grid 
graph.  Assume  for  the  moment  that  G is  Hamiltonian,  and  let 
C be  a Hamiltonian  cycle  in  G.  Let  u and  v be  two  vertices 
in  the  first  and  last  columns  of  G,  respectively.  The 
vertices  u and  v separate  the  cycle  C into  two  disjoint 
paths  P and  Q.  The  graph  G can  now  be  decomposed  into 
segments"  consisting  of  sets  of  adjacent  columns  such  that 
both  P and  Q intersect  each  segment  in  a connected  path. 


o 

o 

• 

o 

Figure  4.9. Columns  with  ( 1 , 3 , 4 ) -pattern  and  ( 1 , 4 ) -pattern, 
respectively. 


Recall  that  G(V,  E)  is  a rectilinear  (4,  n)-grid  graph. 
A pattern  of  vertices  specifies  the  rows  of  a column 
containing  vertices  of  G.  For  example  in  Figure  4.9,  a 
column  with  ( 1 , 3 , 4 ) -pattern  contains  vertices  in  rows  1,  3, 
and  4,  and  a column  with  ( 1 , 4 ) -pattern  contains  vertices  in 
rows  1 and  4.  A segment  of  G is  the  largest  set  of  adjacent 
columns  with  the  same  pattern.  If  G consists  of  more  than 
one  segment,  the  first  and  last  segments  of  G are  referred 
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to  as  end  segments  , while  other  segments  are  referred  to  as 
interior  segments.  For  k = 2,  3,  and  4,  a k-seament  is  a 
segment  with  k vertices  in  each  column.  A segment  S is  said 
to  be  an  elementary  segment  if  one  of  the  following  is  true: 

i)  S is  a 2 -segment  with  one  column, 

ii)  S is  a 3-segment  with  one  or  two  columns, 

iii)  S is  a 4 -segment  with  one  to  five  columns. 

For  k = 2,  3,  and  4,  an  elementary  k-segment  is  an 
elementary  segment  with  k vertices  in  each  column.  Figure 
4.10  shows  all  nineteen  possible  elementary  segments. 


Figure  4.10.  Elementary  k-segments 

(a)  elementary  2 -segments 

(b)  elementary  3-segments 

(c)  elementary  4 -segments 
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Note  that  any  2-connected  rectilinear  (4,  n)-grid  graph 
G can  be  decomposed  into  segments,  since  it  consists  of 
columns  with  2,  3,  or  4 vertices.  Clearly,  if  G consists  of 
a single  segment,  then  G is  disconnected  provided  vertices 
of  G are  in  nonadjacent  rows,  or  it  is  a rectangular 
rectilinear  graph  with  vertices  in  adjacent  rows.  If  G is 
disconnected,  then  it  is  nonhamiltonian.  In  the  latter  case, 
Luccio  and  Mugnai  [1978]  show  that  a rectangular  rectilinear 
graph  is  Hamiltonian  if  and  only  if  the  number  of  vertices 
in  the  graph  is  even.  In  other  words,  a m x n rectangular 
rectilinear  graph  is  Hamiltonian  if  and  only  if  at  least  one 
dimension  is  even.  In  view  of  this  observation,  in  the 
remainder  of  this  section  G is  a 2-connected  nonrectangular 
rectilinear  (4,  n)-grid  graph  with  at  least  two  segments. 

The  next  two  lemmas  show  that  a Hamiltonian  cycle,  if 
it  exists,  intersects  any  segment  S of  G in  at  most  two 
disjoint  paths  that  cover  the  vertices  of  the  segment  S.  A 
path  or  paths  cover  vertices  of  a segment  S if  all  vertices 
of  the  segment  belong  to  the  path  or  paths . 


Lemma — 4 » 4 ; Let  G(V,  E)  be  a Hamiltonian  rectilinear  (4,  n)- 
grid  graph.  A Hamiltonian  cycle  intersects  any  end  segment  S 
of  G in  a single  path  that  covers  all  vertices  of  S. 


Lemma — 4.5;  Let  G(V,  E)  be  a Hamiltonian  rectilinear  (4,  n)- 
grid  graph.  A Hamiltonian  cycle  intersects  any  interior 


51 


segment  S of  G in  exactly  two  disjoint  paths  that  cover  all 
vertices  of  S. 

It  follows  from  the  discussion  above  that  any  2- 
connected  rectilinear  grid  graph  G can  be  divided  into 
segments.  Furthermore,  if  G is  Hamiltonian  and  consists  of 
more  than  one  segment,  then  the  Hamiltonian  cycle  intersects 
each  segment  in  at  most  two  disjoint  paths  that  cover  all 
vertices  in  a segment. 

Let  and  S2  be  two  adjacent  interior  segments  of  G, 
and  let  column  k be  the  end  column  of  Sj  adjacent  to  S2.  Let 
Pi  and  Qi  be  two  disjoint  paths  covering  the  vertices  of  Si 
with  endvertices  u and  v in  column  k,  respectively.  Let  P2 
and  Q2  be  two  disjoint  paths  covering  the  vertices  of  S2 
with  endvertices  x and  y in  column  (k  + 1),  respectively. 

The  paths  P-i  and  Q^,  i = 1,  2,  are  compatible  if  u and  x are 
in  the  same  row,  and  v and  y are  in  the  same  row.  In  other 
words,  the  paths  can  be  connected  by  "horizontal"  edges  of 
G,  so  that  the  vertices  of  the  union  of  the  two  segments  are 
again  covered  by  two  disjoint  paths.  We  say  such  segments 
are  path-compatible . Clearly,  for  G to  be  Hamiltonian  it 
must  be  a union  of  path-compatible  segments  with  "matching" 
end  segments . 

I shall  now  describe  pairs  of  disjoint  paths  that  cover 
vertices  in  the  various  interior  segments  of  G.  I will  show 
that  every  segment  can  be  replaced  by  an  elementary  segment 
to  construct  an  auxiliary  rectilinear  grid  graph  G^,  which 
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is  Hamiltonian  if  and  only  if  G is  Hamiltonian.  There  are 
three  different  kinds  of  segments  to  consider:  segments  with 
two,  three,  or  four  vertices  in  each  column.  Note  that  if  G 
contains  a segment  with  only  one  vertex  in  each  column,  then 
each  of  the  segment  vertices  is  a cutvertex,  and  G is  not 
Hamiltonian.  In  this  discussion  I assume  that  G is  2- 
connected,  and,  hence,  every  segment  in  G is  a k-segment 
with  k at  least  two. 

Given  an  interior  segment  S,  let  P and  Q denote  two 
disjoint  paths  covering  the  vertices  of  S with  one  endvertex 
of  both  P and  Q in  the  first  column  of  S and  the  other 
endvertex  of  P and  Q in  the  last  column  of  S.  The 
endvertices  of  P and  Q in  the  leftmost  column  of  S are 
referred  to  as  entering  vertices , and  the  endvertices  in  the 
rightmost  column  of  S are  referred  to  as  the  leaving 
vertices . If  u and  v are  entering  (or  leaving)  vertices  of  P 
and  Q in  column  p,  define  a function  fp,  that  assigns  a 
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value  to 
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of 

u and  v 
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o o o o 

Figure  4.11.  An  example  of  an  (i,  j) -bipath 

vertices  u and  v in  column  p,  and  leaving  vertices  x and  y 
in  column  q such  that  fp(u,  v)  = i and  fq(x,  y)  = j.  An 
example  of  (i,  j) -bipath  is  illustrated  in  Figure  4.11.  The 
ordered  pair  (i,  j)  representing  the  entering  and  leaving 
vertices  of  an  (i,  j) -bipath  is  referred  to  as  a bipath 
enter-leave  configuration.  Note  that  even  though  (i,  j)- 
bipaths  are  not  unique  in  some  segments,  it  is  the  bipath 
enter-leave  configuration  that  determines  path-compatibility 
of  segments.  For  a given  segment  S,  a segment  enter-leave 
pattern  is  the  set  of  all  bipath  enter-leave  configurations 
(if  j)  possible  in  S.  Two  segments  are  enter-leave 
equivalent  if  they  have  the  same  vertex  pattern  and  the  same 
enter-leave  pattern. 

The  following  lemmas  show  that  every  interior  segment 
is  enter-leave  equivalent  to  an  elementary  segment.  This 
means  that  segments  of  G can  be  replaced  by  smaller 
segments . 
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Figure  4.12.  A 2-segment  and  an  enter-leave  equivalent 
elementary  2-segment 


Lemma  4.6;  Any  interior  2 -segment  in  G is  enter-leave 
equivalent  to  an  interior  elementary  2-segment. 


The  following  Lemmas  4.7  and  4.8  give  similar  results 
for  interior  3-segments  and  4-segments.  These  results  are 
illustrated  in  Figure  4.13  and  4.14. 


Lemma  4.7;  Any  interior  3-segment  in  G is  enter-leave 
equivalent  to  an  interior  elementary  3-segment. 


o o o o o 

• <»»§• 


o o o o o o 


a)  b) 

Figure  4.13.  a)  An  (i,  j) -bipath  in  a 3-segment  with 

vertices  in  nonadjacent  rows 

b)  An  (i,  j) -bipath  in  a 3-segment  with 
vertices  in  adjacent  rows 
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Figure  4.14.  The  construction  of  a (1,  3) -bipath  in  a six 
column  4-segment  based  on  a (1,  3) -bipath  in 
a four  column  elementary  4 -segment 

Lemma  4.8;  Any  interior  4 -segment  in  G is  enter-leave 
equivalent  to  an  interior  elementary  4-segment. 

The  preceding  lemmas  and  Lemma  4 . 5 imply  that  any 
interior  segment  in  a 4 x n rectilinear  graph  G can  be 
replaced  by  an  enter-leave  equivalent  elementary  segment. 

In  a parallel  fashion,  I will  now  develop  the  same 
result  for  the  end  segments  of  G.  Lemma  4.4  shows  that  if  G 
is  Hamiltonian,  the  Hamiltonian  cycle  intersects  an  end 
segment  in  a single  path.  I will  refer  to  this  path  as  a i- 
path,  where  j represents  the  pattern  of  the  entering  or 
leaving  vertices  of  the  right  or  left  end  segment, 
respectively,  and  is  referred  to  as  a path  enter/leave 
configuration . For  a given  end  segment  S,  a segment 
enter /leave  pattern  is  the  set  of  all  path  enter/leave 
configurations  j possible  in  S.  Two  end  segments  are 
enter /leave  equivalent  if  they  have  the  same  vertex  pattern 
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and  the  same  enter/leave  pattern.  Note  that  "enter-leave" 
refers  to  interior  segments,  while  "enter/leave"  refers  to 
end  segments.  The  following  corollaries  show  that  all  end 
segments  are  path  equivalent  to  some  elementary  end  segment. 
In  particular,  an  end  segment  with  a two  vertex  pattern  can 
be  replaced  by  an  enter /leave  equivalent  end  segment  with 
one  column.  An  end  segment  with  a three  vertex  pattern  and 
an  odd  or  even  number  of  columns  can  be  replaced  by  an 
enter /leave  equivalent  end  segment  with  one  or  two  columns, 
respectively.  An  end  segment  with  a four  vertex  pattern  and 
with  more  than  two  columns  can  be  replaced  by  an  enter/leave 
equivalent  end  segment  with  two  columns.  These  results  are 
given  in  the  Corollaries  4.3  through  4.5. 

Corollary  4.3;  Any  end  segment  with  a two  vertex  pattern  is 
enter /leave  equivalent  to  an  elementary  end  segment. 

Corollary  4.4;  Any  end  segment  with  a three  vertex  pattern 
is  enter /leave  equivalent  to  an  elementary  end  segment  with 
one  or  two  columns. 

An  auxiliary  graph  G*  that  consists  only  of  elementary 
segments  can  now  be  constructed  for  any  rectilinear  (4,  n)- 
grid  graph  G.  Given  such  a rectilinear  grid  graph  G,  the  end 
segments  of  G are  replaced  by  enter /leave  equivalent 
elementary  end  segments  as  described  prior  to  Corollary  4.3. 
Every  interior  segment  of  G is  replaced  by  an  enter-leave 
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equivalent  elementary  segment  as  described  prior  to  Lemmas 
4.6,  4.7,  and  4.8.  The  following  Proposition  4.4  shows  that 
GA  is  Hamiltonian  if  and  only  if  G is. 

Proposition  4.4;  The  auxiliary  graph  GA  of  a rectilinear  (4, 
n)-grid  graph  G is  Hamiltonian  if  and  only  if  G is 
Hamiltonian. 

I will  now  present  a necessary  and  sufficient  condition 
for  the  auxiliary  graph  GA  to  be  Hamiltonian.  Let  GA  be  an 
auxiliary  graph  that  consists  of  elementary  end  segments  and 
interior  segments.  Recall  the  function  f that  assigns  a 
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Let  X 

be  an 

elementary  interior  segment 

of 

GA . 

The 

following  defines  a 6 x 6 path  enter-leave  incidence  (PELI^ 
matrix  A(X)  associated  with  the  segment  X.  The  (i,  j)  entry 
of  the  PELI  matrix  A is  1 if  there  exists  an  (i,  j) -bipath 
in  X,  otherwise  the  entry  is  zero.  Thus  A(i,  j)  = 1 if  and 
only  if  there  exists  a pair  of  paths  covering  the  vertices 
of  X,  and  for  the  entering  vertices  u and  v in  column  k,  and 
the  leaving  vertices  x and  y in  column  h,  ffc(u,  v)  = i and 
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fh(x'  Y)  = j*  Figure  4.15  shows  an  interior  segment  L with 
all  possible  (i,  j) -bipaths  and  the  corresponding  PELI 
matrix. 

Let  Y be  a left  elementary  end  segment  of  G* . The 
following  defines  a 6 x 6 path  leave  incidence  (PLI^  matrix 
A(Y)  associated  with  the  segment  Y.  The  j-th  column  of  the 


o o 


o o 


o o 


"0  0 0 0 0 0" 

0 0 0 0 0 0 

0 0 0 0 0 0 

A(L ) = 

0 0 0 1 0 1 

0 0 0 0 1 0 

0 0 0 1 0 1 

Figure  4.15.  An  example  of  an  interior  segment  L with 

possible  (i,  j) -bipaths  and  the  corresponding 
PELI  matrix 
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PLI  matrix  is  a vector  of  ones  if  and  only  if  there  exists  a 
j-path  in  Y.  In  other  words,  the  j-th  column  of  A(Y)  is  a 
vector  of  ones  if  and  only  if  there  exists  a single  path 
covering  the  vertices  of  Y,  and  for  the  leaving  vertices  u 
and  v in  column  k,  fk(u,  v)  = j.  Figure  4.16  shows  a left 
end  segment  P with  all  possible  j-bipaths  and  the 
corresponding  PLI  matrix. 

Let  Z be  a right  elementary  end  segment  of  GA . The 
following  defines  a 6 x 6 path  enter  incidence  fPEI^  matrix 
A(Z)  associated  with  the  segment  Z.  The  j-th  row  of  the  PEI 


A(P)  = 


"10  1 
10  1 
10  1 
10  1 
10  1 
10  1 


10  1" 
10  1 
10  1 
10  1 
10  1 
10  1 


Figure  4.16.  An  example  of  a left  end  segment  P with 

possible  j-paths  and  the  corresponding  PLI 
matrix 
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matrix  is  a vector  of  ones  if  and  only  if  there  exists  a j- 
path  in  Z.  In  other  words,  the  j-th  row  of  A(Z)  is  a vector 
of  ones  if  and  only  there  exists  a single  path  covering  the 
vertices  of  Z , and  for  the  entering  vertices  u and  v in 
column  h,  fft(u,  v)  = j.  Figure  4.17  shows  a right  end 
segment  L with  all  possible  (i,  j) -bipaths  and  the 
corresponding  PEI  matrix. 


O O 


o o 


'0  0 0 0 0 0' 

0 0 0 0 0 0 

0 0 0 0 0 0 

A(L)  = 

111111 
0 0 0 0 0 0 

_1  1 1 1 1 1 

Figure  4.17.  An  example  of  a right  end  segment  L with 
possible  j -paths  and  the  corresponding  PEI  matrix 


61 


A path  incidence  (PI)  matrix  is  a PELI,  PLI,  or  PEI 
matrix.  The  Pi  matrices  for  all  elementary  interior  and  end 
segments  are  given  in  Appendix  B.  Note  that  the  PI  matrices 
describe  the  enter-leave  and  enter/leave  patterns  of 
elementary  segments. 

The  path-compatibility  of  segments  of  G*  can  now  be 
easily  determined  by  simple  matrix  multiplication. 

Lemma — 4.9;  Two  adjacent  segments  of  the  auxiliary  graph 
are  path-compatible  if  and  only  if  the  product  of  the  two 
corresponding  PELI  matrices  is  nonzero. 

The  existence  of  a Hamiltonian  cycle  in  a rectilinear 
(4,  n)-grid  graph  can  now  be  easily  determined.  Let  G be 
separated  into  segments,  and  let  G*  be  the  auxiliary  graph 
associated  with  G.  Then  the  Hamiltonian  property  of  G is 
determined  by  the  following  Theorem  4.1. 

Theorem  4.1;  a rectilinear  (4,  n)-grid  graph  G is 
Hamiltonian  if  and  only  if  the  product  of  the  PI  matrices 
corresponding  to  the  segments  of  the  auxiliary  graph  G"  is 
nonzero. 

Proof;  Let  G be  a rectilinear  (4,  n)-grid  graph.  By  Lemma 
4.4,  graph  G is  Hamiltonian  if  and  only  if  the  auxiliary 
graph  associated  with  G^  is  Hamiltonian.  The  auxiliary  graph 
G"  is  Hamiltonian  if  and  only  if  it  consists  of  path- 
compatible  consecutive  segments.  By  Lemma  4.9,  consecutive 
segments  are  path-compatible  if  and  only  if  the  product  of 
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the  PI  matrices  is  positive.  Hence  G is  Hamiltonian  if  and 
only  if  the  product  of  the  PI  matrices  corresponding  to  the 
consecutive  segments  of  the  auxiliary  graph  G~  is  positive. 


Note  that  the  existence  of  a Hamiltonian  cycle  is 
determined  by  the  nonzero  product  of  the  PI  matrices.  The 
following  example  illustrates  Theorem  4.1.  Suppose  G 
contains  four  interior  segments  A,  G,  S,  and  H (as  labeled 
in  Figure  4.10),  in  that  order.  Since  the  product  of  the 
PELI  matrices  A G S H is  nonzero,  two  disjoint  paths 
spanning  segments  A,  G,  S,  and  H exist.  Figure  4.18 
illustrates  the  product  of  the  four  PELI  matrices  and 
corresponding  bipaths.  If  G contains  a left  end  segment  with 
a 1-path  and  a right  end  segment  with  a 4-path,  6-path,  or 
both,  then  the  product  of  the  PI  matrices  corresponding  to 
all  segments  of  G is  a matrix  of  ones,  and  G is  Hamiltonian. 
Note  that  the  product  of  the  PI  matrices  does  not  indicate 
how  to  construct  the  Hamiltonian  cycle. 

The  Hamiltonian  cycle  in  a rectilinear  (4,  n)-grid 
graph  can  be  constructed  from  the  j-paths  and  (i,  j) -bipaths 
through  the  adjacent  segments  as  follows.  Let  each  segment 
be  represented  by  a bipartite  graph  with  twelve  vertices 
arranged  in  two  columns  with  six  vertices.  Vertices  in  each 
column  are  labeled  from  1 to  6.  Vertex  i in  the  left  column 
is  adjacent  to  vertex  j in  the  right  column  if  and  only  if 
there  exists  a (i,  j) -bipath  in  the  corresponding  segment. 
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An  end  segment  is  also  represented  by  such  a bipartite 
graph.  In  this  case,  however,  vertex  j in  the  left  column  is 
adjacent  to  vertex  j in  the  right  column  if  and  only  if 
there  exists  a j -bipath  in  the  corresponding  segment. 
Construct  a composite  graph  that  consists  of  bipartite 
graphs  corresponding  to  the  segments  of  G and  including 


0 0 10  1' 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

Figure  4.18.  PELI  product  and  a corresponding  bipath 


A G S H = 


edges  (i,  i)  between  every  pair  of  bipartite  graphs 
corresponding  to  the  edges  that  join  the  pairs  of  paths  in 
the  two  segments.  Finally,  add  a source  and  a sink  node  each 
with  six  edges  to  "adjacent"  nodes.  Then  G is  Hamiltonian  if 
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and  only  if  there  exists  a path  between  the  source  and  the 
sink.  The  example  in  Figure  4.19  illustrates  this  procedure 
for  G with  segments  A,  G,  S,  and  H (A  and  H are  end  segments 
in  this  example). 

A simple  labeling  type  algorithm  to  find  a Hamiltonian 
cycle  in  a rectilinear  (4,  n)-grid  graph  is  shown  in  Table 


A G s H 


Figure  4.19.  Auxiliary  graph  indicating  a Hamiltonian  cycle 

4.1  below,  in  this  algorithm,  s and  t are  the  source  and 
sink,  respectively,  of  the  composite  graph  described  above. 
The  variable  Stage  refers  to  the  individual  bipartite 
graphs . 

After  performing  the  algorithm,  there  exists  a path 
from  s to  t corresponding  to  a Hamiltonian  cycle  in  G if  and 
only  if  l(t)  = TRUE.  If  l(t)  = TRUE,  pointers  to  the 
predecessors  can  be  used  to  trace  the  path  from  t to  s . 
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Almost  all  the  work  in  the  algorithm  is  in  the 
procedure  "Label."  In  the  worst  case  it  will  check  six 
possible  predecessor  nodes,  and  it  is  called  at  most  12  n 
times.  Thus  the  effort  to  use  it  is  at  most  72  n or  0(n). 
Checking  for  "Stop"  at  each  stage  is  0(1),  so  checking  for 
Stop  at  all  stages  is  0(n).  Thus,  the  order  of  the 


Table  4.1  Labeling  algorithm  to  find  a Hamiltonian 
cycle  in  a rectilinear  (4,  n)-grid  graph 


PROCEDURE  Label ( j ) ; 

IF  Node  j has  some  Predecessor  i and  l(i)  = TRUE 
THEN  l(j)  :=  TRUE;  {label  j> 

P(j)  :=  i {predecessor  of  j} 

ELSE  l(j)  :=  FALSE; 

{main} 

Stop  :=  FALSE; 

P(s)  :=  0; 

1 ( s ) :=  TRUE;  l(t)  :=  FALSE; 

Stage  :=  1; 

REPEAT 

FOR  every  Node  j in  first  half  Stage  DO 
Label ( j ) ; 

FOR  every  Node  k in  second  half  of  Stage  DO 
Label ( k) ; 

IF 

FOR  every  Node  k in  second  half  of  Stage 
1 ( k ) = FALSE 
THEN  Stop  :=  TRUE; 

Stage  :=  Stage  + 1; 

UNTIL  (Stop  OR  (Stage  = n + 1)); 

IF  NOT  Stop 
THEN  Label ( t ) . 


algorithm  is  proportional  to  the  number  of  stages  (segments 
in  G).  I remark,  that  the  results  for  the  algorithm  can  be 
proven  easily  by  induction  on  the  Stage  index.  What  would 
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have  to  be  proven  is  that  labels  are  assigned  correctly  at 
each  stage,  and  that  the  conclusion  is  correct. 


CHAPTER  5 

THE  PLACEMENT  SEQUENCE  FOR  NONHAMILTONIAN  TCHEBYSHEV 

(4,  N) -GRID  GRAPHS 


In  this  chapter  I consider  the  placement  sequence 
problem  for  a single  part  type  for  which  the  corresponding 
occupancy  graph  is  not  Hamiltonian.  Recall,  the  vertices  of 
an  occupancy  graph  represent  cells  of  a partitioned  PC  board 
containing  locations  where  a part  of  the  given  type  is  to  be 
placed.  The  occupancy  graph  itself  is  a subgraph  of  a 
rectangular  grid  graph  which  represents  the  partitioned  PC 
board.  Edges  in  the  rectangular  grid  graph  indicate  possible 
movements  of  the  Panasonic  Mkl  table  that  can  be  executed  in 
a time  period  less  than  the  time  window  w.  Any  placement 
sequence  that  places  all  parts  in  one  cell  of  the  PC  board 
before  moving  to  another  cell  and  returns  to  the  initial 
location  corresponds  to  a closed  walk  in  the  rectangular 
grid  graph. 

In  the  proposed  model,  the  absence  of  a Hamiltonian 
cycle  in  the  occupancy  graph  means  that  any  placement 
sequence  includes  a move  over  an  unoccupied  cell  or  over  a 
cell  previously  visited  and,  thus,  induces  a delay  time. 
Naturally,  I wish  to  find  a placement  sequence  with  a 
minimum  total  delay  time.  This  is  equivalent  to  solving  a 
traveling  salesman  problem  (TSP)  on  the  occupancy  graph.  The 
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occupancy  graph  discussed  in  the  preceding  sections  was 
defined  with  diagonal  edges  either  included  or  not  included. 
The  inclusion  of  diagonal  edges  corresponds  to  the 
Tchebyshev  metric,  while  the  absence  of  diagonal  edges 
corresponds  to  the  rectilinear  metric.  Since  the  two  metrics 
are  equivalent  under  a 45  degree  rotation  (and  scaling),  the 
following  discussion  is  valid  for  both  metrics.  For  ease  of 
exposition  I will  discuss  the  rectilinear  case,  and, 
whenever  necessary,  provide  results  for  both  metrics. 

Based  on  the  discussion  of  literature  in  Chapter  2 and 
given  the  special  character  of  the  (4,  n)-  grid  graph,  I 
propose  to  develop  the  following  generic  heuristic: 

Step  1:  Construct  a "hull"  of  the  occupancy  graph. 

Step  2:  The  vertices  on  the  "hull"  boundary  form  the 
initial  tour.  Insert  remaining  vertices. 

I will  base  the  specific  details  of  this  generic  heuristic 
on  a combination  of  the  convex  hull  and  1/2  band  insertion 
heuristic  modified  for  the  rectilinear  metric.  In  Section 
5.1  I will  discuss  several  "hulls"  including  the  'efficient 
set',  that  appear  in  the  literature.  In  Section  5.2  I will 
show  properties  of  the  'efficient  set'  that  are  analogous  to 
those  of  the  convex  hull  of  a set  of  points  in  the  Euclidean 
plane.  In  Section  5.3,  I will  summarize  the  efficient  set 
heuristic . 
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5.1  "Hulls"  in  the  Plane  with  Rectilinear  Metric 

In  this  section,  I will  describe  three  rectilinear 
"hull"  concepts  that  appear  in  the  literature  dealing  with 
the  rectilinear  metric.  I will  choose  one  of  these  to 
replace  the  convex  hull  in  the  convex  hull  heuristic  for  the 
traveling  salesman  problem. 

Figure  5.1  below  illustrates  three  different  "hull" 
concepts.  Figure  5.1a)  illustrates  the  hull  for  a set  S 
defined  by  Allison  and  Noga  [1984].  In  this  figure,  the 
hull  is  an  orthogonal  line  segment  seguence.  It  is  shown 
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Figure  5.1.  Rectilinear  "hulls" 

a)  hull  (rectilinear  segment) 

b)  rectangular  hull  (four  points  shown) 

c)  efficient  set 


twice  to  indicate  the  'return'  part  of  the  hull.  Consider 
the  set  S in  the  plane  and  let  ymin,  xmax,  ymax,  and  xmin  be 
the  points  in  S with  the  minimum  y-coordinate,  maximum  x- 
coordinate,  maximum  y-coordinate  and  minimum  x-coordinate. 
The  authors  define  a hull  to  be  "the  union  of  the  four 
paths 
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Pi  = (ymin, 


• • • r 


xmax) 


P2  = (xmax 


• • • » 


ymax) 


P3  = (ymax,  . ..,  xmin) 
P4  = (xmin,  ymin) 


• • • r 


where  each  P-^,  i = 1,  4,  is  an  orthogonal  line  segment 

sequence  that  is  maximal  with  respect  to  the  number  of 
points  it  may  contain  from  S."  The  Li  hull  is  analogous  to 
the  boundary  of  a convex  hull  in  the  sense  that  it  consists 
of  paths  with  a property  that  as  the  path  is  traversed  all 
points  of  S lie  on  the  path  or  to  one  side  of  it.  In  fact, 
each  point  of  S on  the  boundary  of  the  convex  hull  of  S will 
also  be  on  the  boundary  of  the  Li  hull.  Also,  the  Li  hull 
can  be  drawn  so  that  it  contains  the  minimum  amount  of  area, 
and  the  two  points  of  S which  are  farthest  apart  are 
contained  in  the  Li  hull.  On  the  other  hand,  the  authors 
point  out  that  the  order  in  which  points  appear  on  the  Li 
hull  is  not  necessarily  unique.  For  example,  the  hull  in 
Figure  5.1a)  can  be  described  as  A - B - C - D as  well  as  A 
- C - D - B.  Another  source  of  ambiguity  is  the  use  of  paths 
in  the  definition  of  the  L^  hull.  A path  is  described  as  a 
staircase-like  sequence  of  connected  orthogonal  line 
segments.  Since  such  a path  between  two  points  is  not 
unique,  as  the  authors  note,  the  Iq  hull  of  a set  of  two 
points  is  not  unique. 

Figure  5.1b)  shows  a rectangular  hull  defined  in  Love, 
Morris,  and  Wesolowsky  [1988].  Let  S be  a set  of  points  and 
J = {1,  ...,  n}  be  the  index  set  for  the  set  S.  The 
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rectangular  hull  is  defined  as  the  set  of  points  x = (xj, 
x2 ) that  satisfy  the  following  four  conditions 
simultaneously: 

X1  5 ajl  and  X2  a:  aj2  for  some  j e J, 

X1  ^ ajl  and  X2  ^ aj2  for  some  j E J, 

xi  ^ aji  and  X2  ^ aj2  for  some  j E J, 

X1  & ajl  and  x2  ^ aj2  for  some  j E J, . 

Geometrically,  the  rectangular  hull  is  a set  of  points  that 
cannot  be  reached  by  translation  of  any  quadrant  of  the 
coordinate  system  without  including  one  of  the  points  in  S 
as  an  interior  point  of  the  quadrant. 

Finally,  the  efficient  set  is  illustrated  in  Figure 
5.1c)  and  is  defined  in  Chalmet,  Francis,  and  Kolen  [1981]as 
follows.  For  the  set  S of  points  in  the  plane,  let  D(X)  be 
the  vector  of  rectilinear  distances  between  a point  X in  the 
plane  and  points  in  S.  A point  X dominates  point  Y if  each 
entry  in  D(X)  is  no  more  than  corresponding  entry  in  D(Y), 
and  at  least  one  entry  in  D(X)  is  strictly  less  than  the 
corresponding  entry  in  D(Y).  Those  points  that  are  not 
dominated  are  called  efficient  points.  The  efficient  set  is 
the  set  of  all  efficient  points. 

The  efficient  set  is  defined  for  the  rectilinear 
metric.  Furthermore,  there  exist  efficient  algorithms  that 
find  an  efficient  set  for  a given  set  of  points.  The  "arrow" 
algorithm  is  presented  in  Chalmet  et  al.,  who  also  present  a 
lowest  order  "coloring"  algorithm.  The  "coloring"  algorithm 
identifies  the  boundary  of  the  efficient  set.  In  the  case  of 
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the  occupancy  graph  it  identifies  the  leftmost  and  rightmost 
vertices  in  each  row. 

It  is  easy  to  see  that  an  efficient  set  can  be  defined 
analogously  for  the  Tchebyshev  metric.  Furthermore,  if  the 
plane  containing  the  points  is  rotated  45  degrees  and  scaled 
down  by  a factor  of  the  square  root  of  two,  the  efficient 
set  found  by  one  of  the  known  algorithms  under  the  inverse 
transformation  will  be  the  efficient  set  for  the  given  set 
of  points  with  the  Tchebyshev  metric.  A more  rigorous 
exposition  follows  and  is  due  to  Francis  [1993]. 

Let  di  (X,  Y)  and  doo(X,  Y)  denote  the  rectilinear  and 
Tchebyshev  distance  between  two  points  X and  Y, 
respectively,  and  let 


F:  E2  -►  E2:  (a,  b)  -*•  (l/2)(a  - b,  a + b) 

be  a transformation  of  a Euclidean  plane  which  rotates  every 
point  by  forty-five  degrees  followed  by  scaling  of  the 
coordinates  by  a factor  of  1/V2.  The  function  F is 
invertible  and  provides  a distance-preserving  transformation 
between  the  rectilinear  and  Tchebyshev  metrics.  Thus, 

di(X,  Y)  = doo(F(X),  F( Y)  ) (1) 

doo(U,  V)  = d1(F~1(U),  F-1  ( V)  ) (2) 

Let  P^,  i = 1,  ...,  n,  be  a set  of  points  (existing 
locations).  Let  Dr(X)  = (d^(X,  P^))  denote  a vector  of 
rectilinear  distances  of  a point  X from  the  points  p^. 
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Similarly,  let  Dt(X)  = (doo(F(X),  F(P-jJ)  denote  a vector  of 
Tchebyshev  distances  of  these  points  under  the 
transformation  F.  A point  X is  r-ef f icient  if  every  point  Y 
with  Dr(Y)  s Dr(X)  satisfies  Dr(Y)  = Dr(X),  and  a point  U is 
t-ef  f icient  if  every  point  V with  D-j-(V)  £ D-(-(U)  satisfies 
Dt(V)  = Dt(U).  Now,  let  Sr  and  St  be  the  set  of  all  r- 
eff icient  and  t-ef f icient  points,  respectively.  The  proof  of 
the  following  Proposition  5.1  is  straightforward  and  the 
interested  reader  can  find  it  in  Appendix  C. 

Proposition  5.1;  A point  U is  t-efficient  if  and  only  if  the 
point  F“1(U)  is  r-efficient. 

Proposition  5.1  implies  that  the  transformation  F 
preserves  efficiency  from  the  rectilinear  metric  to  the 
Tchebyshev  metric.  In  other  words,  every  point  in  F(Sr)  = 

{U:  there  exists  an  X in  Sr  such  that  F(X)  = U}  is  t- 
eff icient.  The  next  result  shows  the  equivalence  of  the 
efficient  sets  Sr  and  St  under  the  transformation  F. 

Proposition  5.2;  The  sets  Sr  and  St  are  equivalent  under  the 
transformation  F,  i.  e.  F(Sr)  = St. 

Proof : Let  U be  in  F(Sr).  Then  there  exists  a point  X in  Sr 
such  that  F ( X ) = U.  Hence,  X = F-1(U),  so  by  Proposition 
5.1,  it  follows  that  U is  in  S^.  Conversely,  If  U is  in  St, 
by  Proposition  5.1,  F-1(U)  is  in  Sr.  With  X = F~1(U)  in  Sr, 
and  F(X)  = u,  it  follows  that  U is  in  F(Sr).  ■ 


With  Proposition  5.2  in  mind,  I will  develop  the 
remaining  theory  in  this  chapter  using  the  rectilinear 
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metric.  The  rectilinear  metric  allows  easier  exposition,  and 
the  developed  results  will  be  true  for  both  the  rectilinear 
and  Tchebyshev  metrics.  It  follows  now  directly  that  given  a 
set  of  points  Q^,  i=  1,  ...,  n,  I can  find  the  set  of  t- 
efficient  points  using  the  following  procedure: 

1)  compute  the  locations  of  points  P-^  = F-l(Qi), 

2)  find  the  set  Sr  of  r-efficient  points  (using  a known 
algorithm  for  the  points  Pi), 

3)  take  St  = F(Sr) . 

The  efficient  set  has  some  of  the  same  geometric 
properties  as  the  boundary  of  the  convex  hull.  In 
particular,  an  optimal  tour  need  not  leave  a set  closely 
related  to  the  efficient  set  determined  by  the  occupancy 
graph.  Also  there  exists  an  optimal  rectilinear /Tchebyshev 
TSP  tour  in  which  the  relative  order  of  points  on  the 
boundary  of  the  efficient  set  is  preserved.  I will  establish 
these  properties  in  the  next  section. 

5.2  Properties  of  the  Efficient  Set 

The  efficient  set  is  a reasonable  choice  for  the 
determination  of  the  initial  sequence  for  the  TSP  on  an 
occupancy  graph  G for  at  least  two  reasons.  First,  an 
optimal  tour  need  not  leave  a region  that  is  closely  related 
to  the  efficient  set.  Second,  I will  show  that  a key 
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property  holds  for  the  rectilinear  TSP.  This  property  states 
that  there  exists  an  optimal  rectilinear/Tchebyshev  TSP  tour 
in  which  the  relative  order  of  points  on  the  boundary  of  the 
efficient  set  is  preserved. 

The  Euclidean  equivalent  of  these  properties  is  well 
established.  The  optimal  Euclidean  TSP  tour  is  contained  in 
the  convex  hull  formed  by  the  set  of  points.  The  second 
property  is  presented  in  Bellmore  and  Nemhauser  [1968]  as  a 
generalization  of  a result  due  to  Barachet  [1957]. 

To  show  analogous  geometric  properties  of  the  efficient 
set,  I describe  the  difference  between  the  geometric 
representation  of  Euclidean  and  rectilinear  metrics.  Let  dp, 
1 s p s oo,  be  the  lp-distance  metric;  then  L(U,  V)  = { X : 
dp(U,  X)  + dp(X,  V)  = dp(U,  V)}  is  the  geometric 
representation  of  a metric.  If  p = 2,  then  dp  is  the 
Euclidean  distance  metric  and  L(U,  V)  is  the  straight  line 
segment  between  the  points  U and  V.  If  p = 1,  then  dp  is  the 
rectilinear  distance  metric  and  L(U,  V)  is  the  rectangle 
with  vertical  and  horizontal  sides  whose  diagonal  is  the 
straight  line  segment  joining  the  points  U and  V.  In  the 
remainder  of  this  chapter  p = 1 and  L(U,  V)  denotes  the 
above  described  rectangle. 

I can  now  describe  a rectilinear  segment.  A rectilinear 
segment  between  two  points  in  the  plane,  U and  V,  is  a 
collection  of  sequences  of  horizontal  and  vertical  line 
segments  joining  them  so  that  the  total  length  of  the 
sequence  is  d^(U,  V).  The  quantity  d^(U,  V)  is  the  length  of 
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a rectilinear  segment  joining  the  points  U and  V. 

Rectilinear  paths  and  cycles  can  now  be  represented  as 
sequences  of  rectilinear  segments.  Clearly,  L(U,  V) 
represents  the  collection  of  paths  of  the  shortest 
rectilinear  total  length  between  two  points  U and  V,  and 
thus  is  analogous  to  the  straight  line  segment  representing 
the  Euclidean  distance.  Note  that  the  triangle  inequality 
holds  for  the  rectilinear  metric,  in  particular,  di(U,  V)  s 
di(U,  Z)  + d^Z,  V).  Furthermore,  equality  holds  if  and  only 
if  Z lies  within  the  rectangle  L(U,  V),  i.e.  d^U,  V)  = 
df(U,  Z)  + di(Z,  V)  if  and  only  if  Z is  in  L(X,  Y). 

In  Proposition  5.3  below,  I will  show  the  first 
property:  an  optimal  tour  need  not  leave  a region  that  is 
closely  related  to  the  efficient  set.  This  property  holds 
for  both  rectilinear  and  Tchebyshev  versions  of  the 
occupancy  graph.  Let  E denote  the  efficient  set  determined 
by  the  vertices  of  G,  and  let  E'  be  the  union  of  the 
efficient  set  E and  all  rectilinear  segments  L(i,  j)  between 
vertices  i and  j that  are  pairwise  adjacent  on  the  boundary 
of  E. 

In  Figure  5.  2 below,  the  efficient  set  E determined  by 
the  set  of  points  is  shown  shaded  and  includes  the  solid 
line  segments.  The  set  E'  consists  of  the  efficient  set  E 
and  the  regions  bounded  by  the  dashed  line  segments 

Let  S be  a set  of  points  with  an  efficient  set  E,  and 
let  A and  B be  two  points  in  S.  Lemmas  5.1  and  5.2  are  used 
to  prove  Proposition  5.3. 
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Figure  5.2.  The  efficient  set  E and  set  E' 

Lemma  5.1;  If  L(A,  B)  does  not  contain  any  points  of  S other 
than  A and  B,  then  L ( A,  B)  is  contained  in  E'. 

Proof : If  A and  B are  two  consecutive  points  on  the  boundary 
of  E , then  L ( A,  B)  is  in  E ' by  definition.  Otherwise,  L(A, 

B)  is  a either  a 'null  box'  or  a '2 -box'  contained  in  E.  ■ 

Lemma  5.2;  Let  A and  B be  two  points  in  S,  and  let  L(A,  B) 
be  the  rectilinear  segment  joining  points  A and  B of  length 
d^(A,  B).  There  exists  a rectilinear  path  joining  A and  B 
contained  in  E ' of  length  d^(A,  B). 

Proof ; If  L ( A,  B)  is  contained  in  E ' , there  is  nothing  to 
prove.  If  L ( A,  B)  is  not  contained  in  E'  then,  by  Lemma  5.1, 
there  exists  at  least  one  point  C in  S D L(A,  B)  other  than 
A and  B.  Choose  the  point  C that  is  rectilinearly  closest  to 
A.  Then  L ( A,  C)  contains  no  points  from  the  set  S fl  L ( A, 

B),  and  thus  L ( A,  C)  is  contained  in  E ' by  Lemma  5.1.  Since 
C is  in  L ( A,  B),  d;[  ( A,  B)  = d!(A,  C)  + di(C,  B)  so  the 
length  of  the  rectilinear  path  consisting  of  L ( A,  C)  and 
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L(C,  B)  is  di(A,  B).  if  L(C,  B)  is  contained  in  E ' , then  the 
proof  is  complete.  Otherwise,  repeat  the  same  argument  for 
the  rectilinear  segment  L(C,  B).  Since  S is  finite,  there 
exists  a rectilinear  path  joining  A and  B contained  in  E ' of 
length  exactly  d^A,  B).  ■ 

The  Proposition  5.3  now  follows.  Recall  that  a graph 
consists  of  a set  of  points  called  vertices  and  a set  of 
segments  joining  the  vertices  called  edges.  Thus,  I use  the 
terms  "vertex"  and  "point"  interchangeably;  the  same  is  true 
for  the  terms  "edge"  and  "rectilinear  segment". 

Proposition  5.3;  Let  G be  an  occupancy  graph.  There  exists 
an  optimal  traveling  salesman  tour  that  is  contained  in  E ' . 
Proof:  Let  G be  an  occupancy  graph,  and  let  E be  the 
efficient  set  determined  by  the  vertices  in  G.  Let  s be  any 
tour  of  G that  is  not  contained  in  E ' . If  an  edge  L(A,  B)  in 
the  tour  s is  not  contained  in  E ' , by  Lemma  5.2,  it  can  be 
replaced  by  a rectilinear  path  joining  A and  B contained  in 
E'  of  length  d^(A,  B).  if  all  segments  not  contained  in  E' 
are  replaced  in  this  fashion,  the  length  of  the  resulting 
sequence  s'  is  equal  to  the  length  of  s.  If  s is  optimal, 
then  s’  is  also  optimal,  and  it  is  contained  in  E'.  ■ 

The  boundary  of  an  efficient  set  has  a second  even  more 
significant  characteristic  that  makes  it  a good  choice  for 
the  initial  rectilinear  TSP  sequence:  there  exists  an 
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optimal  rectilinear  TSP  tour  in  which  the  relative  order  of 
points  on  the  boundary  of  the  efficient  set  is  preserved. 
The  Euclidean  version  of  this  property  is  given  below. 

Theorem:  Let  H be  a convex  hull  of  a set  of  points  in  2- 
dimensional  Euclidean  space.  There  exists  an  optimal 
(Euclidean)  tour  in  which  the  relative  order  of  points  on 
the  boundary  of  H is  preserved. 

The  proof  of  Barachet ' s Theorem  is  quite  easy  and 
depends  on  another  trivial  result:  "there  exists  an  optimal 
TSP  tour  that  does  not  cross  itself  when  the  distances 
satisfy  the  Euclidean  measure."  (This  result  does  not 
consider  collinearity  of  points,  and  fails  for  a tour  of 
four  collinear  points.)  Consider  a tour  in  Euclidean  space 
that  "crosses"  itself.  Let  AX  and  BY  be  two  intersecting 
straight  line  segments  in  the  tour.  The  length  of  the  tour 
does  not  increase  if  AX  and  BY  are  replaced  by  segments  AY 
and  BX  that  restore  the  tour.  Furthermore,  the  length  of  the 
tour  strictly  decreases  if  the  points  A,  B,  X,  and  Y are  not 
collinear.  Barachet ' s theorem  follows  trivially.  Both  of 
these  results  depend  on  the  distance  metric  satisfying  the 
triangle  inequality  and  serve  to  exclude  from  consideration 
those  sequences  that  do  not  satisfy  these  properties. 

In  the  following  I illustrate  the  concept  of 
intersecting  rectilinear  segments,  and  define  'weakly'  and 
'strongly'  intersecting  rectilinear  tours.  Note  that  a tour 
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is  a cycle.  Thus,  a tour  that  "crosses"  itself  is  a tour 
with  two  or  more  connecting  edges  (or  segments)  that 
intersect.  Consider  the  rectilinear  tours  in  Figures  5.3  and 
5.4  below.  In  each  of  these  figures,  as  well  as  in  the 
subsequent  ones,  a rectilinear  segment  joining  points  A and 
B is  shown  as  a staircase-like  sequence  of  connected 
orthogonal  line  segments  rather  than  as  the  rectilinear 
segment  L ( A,  B). 


A 


<> 

— • 

• B 

Y 


a)  b) 

Figure  5.3. a rectilinear  tour  that  "crosses"  itself  and  can 
be  made  shorter  by  "uncrossing." 
a) tour  A - b - x - Y 
bjtour  A - B - Y - x 


The  rectilinear  tour  in  Figure  5.3a)  has  intersecting 
segments  L(A,  Y)  and  L(B,  X).  The  tour  in  Figure  5.3b)  is 
obtained  by  "uncrossing"  the  intersecting  segments  and  has  a 
shorter  length,  in  contrast,  the  rectilinear  tour  in  Figure 
5.4a)  also  has  intersecting  segments  L ( A,  Y)  and  L(B,  X), 
and  the  tour  in  Figure  5.4b)  is  obtained  by  replacing  the 
segments  L ( A,  Y)  and  L(B,  X)  with  segments  L(A,  X)  and 
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a)  b) 

Figure  5.4.  A rectilinear  tour  that  "crosses"  itself  and 
cannot  be  made  shorter  by  "uncrossing." 

a)  tour  A - B - X - Y 

b)  tour  A - B - Y - X 


L ( B , Y).  Since  these  segments  also  intersect  the  new  tour 
has  the  same  length. 

With  reference  to  Figure  5.3a),  consider  the 
intersecting  segments  joining  vertices  A,  Y and  B,  X, 
respectively.  Note  that  the  rectangles  L(A,  Y)  and  L(B,  X) 
have  a nonempty  intersection.  The  nonempty  intersection  of 
L ( A,  Y)  and  L(B,  X)  is  analogous  to  the  single  point 
intersection  of  two  (straight)  line  segments  AY  and  BX 
representing  the  Euclidean  distances  between  four  (two 
pairs)  of  noncollinear  points.  As  illustrated  in  the  Figure 
5.3b),  another  pair  of  rectilinear  segments  joining  the 
pairs  A,  X and  B,  Y,  respectively,  exists.  In  this  case,  the 
intersection  of  the  corresponding  rectangles  L ( A,  X)  and 
L ( B , Y)  is  empty.  Again,  this  is  analogous  to  the  Euclidean 
case  where  the  line  segments  AX  and  BY  do  not  intersect. 
Note,  that  for  this  arrangement  of  the  points  A,  B,  X,  and 
Y,  the  following  relationship  holds  for  both  the  Euclidean 
and  rectilinear  metrics: 
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dp(A,  X)  + dp(B,  Y)  < dp(A,  Y)  + dp(B,  X),  p = 1 or  2. 
However,  if  the  four  points,  say  A,  B,  X,  and  Y are  located 
as  illustrated  in  the  Figure  5.4a)  and  b),  both  L ( A,  Y)  fl 
L ( B , X)  and  L ( A,  X)  fl  L(B,  Y)  are  nonempty  (and  equal).  In 
this  case, 

dX(A,  X)  + d!(B,  Y)  = dT  ( A,  Y)  + d^B,  X). 

If  I use  the  Euclidean  metric,  d2 , to  measure  the  segment 
lengths  in  this  case,  I have 

d2 (A,  X)  + d2(B,  Y)  < d2 (A,  Y)  + d2 (B,  X). 

This  example  illustrates  a difference  between  the 
rectilinear  and  the  Euclidean  metrics. 

Consider  the  tours  as  shown  in  Figures  5.3  and  5.4.  I 
will  say  that  two  rectilinear  segments  L(A,  Y)  and  L(B,  X) 
weakly  intersect  if 

L ( A,  Y)  fl  L(B,  X)  = L ( A,  X)  fl  L(B,  Y)  * 0. 

Figure  5.5  below  illustrates  a tour  with  two  weakly 
intersecting  rectilinear  segments  in  solid  lines  on  a 


Figure  5.5. A tour  with  two  weakly  intersecting  rectilinear 
segments  L (A,  X)  and  L ( B , Y) 
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dotted'  rectilinear  grid.  (Rectilinear  segments  joining 
points  are  shown  as  a staircase-like  sequence  of  connected 
orthogonal  line  segments  rather  than  as  a single  rectilinear 
segment).  The  rectangle  L(A,  Y)  n L ( B , X)  = L ( A,  X)  n L ( B , 

Y)  is  shaded. 


Figure  5.6. A tour  with  two  strongly  intersecting 

rectilinear  segments  L ( A,  X)  and  L (B,  Y) 


Two  rectilinear  segments  L ( A,  X)  and  L ( B,  Y)  strongly 
intersects  if  only  one  of  the  intersections  L(A,  X)  D L(B, 
Y)  or  L ( A,  Y)  n L ( B , X)  is  nonempty.  Figure  5.6  above 
illustrates  an  example  of  a tour  with  two  strongly 
intersecting  rectilinear  segments  in  solid  lines  on  a 
dotted'  rectilinear  grid.  Rectangle  L ( A,  X)  D L(B,  Y)  is 
shaded;  rectangle  L ( A,  Y)  D L(B,  X)  is  empty. 

A tour  or  a path  may  contain  several  pairs  of 
rectilinear  segments  that  intersect  either  weakly  or 
strongly.  I will  say  that  a tour  (or  a path)  s-intersects  if 
it  has  at  least  one  pair  of  segments  that  strongly 
intersects.  A tour  (or  a path)  w-intersects  if  every  pair  of 
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intersecting  segments  weakly  intersects.  A tour  (or  a path) 
sw- inter sects  if  it  either  s-intersects  or  w-intersects . 

Let  S and  S ' denote  the  intersections  of  rectilinear 
segments  between  points  A,  B,  X,  and  Y;  in  particular,  let 
s = L ( A,  Y)  n L(B,  X)  and  S'  = L ( A,  X)  D L (B,  Y).  The 
following  set  of  properties  illustrate  a simple  relationship 
between  the  rectangles  S and  S ' for  any  arrangement  of  four 
points  in  the  plane.  The  two  rectangles  are  either  equal  or 
at  least  one  is  empty.  The  properties  shown  below  were 
developed  jointly  with  Professor  Francis.  The  properties  are 
summarized  in  Table  5.1.  Recall  that 

S = L(A,  Y)  D L(B,  X)  and  S'  = L ( A,  X)  D L(B,  Y) 
and  let 

rl  = dl(A,  Y)  + dx(X,  B)  and  r2  = dx(A,  X)  + dx(Y,  B). 
The  following  provides  proofs  of  the  Properties  5.1  through 
5.6  given  in  the  Table  5.1  below. 

Property  5.1: 

i)  If  S is  nonempty  and  not  contained  in  S',  then  ri  > r2 

ii)  If  S'  is  nonempty  and  not  contained  in  S,  then  r^  < r2 . 
Proof : i)  if  s is  not  contained  in  S’,  there  exists  a point 
P in  S such  that  P is  not  in  S ' . Since  P is  in  S, 

rl  = di(A,  Y)  + dx (X,  B) 

= di (A,  P)  + dx(P,  Y)  + dx(X,  P)  + dx(P,  B). 

Since  P is  not  in  S ' , it  follows  from  the  triangle 
inequality. 
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Table  5.1.  Summary  of  the  relationship  between  S and  S' 
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r2  = dX(A,  X)  + di(Y,  B) 

< d1(A,  P)  + dx (P,  X)  + d!(Y,  P)  + dx(P,  B). 

Hence  > r2 . 

ii)  The  proof  is  analogous  to  part  i).  ■ 

Property  5.2: 

i)  If  S is  nonempty  and  S is  contained  in  S'  then  r^  = r2- 

ii) If  S'  is  nonempty  and  S'  is  contained  in  S then  r^  = r2 . 
Proof ; i)  Since  S is  not  empty  and  S is  contained  in  S', 
there  exists  a point  P in  the  intersection  of  S and  S ' . 
Since  P is  in  S, 

rj.  = di(A,  Y)  + d!(X,  B) 

= di(A,  P)  + dx(P,  Y)  + d!(X,  P)  + dx(P,  B). 

Since  P is  also  in  S', 

r2  = <h(A,  X)  + di(Y,  B) 

= dx  ( A,  P)  + d]_(P,  X)  + dx  ( Y,  P)  + d1(P/  B). 

Thus  ri  = r2* 

ii)  The  Proof  is  analogous.  ■ 

Property  5.3; 

i)  If  S is  nonempty  and  not  contained  in  S',  then  S'  is 
empty . 

ii) If  S'  is  nonempty  and  not  contained  in  S,  then  S is 
empty . 

Proof ; i)  Suppose  S is  not  contained  in  S ' . By  Property  5.1, 
rl  > r2 • The  set  s'  is  contained  in  S,  else  Property  5.1 


implies  r2  > . It  follows  that  S ' is  empty  else  Property 

5.2  implies  r^  = r2 . 

ii)  Proof  is  analogous.  ■ 
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Property  5.4; 

i)  If  S is  nonempty  and  S is  contained  in  S',  then  S = S'. 

ii) If  S’  is  nonempty  and  S'  is  contained  in  S,  then  S = S'. 
Proof : i)  Since  S is  nonempty  and  S is  contained  in  S', 
Property  5.2  implies  r^  = r2 . If  S is  a proper  subset  of  S', 
then  S'  is  not  contained  in  S,  and  Property  5.1  implies  r^  < 
r2,  which  is  impossible.  Thus  S = S'. 

ii)  Proof  is  analogous.  ■ 

Property  5.5:  S = S’  or  S n S ' is  empty. 

Proof : If  S is  nonempty  and  not  contained  in  S’,  Property 
5.3  implies  S'  is  empty.  If  S'  is  not  contained  in  S,  then, 
by  the  same  Property,  S is  empty.  It  follows  that  SOS'  is 
empty . 

If  S is  nonempty  and  contained  in  S ' , Property  5.4 
implies  S = S'.  If  S'  is  nonempty  and  S'  is  contained  in  S, 
then  also  S = S ' . 

If  either  S or  S ' is  empty,  then  clearly  SDS'  is 
empty . ■ 

Property  5.6; 

i)  If  S is  nonempty,  then  either  S = S'  or  S'  is  empty. 

ii)  If  S'  is  nonempty,  then  either  S = S'  or  S is  empty. 
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Proof ; i)  if  s is  nonempty,  then  we  cannot  have  S'  not 
contained  in  S,  for  then  S is  empty  by  Property  5.3.  Thus  S' 
is  contained  in  S.  If  S'  is  nonempty,  then  S = S'  by 
Property  5.4,  otherwise  S'  is  empty, 
ii)  Proof  is  analogous.  ■ 

Given  a rectilinear  tour  with  two  segments  L(A,  X)  and 
L ( B , Y)  that  sw-intersect , it  follows  from  the  definition 
and  Property  5.6  that  the  tour  either  w- intersects  or  s- 
intersects.  Property  5.1  further  implies  that  if  a tour  s- 
intersects  then  there  exists  a shorter  tour  that  can  be 
obtained  by  "uncrossing"  the  segments  that  strongly 
intersect.  Let  l(s)  denote  the  length  of  a tour  s.  The  Lemma 
5.  3 below  follows  directly  from  the  triangle  inequality. 

Lemma  5.3;  If  a rectilinear  tour  s s-intersects  then  there 
exists  a tour  s'  such  that  l(s)  > l(s'). 

Proof ; Let  s be  a tour  that  s-intersects.  Let  L(A,  Y)  and 
L(B,  X)  be  two  rectilinear  segments  of  s such  that  S is  not 
contained  in  S ' . By  Property  5.1, 

di(A,  Y)  + d1(B,  X)  < dx(A,  X)  + di(B,  Y). 

Thus  the  total  length  of  the  segments  L ( A,  X)  and  L(B,  Y)  is 
less  than  the  total  length  of  the  segments  L(A,  Y)  and  L(B, 
X).  Let  s ' be  a tour  obtained  from  s by  interchanging  the 
two  pairs  of  segments.  Then  l(s)  > l(s').  ■ 
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The  next  Corollary  is  analogous  to  Barachet ' s result 
that  an  "optimal  tour  does  not  intersect"  and  follows 
immediately  from  Lemma  5.3. 


Corollary — 5_,  1 ; An  optimal  tour  does  not  s-intersect. 

Proof:  Suppose  s is  an  optimal  tour  that  s-intersects . By 
Lemma  5.3,  there  exists  a tour  s'  such  that  l(s)  > l(s'). 
This  contradicts  the  assumption  that  s is  optimal.  ■ 

Theorem  5.1  below  establishes  that  there  exists  an 
optimal  tour  that  preserves  the  order  in  which  vertices 
appear  on  the  boundary  of  the  efficient  set.  Let  G be  a 
rectilinear  grid  graph  with  vertex  set  V,  and  let  E~  be  the 
boundary  of  the  efficient  set  of  the  vertex  set  of  G.  The 
next  lemma  and  corollary  are  used  in  the  proof  of  the 
theorem. 

Recall  G is  a rectilinear  occupancy  graph  with  a vertex 
set  V.  Note  that  V is  a finite  set.  Let  E be  the  efficient 
set  of  the  vertex  set  V,  and  let  E~  be  the  subset  of  points 
m V lying  on  the  boundary  of  E.  Assume  the  vertices  Xi  in 
E~  are  labeled  in  the  order  in  which  they  appear  in  E~. 
Recall  for  any  two  points  X and  Y in  the  plane,  L(X,  Y)  is 
the  rectilinear  line  segment  that  can  be  geometrically 
represented  as  a the  set  of  all  points  inside  the  rectangle 
with  opposite  vertices  X and  Y and  with  vertical  and 
horizontal  sides.  Let  E'  = E U {L(X±,  Xi+1):  x±  E E-,  i = 

1,  — , k}  U {L(Xi,  Xk):  X]_,  Xk  E E~} . The  set  E'  is  the 
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union  of  the  efficient  set  E and  all  rectilinear  segments 
L(xi'  Xj)  that  are  pairwise  adjacent  in  E~.  Note  that  E~  is 
a finite  set  of  points  in  the  plane  with  E~  C v,  while  E 
and  E ' are  closed  regions  in  the  plane  with  E C E ' . In 
Figure  5.7  below,  the  efficient  set  E is  shown  shaded;  E'  is 
the  region  inside  the  dashed  polygon,  and  E-  is  the  set  of 
points  shown. 


Figure  5.7.  Sets  E,  E ' , and  E~ 

Let  A,  B,  c,  and  D be  four  points  in  E~  named  in  the 
order  in  which  they  appear  in  E~.  The  next  lemma  shows  that 
two  rectilinear  paths  from  A to  C and  from  B to  D that  are 
contained  in  E'  must  sw-intersect . 

Lemma  5.4:  Let  A,  B,  C,  and  D be  four  points  in  E-  C v 
named  in  the  order  in  which  they  appear  in  E~.  Any  two 
rectilinear  paths  PAC  and  PBD  from  A to  C and  from  B to  D, 
respectively,  that  are  contained  in  E'  sw-intersect. 
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Proof:  Let  PAC  and  PBD  be  two  rectilinear  paths  in  E ' . 
Replace  each  rectilinear  segment  L(Xi/  X j ) between  two 
adjacent  points  in  E~  and  each  rectilinear  segment  in  pAq 

and  PBD  by  a (Euclidean)  straight  line  segment  as  shown  in 
Figure  5.8. 


Figure  5.8.  Paths  PAC,  pbd,  nAC  and  nBD 


Let  r be  the  set  of  straight  line  segments  obtained  by 
the  replacement  of  rectilinear  segments  L(Xi,  Xj)  between 
two  adjacent  points  in  E~,  and  let  nAC  and  nBD  be  the  set 
of  straight  line  segments  obtained  by  the  replacement  from 
PAC  and  PBDf  respectively.  Note  that  r is  a closed  polygonal 
curve,  and  both  nAC  and  nBD  are  polygonal  arcs.  Let  int(T) 
be  the  interior  region  of  T.  it  is  clear  from  the  definition 
of  E'  that  E C r U int(T) . 

Since  A,  B,  C,  and  D are  ordered  in  E-  and,  thus,  in  r, 
points  A and  C decompose  r into  two  polygonal  arcs  and 
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a2  so  that  B is  in  and  D is  in  a 2 (or  vice  versa).  Note 
that  r — ai  U a 2 • Let  and  T2  be  two  closed  polygonal 
curves  such  that  = ct^  U nAc  and  T2  = a2  U nAc.  Consider 
the  straight  line  segments  in  nBD  starting  at  the  vertex  B. 
Since  V is  a finite  set,  there  exists  a straight  line 
segment  in  nBD  that  joins  a pair  of  vertices,  say  and 
z2>  such  that  Z]_  e Tj  U intfTi)  and  z2  e ext(T]_).  The 
notation  ext(r^)  denotes  the  region  exterior  to  the  curve 
rl-  The  Jordan  Curve  Theorem  (Dieudonne  1969)  implies  that 
the  straight  line  segment  from  Zj  to  Z2,  [Z]_,  Z2],  must 
intersect  Tj . Since  [Zj,  Z2]  is  contained  in  F U int(T),  it 
must  intersect  nAc . Let  [Y]_,  Y2  ] be  the  straight  line 
segment  in  nAC  that  intersects  [zlr  Z2].  Recall,  these  line 
segments  are  the  diagonals  of  rectilinear  segments  L(Y]_,  Y2 ) 
and  L(Z]_,  Z2).  Clearly,  if  diagonals  of  two  rectangles 
intersect,  the  rectangles  intersect.  Thus,  L(Y]_,  Y2)  fl 
L(zl/  z2)  * afid  PAq  and  PBD  sw-intersect . ■ 

Corollary  5.2;  Let  A,  B,  C,  and  D be  the  order  in  which  the 
four  points  appear  in  E~.  If  the  points  A and  C and  points  B 
and  D are  adjacent  in  an  optimal  tour  s contained  in  E ' , 
then  the  tour  s sw-intersects . 

Proof:  Since  L ( A,  C)  and  L(B,  D)  sw-intersect  by  Lemma  5.4, 
the  result  is  immediate.  ■ 

To  prove  Theorem  5.1  I show  that  for  any  optimal 
sequence  s that  does  not  preserve  the  relative  order  of 
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points  in  which  they  appear  on  the  boundary  of  the  efficient 
set  there  exists  another  optimal  sequence  that  does.  I prove 
that  s sw-intersects , and  since  an  optimal  sequence  does  not 
s_-'-ntersect • s w-intersects . Thus,  it  can  be  rearranged  into 
a sequence  s ' of  equal  length  with  vertices  in  the  same 

order  in  which  they  appear  on  the  boundary  of  the  efficient 
set. 

Theorem  5.1;  There  exists  an  optimal  rectilinear  TSP  tour  in 
which  the  relative  order  of  vertices  on  the  boundary  of  the 
efficient  set  is  preserved. 

Proof:  Let  s be  an  optimal  tour  in  which  the  relative  order 
of  points  on  the  boundary  of  an  efficient  set  is  not 
preserved.  By  Proposition  5.3,  assume  s is  contained  in  E ' . 

Since  s does  not  preserve  the  order  of  vertices  on  the 
boundary  of  the  efficient  set,  there  exist  at  least  three 
vertices  B,  C,  and  D,  that  appear  on  the  boundary  of  the 
efficient  set  E and  in  s in  different  order.  Assume  the 
order  of  the  vertices  in  s is  B - D - C,  while  the  order  in 
which  they  appear  on  the  boundary  of  E is  B - C - D.  This  is 
illustrated  in  Figure  5.9.  in  the  figure,  the  boundary  of 
the  efficient  set  E is  the  rectangle  shown  and  may  also 
contain  other  vertices  of  G other  than  pictured.  Tour  s is 
partially  shown  passing  through  vertices,  in  order,  B,  i,  j, 
D,  C,  k,  1,  a,  B. 

The  vertex  A is  the  next  boundary  vertex  after  C is 
visited  by  s.  Let  PAC  be  the  rectilinear  path  induced  by  s 
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Figure  5.9. A rectilinear  tour  that  does  not  preserve  the 
order  of  vertices  in  E. 

between  vertices  A and  C,  and  let  PBD  be  the  rectilinear 
path  induced  by  s between  vertices  B and  D.  (Note  that  PAC 
and  PBD  are  shown  as  a staircase-like  sequence  of  connected 
orthogonal  line  segments  rather  than  as  the  sequence  of 
rectilinear  segments.)  By  Lemma  5.4,  the  paths  PAC  and  PBD 
sw- intersect . 

Let  i and  j be  two  vertices  adjacent  in  PAC,  and  let  k 
and  1 be  two  vertices  adjacent  in  PBD  such  that  the  segments 
L(i,  j)  and  L(k,  1)  contain  the  sw-intersection  of  PAc  and 
PBD  as  illustrated  in  Figure  5.11.  Consider  the  tour  s' 
partially  described  with  vertices  in  order 

a ...  B ...  l,  k,  ...  C ......  D ...  j , 1 ...  A ... 

obtained  from  uncrossing  L(i,  j)  and  L(k,  1).  Corollary  5.1 
implies  s weakly  intersects,  hence  the  length  of  the  tour  is 
preserved. 

Thus,  a new  sequence  of  the  same  length  can  be 
constructed  from  the  current  sequence  for  any  triple  of 
vertices  that  appear  in  different  order  on  the  boundary  than 
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in  the  current  sequence.  Since  s is  optimal,  the  final 
sequence  is  optimal  and  preserves  the  relative  order  of 
vertices  in  which  they  appear  on  the  boundary  of  the 
efficient  set.  ■ 

The  following  result  for  Tchebyshev  tours  follows 
easily. 

Corollary  5.3:  There  exists  an  optimal  Tchebyshev  TSP  tour 
in  which  the  relative  order  of  vertices  on  the  boundary  of 
the  efficient  set  is  preserved. 

Proof : Since  the  rectilinear  and  Tchebyshev  efficient  sets 
are  equivalent  by  Proposition  5.2,  the  proof  follows 
directly  from  Theorem  5.1.  ■ 

The  results  in  this  section  show  that  the  efficient  set 
has  some  of  the  same  properties  as  the  convex  hull,  and, 
thus,  the  vertices  on  the  boundary  of  the  efficient  set 
serve  as  a reasonable  initial  TSP  tour. 

5.3  The  Efficient  Set  Heuristic 

The  properties  of  the  efficient  set  shown  in  the 
preceding  section  establish  that  the  boundary  of  an 
efficient  set  is  a reasonable  initial  sequence  for  a 
traveling  salesman  tour  in  a Tchebyshev  grid  graph.  The 
efficient  set  of  an  occupancy  graph  consists  of  certain 
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boxes  and  connecting  edges  formed  by  vertical  and  horizontal 
line  segments  passing  through  the  vertices  of  the  graph 
(Chalmet  et  al.  [1981]).  in  the  case  of  an  occupancy  graph 
that  is  a Tchebyshev  (4,  n)-grid  graph,  it  is  easy  to  see 
that  the  vertices  in  the  highest  and  the  lowest  rows  are  on 
the  boundary  of  the  efficient  set.  Thus  selecting  the 
boundary  of  the  efficient  set  followed  with  the  widest 
insertion  heuristic  is  similar  in  nature  to  the  1/2  band 
insertion  heuristic  studied  by  Bozer,  Schorn,  and  Sharp 
[1990].  Based  on  the  reported  performance  of  1/2  band 

Table  5.2.  Efficient  set  heuristic 

Step  1:  Construct  the  efficient  set  E of  the  occupancy 

graph  and  let  the  boundary  of  E be  the  initial 
seguence  s. 

Step  2:  Insert  the  remaining  vertices  into  s using  the 

widest  angle  insertion  heuristic. 

Step  3:  (Optional)  Use  the  2-opt  heuristic  to  improve  s. 


insertion  heuristic,  I propose  the  heuristic  outlined  in 
Table  5.3  for  the  placement  problem  for  a single  part  type 
if  the  associated  occupancy  graph  is  not  Hamiltonian. 

Based  on  the  theory  discussed  in  Chapters  4 and  5,  I can 
now  determine  the  placement  seguence  for  a single  part  type. 
Given  the  occupancy  graph  associated  with  the  part  type,  I 
can  check  whether  the  occupancy  graph  is  2-connected,  and 
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hence  Hamiltonian.  If  it  is  Hamiltonian,  I use  the 
Hamiltonian  cycle  as  the  (optimal)  placement  sequence. 
Otherwise,  I use  the  Efficient  Set  Heuristic  to  find  a 
(suboptimal)  placement  sequence.  By  doing  so,  I hope  to 
minimize  the  total  time  required  to  place  all  the  parts  of 
single  type  on  the  PC  board. 


CHAPTER  6 

SEQUENCE  MENDING  PROBLEM  AND  PC  BOARD  ASSEMBLY  HEURISTIC 

In  this  chapter  I will  describe  a greedy  heuristic  to 
solve  the  sequence  mending  subproblem  of  the  PC  board 
assembly  problem.  A solution  to  this  subproblem  determines 
the  allocation  of  reels  containing  the  parts  to  the  feeder 
mechanism.  Each  reel  in  the  carrier  of  the  Panasonic  Mkl 
machine  contains  a tape  with  parts  of  a single  type.  The 
greedy  heuristic  I propose  in  this  chapter  places  all  parts 
of  one  type  on  the  PC  board  before  it  places  a part  of 
another  type.  The  heuristic  described  combines  the  placement 
sequences  found  by  the  heuristic  described  in  the  previous 
chapter  into  a grand  assembly  sequence.  The  sequence  of  part 
types  in  the  grand  assembly  sequence  determines  the  order  in 
which  reels  with  parts  are  to  be  placed  in  the  carrier 
slots.  Thus  if  the  mending  causes  part  of  type  A,  B,  C,  and 
D to  be  placed  sequentially,  then  the  reels  with  these  parts 
will  be  placed  in  the  carrier  in  this  order  also.  The  grand 
assembly  sequence  itself  determines  the  order  in  which  parts 
are  to  be  placed  on  the  PC  board.  In  the  remainder  of  this 
chapter  I present  an  overview  of  my  proposed  method  to  solve 
the  PC  board  problem  and  describe  its  computational 
performance . 
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6 . 1 Sequence  Mending  Problem 

In  this  section  I will  describe  the  sequence  mending 
problem.  Recall  that  a placement  sequence  for  a given  part 
type  is  a cycle  of  vertices  in  the  occupancy  graph 
representing  part  locations  on  the  PC  board.  The  length  of  a 
sequence  is  given  by  the  total  distance  between  consecutive 
elements  of  the  sequence  in  the  occupancy  graph.  The 
sequence  mending  problem  can  be  described  as  follows.  Given 
a set  of  cyclical  sequences  of  part  locations  on  the  PC 
board,  I wish  to  order  the  sequences  in  such  a way  - ideally 
- that  the  total  length  of  the  grand  assembly  sequence  is 
minimized.  For  illustration,  consider  the  set  of  cycles  C^, 
i = 1 , . . . , 4 , representing  the  placement  sequences . The 
cycles  can  be  ordered  into  a grand  assembly  sequence  by 
removing  an  arbitrary  edge  from  each  cycle  and  connecting 
(mending^  them  into  a single  cycle  in  some  order  as 
illustrated  in  Figure  6.1. 


Figure  6.1.  Mending  of  four  cycles 
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An  adaptation  of  the  nearest  neighbor  heuristic  for 
solving  the  sequence  mending  subproblem  proceeds  as  follows. 


Let  oi,  i = 1,...,  p,  be  the  set  of  placement  sequences  that 
need  to  be  mended.  Since  is  a permutation  of  a finite  set 
of  indices  representing  the  locations  of  parts  of  type  i, 
oi-1(k)  and  o-jjk)  denote  the  elements  in  the  sequence 
immediately  preceding  and  following  the  element  k in  a^. 
Select  a 'seed'  sequence,  say  o\,  and  determine  the  element 
p in  it  with  the  shortest  (Tchebyshev)  distance  to  an 
element  q in  any  other  sequence,  say  02 . The  elements  p and 
q are  referred  to  as  an  exit  element  of  and  an  entry 
element  of  02  for  this  mending  operation,  respectively.  The 
two  sequences  03  and  02  are  now  considered  as  mended  and 
become  the  new  seed  sequence.  To  determine  the  exit  element 
of  the  new  seed  sequence,  consider  the  elements  02 (q)  or 
°2-1(c3)*  Choosing  one  of  02  (q)  or  02-1(q)  as  an  exit  element 
allows  the  sequence  02  to  remain  intact  in  the  grand 
assembly  sequence.  Select  that  element,  say  a2(q),  which  is 
closest  to  a sequence  that  has  not  yet  been  mended  to  become 
the  new  exit  element.  The  procedure  is  repeated  until  all 
sequences  are  mended.  The  final  step  is  to  determine  the 
entry  element  of  aj_.  For  this  purpose,  03  (p)  or  o^“l(p)  is 
chosen,  whichever  is  closer  to  the  exit  element  of  the  last 
mended  sequence. 

Consider  the  following  example.  Let  o\,  02,  and  03 
denote  the  placement  sequences  for  three  part  types . 
Individual  part  locations  for  each  sequence  are  shown  in 
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Figure  6.2.  The  placement  sequences  are  shown  as  'solid  - 
dotted'  cycles.  Select  sequence  03  as  the  'seed'  sequence. 
The  element  is  closest  to  a sequence  not  yet  mended,  03 . 
Element  p^  is  the  exit  element  of  03,  and  element  q^  is  the 
entry  element  in  02-  Next,  examine  the  elements  adjacent  to 
qi  in  «2,  q2  and  . The  element  q2  is  closest  to  a sequence 
not  yet  mended,  03.  The  element  q2  becomes  the  exit  element 
of  the  new  ' seed ' sequence  03  U 02 , and  the  element  r^ 
becomes  the  new  entry  element  in  03.  Next,  examine  elements 
adjacent  to  r^  in  o 3,  r2  and  . Since  all  sequences  are 
mended  now,  one  of  r2  and  will  be  mended  to  one  of  the 
elements  adjacent  to  P3  in  03,  P2  and  P4.  Since  r2  in  03  is 
closest  to  P2  in  03,  r2  is  the  new  exit  element  and  P2  is 
the  new  entry  element.  The  resulting  grand  assembly  sequence 


Figure  6.2.  An  example  of  grand  assembly  sequence 

(°l)  P2/  P3  • P4  / Plf  (°2 ) <3lr  33'  <32  / (°3)  rl>  r4  / r3/ 
r2  is  shown  in  Figure  6.2  as  the  'solid-dashed'  cycle.  The 
parts  should  be  placed  in  the  carrier  in  order  o^,  03,  and 
°3  • 


102 


In  the  next  section,  I incorporate  this  greedy  heuristic 
into  a PC  board  assembly  problem  heuristic  that  provides  the 
grand  assembly  sequence  as  well  as  indicates  how  the  reels 
with  parts  should  be  allocated  on  the  carrier. 

6.2  PC  Board  Assembly  Heuristic 

In  a class  of  manufacturing  problems  discussed  here,  N 
parts  are  to  be  placed  on  a rectangular  PC  board  during 
assembly  at  prespecified  locations  by  an  automatic  robot. 
During  manufacture  of  large  quantities  of  PC  boards  it  is 
desirable  to  minimize  the  total  assembly  time  per  PC  board. 
The  assembly  consists  of  a series  of  pick-and-place 
operations  during  which  a part  is  retrieved  from  the  supply 
and  another  part  is  placed  on  the  PC  board.  In  order  to 
reduce  the  complexity  of  this  problem  I formulated  the  PC 
board  assembly  problem  in  Chapter  2 as  a clustered  traveling 
salesman  problem.  The  'cluster'  constraint  forces  the 
placement  of  all  parts  of  one  type  before  parts  of  any  other 
type  are  placed.  This  eliminates  any  potential  delay  due  to 
a late  carrier.  Further  analysis  of  a class  of  automatic 
robots  leads  to  a breakup  of  the  clustered  traveling 
salesman  problem  into  a set  of  Hamiltonian  cycle  problems  in 
Tchebyshev  (4,  n)-grid  graphs.  In  the  context  of  the  PC 
board  assembly  problem  the  Tchebyshev  (4,  n)-grid  graphs  are 
also  referred  to  as  occupancy  graphs.  I proposed  the 
following  PC  board  assembly  heuristic  to  solve  this  problem. 
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For  each  part  type  determine  whether  the  associated 
occupancy  graph  is  2 -connected.  If  it  is,  then  it  is 
Hamiltonian.  In  this  case,  the  Hamiltonian  cycle 
corresponding  to  a placement  sequence  can  be  found  using  an 
adaptation  of  the  convex  hull  heuristic.  If  the  occupancy 
graph  is  not  2-connected,  then  a traveling  salesman  tour 
corresponding  to  a placement  sequence  can  be  determined  by 
another  adaptation  of  the  convex  hull  heuristic.  Both  of 
these  techniques  are  based  on  the  properties  of  the 
occupancy  graphs,  and  they  are  described  in  Chapters  4 and 
5.  After  the  placement  sequences  are  determined  for  all  part 
types,  the  grand  assembly  sequence  is  determined  by  sequence 
mending,  an  adaptation  of  the  nearest  neighbor  heuristic. 

The  sequence  mending  heuristic  is  described  in  the  previous 
section.  The  PC  board  assembly  heuristic  is  summarized  in 
Table  6.1  below. 

Table  6.1.  The  PC  board  assembly  heuristic 

Procedure  Biconnected; 

IF  2 -Connected THEN  Find  HamiltonianCycle 

ELSE  Find  Eff  Set  Tour 

Procedure  Mend. 


Next,  I discuss  the  implementation  of  the  PC  board 
assembly  heuristic.  The  Biconnected  procedure  determines 
which  of  the  occupancy  graphs  are  2-connected.  This  property 
is  determined  by  checking  for  each  vertex  in  the  occupancy 
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graph  whether  it  is  a cut  vertex  or  not.  It  is  well  known 
that  a graph  is  2 -connected  if  and  only  if  it  has  no  cut 
vertex.  There  are  other,  perhaps  more  elegant  algorithms, 
however,  they  are  more  complex  to  code. 

Once  an  occupancy  graph  for  a part  type  is  labeled  2- 
connected,  the  Hamiltonian_Cycle  procedure  computes  a 
placement  seguence.  This  procedure  follows  the  outline  given 
in  detail  in  Section  4.1.3.  If  the  occupancy  graph  is  not  2- 
connected  and,  therefore,  not  Hamiltonian,  then  the 
placement  sequence  is  computed  in  the  Eff_Set_Tour  procedure 
described  in  Chapter  5 . 

To  implement  this  procedure,  I have  adapted  an  existing 
code  by  Dale  Kostamo  for  the  "coloring"  algorithm  (Chalmet 
et  al.  [1981]).  The  coloring  algorithm  finds  the  efficient 
set  for  a given  set  of  points  with  rectilinear  distance 
measure.  Since  the  Tchebyshev  and  rectilinear  metric  are 
equivalent  under  rotation  and  scalar  multiplication,  the 
Eff_Set_Tour  procedure  transforms  the  Tchebyshev  occupancy 
graph  to  an  equivalent  grid  graph  with  the  rectilinear 
metric.  The  coloring  algorithm  then  determines  the  boundary 
of  the  efficient  set  by  identifying  the  leftmost  and 
rightmost  vertices  in  each  row  of  the  graph,  and  the 
solution  is  transformed  back  to  the  original  occupancy 
graph.  The  placement  sequence  is  completed  by  inserting  the 
remaining  vertices  from  the  second  and  third  row  into  the 
boundary  using  the  widest  angle  insertion  heuristic. 
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After  a placement  sequence  is  found  for  every  part  type, 
procedure  Mend  mends  the  sequences  into  a single  grand 
assembly  sequence.  This  procedure  is  an  implementation  of  a 
greedy  heuristic.  The  procedure  initializes  the  grand 
assembly  sequence  by  setting  it  equal  to  the  placement 
sequence  of  the  first  part  type.  The  placement  sequence 
closest  to  the  last  location  in  the  current  grand  assembly 
sequence  is  added  to  the  current  grand  assembly  sequence. 

The  order  (clockwise  or  counterclockwise)  in  which  the 
sequence  is  added  is  determined  by  the  closest  distance  to 
the  next  placement  sequence.  This  procedure  is  repeated 
until  all  placement  sequences  are  added  to  the  grand 
assembly  sequence.  After  the  grand  assembly  sequence  is 
known,  the  total  assembly  time  is  computed. 

In  the  next  Section,  I compare  the  total  assembly  time 
obtained  by  this  heuristic  to  the  same  measure  obtained  from 
a clustered  traveling  salesman  heuristic. 

6_^3 — Performance  of  the  Heuristic 

In  this  section  I present  some  computational  results. 
Tables  6.2  through  6.4  contain  a summary  of  the  solution 
times  obtained  by  the  PC  board  assembly  heuristic. 

The  PC  board  assembly  problem  solution  I propose  is  a 
heuristic  and  does  not  guarantee  an  optimal  solution. 

However  the  computational  results  show  it  gives  near-optimal 
solutions.  Furthermore,  the  heuristic 


approach  is  easily 
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implemented  in  a personal  computer  environment.  The  code  I 
developed,  written  in  Turbo  Pascal,  is  included  in  Appendix 
D.  It  was  implemented  on  386SX  (25MHz)  machine. 

For  the  evaluation  of  the  PC  board  assembly  heuristic,  I 
consider  PC  boards  with  varying  dimensions  and  number  of 
part  types.  For  each  chosen  PC  board  specifications,  I 
generate  random  part  locations  on  the  PC  board  and  collect 
data  on  the  total  length  of  the  grand  assembly  sequence 
provided  by  the  heuristic.  The  average  values  are  shown  in 
the  tables  below.  The  lower  bound  described  in  Chapter  3 is 
also  shown  in  each  table  as  well  as  the  relative  error. 

Each  table  also  contains  average  lengths  of  tours 
obtained  from  equivalent  clustered  traveling  salesman 
problems.  For  each  PC  board,  the  distances  between  parts  are 
computed  from  the  occupancy  graph  representing  the  PC  board. 
To  solve  the  clustered  TSP,  I add  to  the  distances  between 
locations  of  parts  of  different  types  a 'sufficiently  large' 
quantity.  This  allows  an  ordinary  TSP  algorithm  to  solve  a 
clustered  TSP  problem  (Chisman  [1975]).  To  solve  the 
clustered  TSP  test  problem  I used  a branch  and  bound  TSP 
algorithm  encoded  by  T.  Volgenant  and  W.B.  van  den  Hout 
[1990] . 

I now  discuss  the  performance  of  the  PC  board  assembly 
heuristic.  In  the  tables  below  , I report  the  following 
average  values:  the  length  of  the  grand  assembly  tour,  the 
lower  bound,  the  relative  error  and  the  length  of  the  best 
TSP  tour  obtained  from  the  branch  and  bound  algorithm. 
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In  the  first  set  of  results  I tested  the  heuristic  on  PC 
boards  of  varying  dimensions.  I chose  PC  boards  with  side 
ratios  4:4  (square),  4:6  and  4:8  (rectangles).  The  number  of 
part  types  and  the  fill  ratio  (average  number  of  parts  per 
PC  board  cell)  is  held  constant  as  much  as  possible  with 
three  parts  per  cell.  Locations  of  individual  parts  on  the 
PC  board  were  generated  randomly,  and  a sample  size  of  20  PC 
boards  of  each  type  was  generated.  Table  6.2  contains  PC 
board  specifications  chosen  in  each  case. 

Table  6.2  Experimental  PC  board  specifications 


PC  Board 
Dimensions 
Ratio 

Number 

of 

cells 

Number 
of  part 
types 

Number  of 
parts  per 

type 

Number  of 
parts  per 
cell 

Total  number 
of  parts 

4 : 4 

16 

3 

16 

3 

48 

4 : 6 

24 

3 

25 

3 

75 

4 : 8 

32 

3 

30 

3 

90 

Table  6.3  below  shows  the  summary  of  results.  The  first 
column  shows  the  PC  board  dimension  ratios.  The  heuristic 
performed  satisfactorily,  with  worst  average  error  only 
12.1  percent  for  the  rectangular  PC  board  with  side  ratio 
4:8.  In  comparison,  the  branch  and  bound  algorithm  reduced 
the  gap  between  a lower  bound  and  the  'best  sequence  found' 
in  less  than  300  iterations  and  gave  no  additional 
improvement  for  another  1000  iterations.  At  this  point 
(after  about  5 minutes  of  running  time),  I aborted  the 
program  since  the  relative  gap  was  small,  as  is  shown  in  the 
fourth  column  of  the  table.  This  result  illustrates  a 
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definite  advantage  of  the  PC  board  assembly  heuristic.  Its 
running  times  are  negligible  for  larger  problems  when 
compared  to  the  running  time  of  TSP  algorithms;  furthermore, 
the  small  relative  error  is  reasonable  even  for  practical 
applications . 


Table  6.3.  Performance  of  heuristic  solution  on  PC  board 
with  varying  dimensions 


PC  Board 
Dimensions 
Ratio 

Length  of  grand 
assembly  tour 

Lower 

Bound 

Average 

error 

(*) 

Worst 

error 

(*) 

Length  of 
best  TSP  tour 

4 : 4 

516 

480 

8.1 

12.5 

511 

4 : 6 

816 

750 

8.8 

12.0 

806 

4 : 8 

1009 

900 

12.1 

21.1 

984 

In  the  next  set  of  results  I tested  the  PC  board 
assembly  heuristic  on  rectangular  and  square  PC  boards  with 
increasing  number  of  part  types.  The  rectangular  PC  board 
dimensions  were  set  at  ratio  4:6.  The  number  of  part  types 
varies  form  one  to  eight.  The  total  number  of  parts  on  the 
PC  board  increases  with  the  number  of  part  types.  However,  I 
hold  the  number  of  parts  per  type  constant  at  25.  Hence  the 
number  of  parts  varies  with  the  number  of  part  types  from  75 
to  200.  The  rectangular  PC  board  has  24  cells,  and  the 
average  number  of  parts  per  cell  varies  from  3.1  to  8.3;  the 
results  are  summarized  in  Table  6.4.  The  sguare  PC  board  has 
16  cells,  and  the  average  number  of  parts  per  cell  varies 
from  4.7  to  12.5;  the  results  are  summarized  in  Table  6.5. 
Again,  locations  of  individual  parts  on  the  PC  board  were 
generated  randomly,  and  a sample  size  of  20  PC  boards  of 
each  type  was  generated. 
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Table  6.4.  Performance  of  heuristic  solution  on  rectangular 
PC  board  with  varying  number  of  parts 


Number  of 
part  types 

Length  of 
grand 

assembly  tour 

Lower 

Bound 

Average 

error 

(») 

Worst 

error 

(*) 

Length  of 
best  TSP  tour 

3 

796 

750 

6.1 

9.3 

802 

4 

1089 

1000 

8.9 

12.0 

1082 

5 

1358 

1250 

8.7 

14.4 

1356 

6 

1618 

1500 

7.9 

11.3 

1630 

7 

1930 

1750 

10.3 

12.6 

1950 

8 

2158 

2000 

7.9 

10.5 

2244 

Table  6.5  summarizes  similar  results  for  the  square  PC 
board.  These  results  suggest  that  the  PC  board  assembly 
heuristic  outperforms  the  branch  and  bound  algorithm  on  a 
rectangular  PC  board  with  75  or  more  parts.  In  most  cases, 
the  heuristic  provides  tours  of  shorter  length  than  the 
branch  and  bound  algorithm,  and  it  provides  them  in  only  a 
few  seconds  on  a 386SX  machine.  In  comparison,  the  branch 
and  bound  algorithm  often  terminated  for  one  of  the  reasons 
described  in  the  previous  case.  Note  that  the  average  worst 
possible  error  of  the  heuristic  solution  does  not  change 
significantly  with  the  increasing  number  of  parts.  The 
square  PC  board  is  more  densely  populated  than  the 
rectangular  board  in  the  previous  case.  The  smaller  error 
reflects  this  property. 

The  computational  results  presented  in  this  section 
illustrate  the  advantages  of  the  PC  board  assembly 
heuristic.  The  most  significant  advantage  is  its  running 
time.  For  larger  problems  with  number  of  parts  exceeding  one 
hundred,  it  provides  a good  suboptimal  solution  in  matter  of 


Table  6.5.  Performance  of  heuristic  solution  on  square 
PC  board  with  varying  number  of  part  types 
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Number  of 
part 
types 

Length  of 
grand  assembly 
tour 

Lower 

Bound 

Average 

error 

(%> 

Worst 

error 

(%> 

Length  of  best 
TSP  tour 

3 

763 

750 

1.7 

2.7 

776 

4 

1040 

1000 

4.0 

7.0 

1036 

5 

1290 

1250 

3.2 

4.8 

1304 

6 

1535 

1500 

2.3 

2.7 

1558 

7 

1798 

1750 

2.7 

3.4 

1824 

8 

2043 

2000 

2.1 

3.0 

2106 

seconds  on  a personal  computer.  By  comparison,  a branch  and 
bound  TSP  algorithm  runs  for  a long  time  and,  in  some  cases, 
even  without  reaching  a comparable  suboptimal  solution 
before  aborting.  The  worst-case  errors  found  for  the  PC 
board  assembly  heuristic  are  small,  not  exceeding  twelve 
percent . 


6.4  Conclusions  and  Future  Research 

The  manual  scheduling  of  the  operations  required  to 
assemble  a PC  board  is  very  time  consuming  for  the  PC  board 
industry.  Automated  PC  board  assembly  software  is  needed  in 
this  industry. 

In  this  work  I study  the  PC  board  assembly  problem  on  a 
class  of  automated  component  placement  machines.  I present  a 
heuristic  algorithm  which  provides  the  grand  assembly 
sequence  for  the  PC  board  assembly.  This  sequence  provides 
the  order  in  which  the  parts  should  be  placed  on  the  PC 
board  as  well  as  the  assignment  of  part  types  to  the  feeder 
mechanism  of  the  machine.  Computational  results  indicate 
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that  the  sequence  obtained  is  very  close  to  the  optimal 
solution. 

The  real  advantage  of  this  approach  to  automating  the  PC 
board  assembly  is  its  application  on  a personal  computer. 
While  TSP  algorithms  may  provide  optimal  solutions,  the  PC 
board  assembly  heuristic  presented  in  this  work  runs  much 
faster,  providing  adequate  suboptimal  solutions.  An  operator 
could  use  a personal  computer  on  the  assembly  floor,  and, 
given  specifications  of  a PC  board,  he  could  find  and 
implement  a good  assembly  sequence  quickly. 

The  future  work  based. on  my  findings  is  both  in  the 
practical  and  theoretical  area.  The  PC  board  assembly 
heuristic  should  be  further  tested  and  then  implemented  in 
real  PC  board  assembly  environment.  Its  impact  on 
productivity  can  then  be  evaluated. 

There  are  many  opportunities  for  future  investigation  of 
the  theoretical  aspects  of  this  work  as  well.  I have  studied 
interesting  combinatorial  problems  (Hamiltonian  cycles, 
ordinary  and  clustered  traveling  salesman  tours)  in  (4,  n)- 
grid  graphs.  It  would  be  interesting  to  investigate  these 
problems  in  general  (m,  n)-grid  graphs. 

Another  problem  that  can  be  further  investigated  is  the 
problem  of  mending  several  placement  sequences . I have 
presented  a greedy  heuristic  for  this  purpose.  However,  the 
special  structure  of  the  (4,  n)-grid  graphs  may  allow  for  a 
better  approach. 


APPENDIX  A 

CHAPTER  4 RESULTS  AND  PROOFS 


This  Appendix  contains  the  results  and  proofs  from 
Section  4.2 

Lemma  4.4:  Let  G(V,  E)  be  a Hamiltonian  rectilinear  (4,  n)- 
grid  graph.  A Hamiltonian  cycle  intersects  any  end  segment  S 
of  G in  a single  path  that  covers  all  vertices  of  S. 

Proof : Let  G(V,  E)  be  a Hamiltonian  rectilinear  (4,  n)-grid 
graph  with  Hamiltonian  cycle  C,  and  let  S be  an  end  segment 
in  G.  Clearly  C covers  S.  Suppose  that  C intersects  S in 
more  than  one  path.  A Hamiltonian  cycle  in  a rectilinear 
grid  graph  is  a Jordan  curve,  and,  as  a consequence  of  the 
Jordan  Closed  Curve  Theorem,  any  cut  from  "outside"  to 
"outside"  of  the  Hamiltonian  cycle  cuts  the  cycle  in  an  even 
number  of  places.  In  particular,  if  column  k is  the  column 
of  S adjacent  to  an  interior  segment  or  the  other  end 
segment  of  G,  then  a cut  through  column  k must  intersect  C 
in  an  even  number  of  vertices.  Thus  (since  C intersects  S in 
at  least  two  paths ) column  k of  S contains  at  least  four 
vertices  of  C that  are  not  adjacent  in  C else  C is  not  a 
Hamiltonian  cycle.  Hence  S is  a 4 -segment.  Furthermore, 
since  C is  a cycle,  column  k + 1 adjacent  to  S in  the  next 
segment  must  contain  four  vertices  of  G that  are  adjacent  in 
C to  the  vertices  in  column  k.  This  contradicts  the 
assumption  that  S is  a segment,  i.e.  the  largest  set  of 
adjacent  columns  with  vertices  of  G in  the  same  rows.  ■ 
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Lemma  4.5:  Let  G(V,  E)  be  a Hamiltonian  rectilinear  (4,  n)- 
grid  graph.  A Hamiltonian  cycle  intersects  any  interior 
segment  S of  G in  exactly  two  disjoint  paths  that  cover  all 
vertices  of  S. 

Proof : Let  G(V,  E)  be  a Hamiltonian  rectilinear  (4,  n)-grid 
graph  with  Hamiltonian  cycle  C,  and  let  S be  an  interior 
segment  in  G.  Clearly  C covers  S.  Suppose  that  C intersects 
S in  more  than  two  disjoint  paths.  A Hamiltonian  cycle  in  a 
rectilinear  grid  graph  is  a Jordan  curve,  and,  thus  as  a 
consequence  of  the  Jordan  Closed  Curve  Theorem,  any  cut  from 
"outside"  to  "outside"  of  the  Hamiltonian  cycle  cuts  the 
cycle  in  an  even  number  of  places.  In  particular,  if  columns 
k and  1 are  columns  of  S adjacent  to  other  segments  of  G, 
then  a cut  through  column  k (and  also  cut  through  column  1) 
must  intersect  C in  an  even  number  of  vertices.  Thus  (since 
C intersects  S in  at  least  three  disjoint  paths)  the  columns 
k and  1 contain  at  least  six  vertices  of  C that  are  not 
adjacent  in  C else  C is  not  a Hamiltonian  cycle.  Therefore, 
at  least  one  end  column,  say  the  right  end  column  k,  of  S 
contains  four  vertices  that  are  not  adjacent  in  C.  Hence  S 
is  a 4-segment.  Furthermore,  since  C is  a cycle,  column  k + 

1 adjacent  to  column  k in  the  next  segment  must  contain  four 
vertices  of  G that  are  adjacent  in  C to  the  vertices  in 
column  k.  This  contradicts  the  assumption  that  S is  a 
segment  i.e.  the  largest  set  of  adjacent  columns  with 
vertices  of  G in  the  same  rows . ■ 
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Lemma  4.6;  Any  interior  2-segment  in  G is  enter-leave 
equivalent  to  an  interior  elementary  2-segment. 

Proof ; Interior  2-segments  contain  only  one  pair  of  disjoint 
paths  covering  the  vertices  of  the  segment  regardless  of  the 
number  of  columns,  an  (i,  i) -bipath  for  some  i.  ■ 

Lemma  4.7;  Any  interior  3-segment  in  G is  enter-leave 
equivalent  to  an  interior  elementary  3-segment. 

Proof ; Let  S be  an  interior  3-segment.  Consider  the 
following  cases. 

Case  1:  If  the  vertices  of  S are  not  in  adjacent  rows, 
then  an  (i,  j) -bipath  consists  of  paths  P and  Q;  one  of 
these  paths,  say  P,  contains  vertices  in  the  isolated  row, 
and  the  path  Q can  only  be  a "square-wave"  path  that 


Figure  A.l.  An  (i,  j) -bipath  in  a 3-segment  with  vertices 
in  nonadjacent  rows 

contains  the  remaining  vertices.  Due  to  its  square-wave 
characteristic,  the  path  Q repeats  its  shape  every  two 
columns.  Clearly  then,  all  interior  3-segments  with  vertices 
in  nonadjacent  rows  and  an  even  number  of  columns  have  the 
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same  enter-leave  pattern,  and  the  same  is  true  about  all 
such  3-segments  with  an  odd  number  of  columns.  It  follows, 
that  any  interior  3-segment  with  vertices  in  nonadjacent 
rows  is  enter-leave  equivalent  to  an  interior  elementary  3- 
segment  with  two  columns  or  one  column. 

Case  2 : Let  S be  an  interior  3-segment  with  vertices  in 
adjacent  rows.  If  P and  Q are  two  disjoint  paths  covering 
the  vertices  of  S,  then  one  or  the  other  starts  from  the 
entering  vertices  as  a "square-wave"  path,  however,  in  any 
column  the  other  path  may  become  the  "square-wave"  path 
provided  S has  at  least  two  columns.  For  example,  consider  a 
segment  S with  (1,  2,  3) -pattern  columns  shown  in  Figure 
A. 2.  If  the  entering  vertices  of  P and  Q are  in  the  first 
and  third  row,  respectively,  P leaves  the  first  column 
through  the  second  row  vertex  and  Q leaves  the  first  column 
through  the  third  row  vertex.  Path  P is  the  "square-wave" 
path.  In  the  second  column,  path  P has  to  include  the  first 
row  vertex  and  leave  the  column  through  the  first  row.  Path 
Q continues  through  the  third  row  vertex.  In  the  third 
column  either  P or  Q can  include  the  second  row  vertex.  If  P 
does,  then  P continues  as  the  "square-wave"  path.  If  the 
path  Q includes  the  second  row  vertex,  then  Q becomes  the 
"square-wave"  path.  (In  Figure  A. 2,  this  occurs  in  the  fifth 
column.)  In  either  case,  P and  Q leave  the  third  column 
through  adjacent  vertices.  In  general,  if  the  entering 
vertices  of  P and  Q are  adjacent,  then  the  paths  P and  Q 
leave  an  odd  column  through  nonadjacent  vertices  and 
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Figure  A.  2.  An  (i,  j) -bipath  in  a 3-segment  with  vertices 
in  adjacent  rows 

leave  an  even  column  through  a pair  of  adjacent  vertices. 
Similarly,  if  the  entering  vertices  of  P and  Q are 
nonadjacent,  then  the  paths  P and  Q leave  an  odd  column 
through  a pair  of  adjacent  vertices  and  leave  an  even  column 
through  nonadjacent  vertices.  Clearly  then,  all  interior  3- 
segments  with  vertices  in  adjacent  rows  and  an  even  number 
of  columns  have  the  same  enter-leave  pattern,  and  the  same 
is  true  about  such  3-segments  with  an  odd  number  of  columns. 
It  follows,  that  any  interior  3-segment  with  vertices  in 
adjacent  rows  is  enter-leave  equivalent  to  an  interior 
elementary  3-segment  with  two  columns  or  one  column.  ■ 

Lemma  4.8;  Any  interior  4 -segment  in  G is  enter-leave 
equivalent  to  an  interior  elementary  4 -segment. 

Proof ; If  S is  an  interior  4-segment  with  less  than  six 
columns,  then  S is  an  elementary  4-segment. 

Let  S be  an  interior  4-segment  with  at  least  six 
columns.  It  is  easy  to  see  that  {(1,  1),  (1,  3),  (1,  4),  (1, 
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6), 

(3, 

1), 

(3, 

3), 

(3, 

4), 

(3, 

6), 

(4, 

1), 

(4, 

3), 

(4, 

4), 

(4, 

6), 

(6, 

1), 

(6, 

3), 

(6, 

4), 

(6, 

6), 

(2, 

2), 

(5, 

5)> 

is 

the  bipath  enter-leave  pattern  of  a 4 -segment  with  four 
columns,  and  {(1,  1),  (1,  3),  (1,  4),  (1,  6),  (3,  1),  (3, 

3),  (3,  4),  (3,  6),  (4,  1),  (4,  3),  (4,  4),  (4,  6),  (6,  1), 
(6,  3),  (6,  4),  (6,  6),  (2,  5),  (5,  2)}  is  the  enter-leave 
pattern  of  a 4-segment  with  five  columns.  Note,  that  the 
segment  enter-leave  patterns  do  not  include  bipath  enter- 
leave  configurations  with  i or  j equal  to  2 or  5 with  two 
exceptions  in  each  pattern.  I will  demonstrate  that  the 
enter-leave  pattern  of  the  segment  S contains  all  bipath 
enter-leave  configurations  listed  above  for  a 4-segment  with 
four  or  five  columns  depending  on  whether  S has  an  even  or 
odd  number  of  columns,  respectively.  Furthermore,  I will 
show  that  the  enter-leave  pattern  of  the  segment  S contains 
no  other  bipath  enter-leave  configurations  (i.e.  those  with 
i or  j equal  to  2 or  5 that  are  not  included  in  the  enter- 
leave-patterns  above).  Thus,  by  showing  that  S has  the  same 
segment  enter-leave  pattern  as  a 4 -segment  with  four  or  five 
columns,  I will  prove  that  six  or  more  column  4-segments 
with  an  even  or  odd  numbers  of  columns  are  enter-leave 
equivalent  to  an  elementary  4 -segment  with  four  or  five 
columns,  respectively. 

Since  an  elementary  4 -segment  with  2 columns  contains 
an  (i,  i) -bipath  for  all  i = 1,...,  6,  it  can  be  attached  to 
any  4-segment  to  obtain  a larger  4-segment  so  that  the  new 
segment  contains  an  (i,  j) -bipath  if  the  original  4-segment 
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does.  It  follows  that  the  enter-leave  pattern  of  a 4-segment 
S with  an  even  (odd)  number  of  columns  includes  all  bipath 
enter-leave  configurations  in  the  pattern  of  an  elementary 
4-segment  with  four  (five)  columns.  The  construction  of  a 
(1,  3) -bipath  in  a six  column  4-segment  based  on  a (1,  3)- 
bipath  in  an  four  column  elementary  4 -segment  is  illustrated 
in  Figure  A. 3. 


Figure  A. 3.  The  construction  of  a (1,  3) -bipath  in  a six 

column  4-segment  based  on  a (1,  3) -bipath  in  an 
four  column  elementary  4 -segment 

Let  u and  v be  the  entering  vertices  of  an  (i,  j)- 
bipath  in  S with  i = 2.  The  case  when  i = 5 is  symmetric. 
Since  the  entering  vertices  of  the  (i,  j) -bipath  are  in  rows 
1 and  3,  the  (i,  j) -bipath  can  only  be  two  "square-wave" 
paths  as  illustrated  in  the  Figure  A. 4.  It  follows  that  the 
only  possible  value  of  j is  2 or  5,  depending  on  whether  S 
has  an  even  or  odd  number  of  columns,  respectively. 
Similarly,  if  x and  y are  the  leaving  vertices  of  an  (i,  j)- 
bipath  in  S with  j = 2,  the  "square-wave"  pattern  is  forced 
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Figure  A. 4.  The  "box-wave"  (5,*) -bipath  in  an  interior  4- 
segment 

on  the  two  disjoint  paths  covering  the  vertices  of  S,  and 
the  only  possible  values  of  i are  2 and  5 depending  on 
whether  S has  an  even  or  odd  number  of  columns.  Again,  the 
case  when  j = 5 is  symmetric.  Hence  the  enter-leave  pattern 
of  S contains  no  bipath  enter-leave  configurations  with  i or 
j equal  to  2 or  5 except  those  included  in  the  enter-leave 
pattern  of  a corresponding  elementary  4 -segment  with  four  or 
five  columns. 

Since  S has  the  same  enter-leave  pattern  as  an  elementary  4- 
segment  with  four  or  five  columns,  S is  enter-leave 
equivalent  to  an  elementary  4-segment.  ■ 

Corollary  4.3;  Any  end  segment  with  a two  vertex  pattern  is 
enter /leave  equivalent  to  an  elementary  end  segment. 

Proof : The  proof  follows  directly  from  Lemma  4.6.  A j-path 
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can  be  easily  obtained  from  a (j,  j) -bipath  if  and  only  if  j 
represents  a pattern  with  vertices  in  adjacent  rows.  ■ 

Corollary  4.4;  Any  end  segment  with  a three  vertex  pattern 
is  enter/leave  equivalent  to  an  elementary  end  segment  with 
one  or  two  columns . 

Proof ; Let  S be  an  end  segment  with  a three  vertex  pattern. 
Clearly,  a j-path  can  be  obtained  from  an  (i,  j) -bipath  if 
and  only  if  the  vertices  of  S are  in  adjacent  rows.  The 
result  follows  from  Lemma  4.7.  ■ 

Corollary  4.5:  Any  end  segment  with  a four  vertex  pattern  is 
path  equivalent  to  an  end  segment  with  one  or  two  columns . 
Proof ; Let  S be  a left  end  segment  with  a four  vertex 
pattern.  The  proof  for  right  end  segments  is  symmetric.  It 
is  easy  to  see  that  an  end  segment  with  a four  vertex 
pattern  and  two  columns  contains  a j-path  with  leaving 
vertices  in  every  possible  combination  of  rows  except  in 
rows  1 and  3 , and  rows  2 and  4 . In  other  words , such  an  end 
segment  contains  no  2-path  or  5-path.  From  the  proof  of 
Lemma  4.8  it  follows  that  any  end  segment  S with  three  or 
more  columns  and  a four  vertex  pattern  contains  the  same  j- 
paths.  Any  such  segment  also  contains  no  2-path  or  5-path, 
since  the  leaving  vertices  force  a double  "square-wave" 
which  cannot  cover  all  vertices  in  the  segment  in  a single 
path.  Hence  all  end  segments  with  more  than  two  columns  and 
a four  vertex  pattern  have  the  same  enter/leave  pattern.  It 
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follows  that,  if  S is  an  end  segment  with  a four  vertex 
pattern,  then  S is  enter /leave  equivalent  to  an  end  segment 
with  one  or  two  columns . ■ 


Proposition  4.4;  The  auxiliary  graph  GA  of  a rectilinear  (4, 
n)-grid  graph  G is  Hamiltonian  if  and  only  if  G is 
Hamiltonian. 

Proof ; Let  G be  a Hamiltonian  rectilinear  (4,  n)-grid  graph. 
The  Hamiltonian  cycle  determines  the  enter-leave 
configuration  of  each  interior  segment  and  the  enter/leave 
configuration  of  the  end  segments.  Since  each  segment  of  GA 
is  either  enter-leave  or  enter/leave  equivalent  to  a 
corresponding  segment  of  G,  a Hamiltonian  cycle  in  GA  can  be 
constructed  using  (i,  j) -bipaths  and  j-paths  in  the  segments 
of  GA  with  the  same  configurations  as  in  the  corresponding 
segments  of  G. 

Conversely  and  similarly,  if  GA  is  Hamiltonian,  a 
Hamiltonian  cycle  in  G can  be  constructed  using  (i,  j)- 
bipaths  and  j-paths  in  the  segments  G that  have  the  same 
configurations  as  those  given  by  the  Hamiltonian  cycle  in 
the  segments  of  GA . The  enter-leave  configuration  preserving 
construction  for  extending  (i,  j) -bipaths  is  described  in 
the  proof  of  Lemma  4.8.  ■ 

Lemma  4.9;  Two  adjacent  segments  of  the  auxiliary  graph  GA 
are  path-compatible  if  and  only  if  the  product  of  the  two 
corresponding  PELI  matrices  is  nonzero. 
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Proof : The  proof  follows  from  the  definition  of  matrix 
multiplication.  Given  two  PI  matrices,  A and  B, 
corresponding  to  two  segments  SA  and  S3  their  product  C is 
defined  as  follows: 

C(i,  j)  = 2k=l,...,  n A(i,  k)  B(k,  j). 

Clearly,  C(i,  j)  is  positive  (at  least  one)  if  and  only  if 
both  A(i,  k)  and  B(k,  j)  are  positive  for  some  k. 

Thus,  if  C is  a nonzero  product  of  PI  matrices  A and  B, 
then  there  exist  integers  i and  j such  that  C(i,  j)  s 1,  and 
thus  A(i,  k)  and  B(k,  j)  are  both  positive  for  some  k.  This 
implies  that  the  segment  associated  with  the  matrix  A,  SA, 
contains  a k-path  or  a (k,  j) -bipath  P with  leaving  vertices 
u and  v in  column  h such  that  fft(u,  v)  = k.  Similarly,  there 
exists  a k-path  or  a (k,  j) -bipath  Q contained  in  the 
segment  corresponding  to  the  matrix  B,  S3,  with  the  entering 
vertices  x and  y in  column  h + 1 such  that  fh+i(x,  y)  = k. 
Hence  P and  Q are  compatible,  and  SA  and  S3  are  path- 
compatible  . 

Conversely  and  clearly,  if  P and  Q is  a combination  of 
two  paths  or  bipaths  in  the  path-compatible  segments  SA  and 
SB,  joined  through  the  leaving  vertices  u and  v in  column  h 
of  P and  entering  vertices  x and  y in  column  h = 1 of  Q, 
then  the  entries  A(i,  k)  and  B(k,  j)  are  positive  for  some 
k = fh(u,  v)  = fh+i(x,  y)  and  some  i and  j.  This  implies  the 
product  of  A and  B has  a positive  entry  in  cell  (i,  j)  , and 
is,  therefore,  nonzero.  ■ 
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Theorem  4.1:  A rectilinear  (4,  n)-grid  graph  G is 

Hamiltonian  if  and  only  if  the  product  of  the  PI  matrices 
corresponding  to  the  segments  of  the  auxiliary  graph  is 


nonzero. 

Proof ; Let  G be  rectilinear  (4,  n)-grid  graph.  By  Lemma 
4.4,  graph  G is  Hamiltonian  if  and  only  if  the  auxiliary 
graph  associated  with  G~  is  Hamiltonian.  The  auxiliary  graph 
G^  is  Hamiltonian  if  and  only  if  it  consists  of  path- 
compatible  consecutive  segments . By  Lemma  4.9,  consecutive 
segments  are  path-compatible  if  and  only  if  the  product  of 
their  PI  matrices  is  positive.  Hence  G is  Hamiltonian  if  and 
only  if  the  product  of  the  PI  matrices  corresponding  to  the 
segments  of  the  auxiliary  graph  G*  is  positive.  ■ 


APPENDIX  B 
PI  MATRICES 


This  Appendix  contains  the  PELI  matrices  associated  with 
interior  elementary  segments.  The  matrices  A,  B,  C,  etc. 
correspond  to  the  segments  illustrated  in  Figure  4.10. 
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The  PLI  matrices  of  those  elementary  segments  that 
contain  a j-path  for  some  j are  given  below.  The  PI  matrices 
are  transposed  PLI  matrices. 
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APPENDIX  C 

EFIICIENT  POINTS  IN  RECTILINEAR  AND  TCHEBYSHEV  METRICS 


This  Appendix  contains  Proposition  5.1  and  its  proof 
from  Section  5.1 

Proposition  5.1;  A point  U is  T-efficient  if  and  only  if  the 
point  F-1(U)  is  r-efficient. 

Proof:  Let  U be  T-efficient,  and  let  X = F“1(U).  To  show 
that  X is  r-efficient,  for  every  Y with  rY  ^ rx,  it  must  be 
shown  that  rY  = rx  Now, 
rY  s rx 

implies 

di(Y,  Pi)  s dx (X,  P±),  i = 1,  ...,  n. 

It  follows  by  (2) 

doo(F(Y),  F(P±))  <;  doo(F(X),  F(P±)),  i = 1,  ...,  n 

or,  equivalently, 

d«(F(Y),  F(Pi))  <;  doo(U,  F(Pi)),  i = 1,  ...,  n. 
Since  U is  T-efficient, 

doo(F(Y),  F(P±))  = doo(U,  F(P±)),  i = 1,  ...,  n, 

and,  it  follows 

di(Y,  Pi)  = di(F-l(U),  P±),  i = 1,  ...,  n 
or 

di(Y,  P±)  = di(X,  Pi),  i = 1,  ...,  n. 

Thus  X is  r-efficient. 

Conversely,  let  X = F-1(U)  be  r-efficient.  For  every  V 
with  tv  ss  ty,  it  must  be  shown  that  tv  = ty.  Now, 

doo(V,  F(Pi))  s doo(U,  F(P±)),  i = 1,  ...,  n 

implies 

di(F-!(V),  P±)  ^ d1(F-1(U),  P±),  i = 1,  ...,  n. 


130 


However,  since  F~1(U)  is  r-efficient. 


di(F"l(V),  P±)  = d!(F-l(U),  P±) 

and 

doo(V,  F ( P± ) ) = da,(U,  F ( P± ) ) , i 
or 

■tv  = ty. 

Thus  U is  T-efficient.  ■ 


APPENDIX  D 

PASCAL  CODE  FOR  THE  PC  BOARD  ASSEMBLY  HEURISTIC 


{S+,R+,D+} 

USES  crt,  dos; 

LABEL  1; 

CONST 

N _prime  =6  ; {*  Number  of  columns  in  Tch.  graph  *} 

M_prime  = 4 ; { Number  of  rows  in  Tch.  graph} 


MaxP  - 10  ; {#  of  different  parts } 


Size=35; 

MaxN=20; 

MaxM=20; 

Length  = 10; 
Width  = 10; 


{ Total  # of  parts  } 

{ number  of  boxes  in  each  row  i.e.  # of  cols  } 

{ number  of  boxes  in  each  column  i.e.  # of  rows} 
{ Requirement:  MaxM  & MaxN  > N'  + M'  -1} 

{ length  of  a box  } 

{ width  of  a box  } 


TYPE 

Circ  = array[l..Size+l]  of  byte; 

WA  = array[l..MaxM,l..MaxN]  of  byte; 

Stack  = array[l..MaxP]  ofWA; 

Patch  Matrix  = array[l..MaxP,O..Size]  ofbyte; 

Point  = RECORD  { For  the  X and  Y coodinates  of  the  EPs  } 

x,y  : real 
END; 

Points  = array  [1..  Size]  of  Point;  { 100  X and  Y values  of  EFs  } 

Points  Array  = array  [1.. MaxP]  of  Points;  { STORES  THE  PARTS  } 

Directs  = array[0..100]  of  real;  { 101  possible  rows  } 

Hundredints  = array[1..100]  of  integer; 

Directions  = RECORD 
W,  E, 

Wstar,  Estar, 

NW, 

NE, 

SW, 

SE, 

U, 

V : Directs; 

IP,  { for  conversion  to  and  from  screen  } 

PI : Hundredints  { locations  and  facility  coordinates  } 

END; 

NType  = string[10]; 

FlagArr  = array  [1 . . M axM  * MaxN  ] ofbyte; 

VAR  occ,  rocc  : STACK; 

NR.  { Number  ofRows  } 

P,  {#  of  part  types} 

N,  M,  iter,PartType,iter2,iter3, 
i,  j,  counter  : integer; 
ch,  yes,  More  : char; 

EffSet  : Directions;  { Contains  the  Efficient  Set  values  } 

LatPt,  ChebLat  : Points  Array;  { Contains  (col, row)  tuples  of  occupancy  points  } 
Location  : Array[l..MaxP*Size,1..3]  of  real; 

Done,  Plots  : boolean; 

Outfile  : text; 
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dummy  : directs; 

CycleLength, 

Sizes,  quantity  : array  [L.MaxP]  of  integer; 

{ counts  the  number  of  parts  for  each  part  type  } 
Cycles  : array  [L.MaxP]  of  Circ; 

Patch  Array  : Patch  Matrix; 

DL  : Circ; 

Flag  : FlagArr; 

dummyarray  : WA; 

TSP  index  : integer; 

Min  X,  { the  minimum  value  of  the  x-coordinates  } 

Max  X,  { the  maximum  value  of  the  x-coordinates  } 

A : real;  { scaling  factor  value  } 

TLength  : real; 


{$1  conn.pas} 


I****************************  FUNCTIONS  ****************************** 


FUNCTION  Max(  A,B:  real):  real; 
BEGIN 

IF  ( A > B ) THEN 
MAX  :=  A 
ELSE 
MAX  :=  B 


END;  { of  Max  } 


FUNCTION  Min(A,B:  real):real; 
BEGIN 

IF  ( A < B ) THEN 
MIN  :=  A 
ELSE 
MIN  :=  B 


END;  { of  Min  } 


FUNCTION  dist(ptl,pt2  : integer) :real; 
VAR  xl,x2,yl,y2  : integer; 


BEGIN 


xl:=(ptl  modN); 

IF  (xl=0)  THEN  BEGIN 
xl:=N; 

y l:=(ptl  div  N); 

END 

ELSEyl:=(ptl  divN)+l; 
x2:=(pt2  modN); 

IF  (x2=0)  THEN  BEGIN 
x2:=N; 

y2:=(pt2  div  N); 

END 

ELSE  y2:=(pt2  div  N)+l; 
dist:=sqrt((xl-x2)*(xl-x2)+(y  l-y2)*(y  l-y2)); 
END; 


FUNCTION  chdist(ptl,pt2  : integer): real; 

VAR  xl,x2,yl,y2  : integer; 

BEGIN 

xl:=(ptl  modNjprime); 

IF  (x  1=0)  THEN  BEGIN 
xl:=N_prime; 
y 1 :=(ptl  div  Njprime); 

END 

ELSEyl:=(ptl  div  Njrime)+1; 

x2:=(pt2  mod  Njrime); 

IF  (x2=0)  THEN  BEGIN 
x2:=N  prime; 
y2:=(pt2  div  N_prime); 

END 

ELSE  y2:=(pt2  div  Njprime)+1; 

chdist:=max(abs(xl-x2),abs(y  l-y2)); 

END; 

{*******************  END  OF  FUNCTION  DECLARATIONS  ******************** 


i**************************  procedures  ******************************* 

PROCEDURE  Initialize; 

VAR  i,j,k  : integer; 

BEGIN 

TSP  index:=0; 

FOR  k:=l  TO  MaxP  DO 
FOR  j:=l  TON  DO 
FOR  i:=l  TO  M DO  Occ[k][i,j]:=0; 


FOR  K:=l  TO  MaxP  DO 
FOR  J:=l  TO  Size  DO  BEGIN 
LatPt[K][J].X:=0; 

LatPt[K][J].Y  :=0; 

C hebLat  [K]  [J]  .X : =0 ; 
ChebLat[K][J].Y:=0; 

END; 

FOR  k:=l  TO  MaxP*Size  DO  BEGIN 


Location[k,l]:=0; 

Location[k,2]:=0; 

Location[k,3]:=0; 

END; 

END;  {**  Procedure  Initialize  **} 


PROCEDURE  Generate_0(pt : integer); 
VAR  row,  col : integer; 
z,w  : real; 

BEGIN 

FOR  iter:=l  TO  quantity[pt]  DO  BEGIN 
z:=random; 
z:=z*N_prime*width; 
w:=random; 
w:=w*M_prime*length; 
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tsp  index:=tsp  index+ 1 ; 

Location[tsp  index,  1 ] : = (int(z/width)+  (width/2))  * w idth ; 
Location[tsp  index, 2]:=(int(w/width)+(width/2))*width; 
Location  [tsp  index,  3 ] : =pt; 

col:=tranc(z/width)  +1; 
row:=trunc(w/length)+ 1 ; 

ChebLat[PT]  [iter]  .x:=col; 

ChebLat  [PT]  [iter]  .y  :=row ; 

rocc  [PT]  [col,row]  :=  1 ; 


{**  tranform  Chebyshev  plane  to  rectilinear  plane  **} 
LatPt[PT][iter].x:=col-row+M_prime; 

LatPt  [PT]  [iter]  .y  :=col+row- 1 ; 

occ[PT]  [col-row+M_prime,col+row- 1 ] := 1 ; 

END; 

END;  {**  Procedure  Generate  0 **} 

^**************>|!*****!(i!(c**s|,****^****  + ^ + ^ + 4:  + + + + + + ^*  + *^**  + + + * + + + + + + + + + + + + ^ + | 


PROCEDURE  Generate  l(pt  : integer); 
VAR  row,  col : integer; 
z,w,xl,yl,x2,y2  : real; 


BEGIN 

writeOenter  coords  of  rectangle  (LL,  UR)  : '); 

readln(yl,xl,y2,x2); 

FOR  iter;=l  TO  quantity[pt]  DO  BEGIN 
z:=random; 
z:=z*(x2-xl)+xl; 
w:=random; 
w:=w*(y2-yl)+yl; 
tsp  index  :=tsp  index+ 1 ; 


Location  [tsp  index,  1 ] : =z; 
Lo  cation  [tsp  index,  2 ] : =w ; 
Location[tsp  index, 3]:=pt; 

col:=trunc(z/width)  +1; 
row  :=trunc(w/length)+ 1 ; 
ChebLat[PT]  [iter]  .x:=col; 
ChebLat[PT]  [iter]  y :=row; 


rocc  [PI]  [col, row]  ;= 1 ; 


{**  tranform  Chebyshev  plane  to  rectilinear  plane  **} 
LatPt[PT][iter].x:=col-row+M_prime; 

LatPt  [PT]  [iter]  y :=col+row- 1 ; 

occ[PT]  [col-row+M_prime,col+row- 1 ] := 1 ; 

END; 

END;  {**  Procedure  Generate  1 **} 


PROCEDURE  Generate  2 (pt : integer); 

VAR  row,  col : integer; 

choice, Z,w,yl,xl,y2,x2,y3,x3,y4,x4  : real; 

BEGIN 

writeCenter  coords  of  rectangle  (LL1,  UR1,  LL2,  UR2) : '); 
readln(yl,xl,y2,x2,y3,x3,y4,x4); 

FOR  iter:=l  TO  quantity  [pt]  DO  BEGIN 
z:=random; 
choice:=random; 

IF  (choice  < 0.5)  THEN  BEGIN 
z:=z*(x2-xl)+xl; 
w:=random; 
w:=w*(y2-yl)+yl; 

END 

ELSE  BEGIN 
z:=z*(x4-x3)+x3; 
w:=random; 
w:=w*(y4-y3)+y3; 

END; 

TSP  index:=tsp  index  + 1 ; 

Location[tsp  index, l]:=z; 

Lo  cation [tsp  index, 2 ] : =w ; 

Location[tsp  index, 3]:=pt; 

col:=trunc(z/width)  +1; 
row:=trunc(w/length)+ 1 ; 

ChebLat[PT]  [iter]  .x:=col; 

ChebLat[PT]  [iter]  .y  :=row; 

rocc  [PT]  [col,row]  :=  1 ; 

{**  tranform  Chebyshev  plane  to  rectilinear  plane  **} 
LatPt[PT][iter].x:=col-row+M_prime; 

LatPt[PT]  [iter]  .y  :=col+row- 1 ; 

occfPT]  [col-row+M  prime, col  t row- 1 ] :=1 ; 

END; 

END;  {**  Procedure  Generate  2 **} 

PROCEDURE  Instance; 

VAR  i,  j,  PT,  code  number  : integer; 

BEGIN 

writeCenter  # of  part  types  (0-',MaxP,') : '); 
readln(P); 

{**  random  generation  **} 

Randomize; 

FOR  PT:=1  TO  P DO  BEGIN 
FOR  i:=l  TO  N prime  DO 
FOR  j:=l  TO  M_prime  DO  rocc[PT][i,j]:=0; 
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writeCenter  # of  parts  for  parttype  \pt,'  (0-',size,') : '); 
readln  (quantity  [pt]) ; 

writef  enter  code  number  for  input  (0-2): '); 
readln(code  number); 

IF  (code  number  = 0)  THEN  Generate  0(pt); 

IF  (code  number  =1)  THEN  Generate  1 (pt); 

IF  (code  number  = 2)  THEN  Generate  2 (pt); 

{if  (code  number  = 3)  THEN  Generate  3 (pt);} 

END;  {of  for} 

FOR  PT:=1  TO  P DO  BEGIN 
Sizes  [PT]:=0; 

FOR  i:=l  TO  M DO 
FOR  j:=l  TON  DO 

IF  (occ[PT][i,j]=l)  THEN  Sizes|TT]:=Sizes[PT]+l; 

FOR  i:=l  TO  sizes[pt]  DO  Cycles[pt][i]:=0; 

END; 

END;  {**  Procedure  Instance  **} 

/*********************************************************************** 


{**********  SORTING  *****  Using  Quick  Sort  *****************************} 

{**  **j 

{**  This  procedure  is  Quick  Sort  which  is  one  of  the  fastest  if  **} 

{**  not  the  fastest  known  sort  in  Pascal.  The  Quick  Sort  is  **} 

{**  recursive  and  of  order  n log(n)  **} 

PROCEDURE  SortData(VAR  EF:  Points;  Start,  Finish  : integer); 

VAR  Temp:  Point; 

StarterValue:  real; 

Left,  Right : integer; 

BEGIN 

Left:=Start; 

Right:=Finish; 

StarterValue:=EF[(Start+Finish)  div  2 ].Y; 

REPEAT 

WHILE  EF[Left].Y  > StarterValue  DO 
Left:=  Left  + 1; 

WHILE  StarterValue  > EF[Right].Y  DO 
Right  :=  Right  - 1; 

IF  Left  <=  Right  THEN  BEGIN 
Temp  :=  EF[Left]; 

EF[Left]  :=  EFfRight]; 

EF  [Right]  :=  Temp; 

Left  :=  Left  + 1; 

Right  :=  Right  - 1 
END 

UNTIL  Right  <=  Left; 

IF  Start  < Right  THEN  SortData  (EF,Start, Right); 

IF  Left  < Finish  THEN  SortData  (EF,Left,Finish) 

END;  {**  Procedure  SortData  **} 

I*********************************************************************** 
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j**********  efficient  SET  CALCULATIONS  ********************************} 
{**  **} 

{**  This  procedure  performs  Steps  1 and  2 of  the  Row  Algorithm  as  **} 

{**  presented  in  the  paper  by  Chalmet,  Francis  and  Kolen.  It  **} 

{**  finds  the  values  of  W,  E,  NW,  NE,  SW,  SE,  E*.  and  W*.  **} 

PROCEDURE  FindSet(EF:Points;N integer;  VAR  EffSet:  Directions;VARNR:integer); 
VAR  I.  P:  integer; 

BEGIN 

WITH  EffSet  DO  BEGIN 
{*  Finding  of  Wi 's  & Ei 's  *} 

P:=l; 

W[1]:=EF[1].X; 

E[1]:=EF[1].X; 

IP[1]  :=  1; 

PI[1]  :=  1; 

FOR  I:=2  TO  N DO  BEGIN 
IF  ( EF[I].Y  = EF[I-1].Y  ) THEN  BEGIN 
IP[I]  :=  P; 

PI[P]  :=I; 

IF  (EFPJ.X  < W[P])  THEN  W[P]  :=  EFp].X; 

IF  (EFPJ.X  > E[P])  THEN  E[P]  :=  EFp].X 
END 

ELSE  BEGIN 
P :=P+  1; 

IPP]  :=  P; 

PIP>]  :=  I; 

Wp>]:=  EFp].X; 

E[P]:=  EF[I].X 
END 
END; 

NR:=  P - 1; 

{*  Finding  NW  & NE's  *} 

SW[1]:=  W[l]; 

SE[1]:=  E[l]; 

FOR  P:=2  TO  NR  DO  BEGIN 
SW[P]  :=  Min  ( wp»],  SWp>- 1] ); 

SE[P]  :=  Max  ( E[P],  SE[P-1] ); 

END; 

{*  Finding  SW  & SE's  *} 

NWp\[R]:=W[NR+l]; 

NE  P'JR] : =Ep>JR+ 1 ] ; 

FOR  P:=NR-1  DOWNTO  1 DO  BEGIN 
NW[P]  :=  Min  ( Wp>+1],  NW[P+1] ); 

NEP>]  :=  Max  ( E[P+1],  NE[P+1] ) 

END; 

{*  Finding  W*  and  E*  *} 

FOR  P:=l  TO  NR  DO  BEGIN 
WstarP1]  :=  Max  ( SWp1],  NW[P] ); 

EstarP1]  :=  Min  ( SEp>],  NE[P] ); 

END; 

{*  Finding  U and  V *} 

NW[0]:=  -l*maxint; 


{ - infinity  } 
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NE[0]:=  maxint;  { infinity  } 

SW[NR+1]  := -l*maxint;  {-infinity} 

SE[NR+1]  :=  maxint;  { infinity  } 

FOR  P:=l  TO  NR+1  DO  BEGIN 
U[P]  :=  MAX  ( SW[P],  NW[P-1] ); 

V[P]  :=  MIN  ( SE[P],  NE[P-1] ) 

END 

END  {of  with  } 

END;  {**  Procedure  FindSet  **} 

r *********************************************************************** t 


/**********  Resujts  *****  prints  the  Efficient  Set  results  **************} 

{**  **j 

{**  This  long  routine  writes  out  the  numerical  results  to  the  screen  **} 

{**  and  to  a file  if  desired.  **} 


PROCEDURE  Results (EffSet directions;  N,  NR:  integer); 
VAR  x,P  : integer; 

XI,  X2,  Yl,  Y2  : real; 

Yes:  char; 

Name : Ntype; 


BEGIN 

WITH  EffSet  DO  BEGIN 
FOR  x:=l  TO  Nr  DO  dummy[x]:=wstar[NR-x+l]; 
FOR  x:=l  TO  Nr  DO  Wstar[x]:=dummy[x]; 

FOR  x:=l  TO  Nr  DO  dummy [x]:=estar[NR-x+l]; 
FOR  x:=l  TO  Nr  DO  estar[x]:=dummy[x]; 

END; 

END;  {**  Procedure  Results  } 


PROCEDURE  Boundary  (k:integer); 

VAR 

counter,elt,inc,x,i,j,l : integer; 

Leftmost,TOE  : longint; 

rightmost, BOE  : longint; 

row  : array  [L.MaxM]  of  integer; 

BEGIN 

WITH  EffSet  DO  BEGIN 
FOR  x:=l  TO  Nr  DO  dummy  [x]:=wstar[NR-x+l]; 
FOR  x:=l  TO  Nr  DO  Wstar[x]:=dummy[x]; 

FOR  x:=l  TO  Nr  DO  dummy [x]:=estar[NR-x+l]; 

FOR  x:=l  TO  Nr  DO  estar[x]:=dummy[x]; 

END; 

{*  initialize  the  fist  of  circuits:  cycles  *} 

FORj:=l  TO  Size+1  DO  Cycles[k][j]:=0; 

inc:=0;  { # of  nonzero  rows  in  occ  } 

FOR  j:=l  TO  M DO  BEGIN 
counter:=0; 

FOR  1:=1  TO  N DO  counter:=counter+occ[k][j,l]; 
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IF  (counter  <>  0)  THEN  BEGIN 
inc:=inc+l; 

row[inc]:=j;  {:  stores  inc-th  nonzero  row  of  occ} 

END; 

END; 

FOR  j:=l  TO  Size+1  DO  DL[j]:=0; 

FOR  j:=l  TO  N*M  DO  Flag[j]-0; 

elt:=0;  { row  1 of  boundary  of  Efset  } 

j:=0; 

REPEAT  j:=j+l; 

UNTIL  (occ[k][row[l],j]  = 1); 

Leftmost:=j; 

TOE— trunc(Effset.Wstar[l]); 

FOR  j -Leftmost  TO  TOE  DO  BEGIN 
IF  (occ  [k]  [row  [ 1 ] , j] = 1 ) TFIEN  BEGIN 
elt:=elt+l; 

DL[elt]  — j+N*  (row  [ 1 ]- 1 ); 

Flag[j+N*(row[l]-l)]— 1; 

END; 

END; 

FOR  i:=2  TO  inc-1  DO  BEGIN  {i=row  number  in  Francis's  Efset  alg} 

j:=0; 

REPEAT  j:=j+l; 

UNTIL  (occ[k][row[i],j]  = 1); 

Leftmost— j; 

TOE:=trunc(Effset.wstar[i]); 

BOE : =trunc  (Effset.  wstar  [i- 1 ] ); 

IF  (BOE  >=  TOE)  THEN  FOR j -BOE  DOWNTO  Leftmost  DO 
IF  (occ[k][row[i],j]=l)  THEN  BEGIN 
elt— elt+1; 

DL[elt]:=j  + N*(row[i]-l); 

Flag[j+N*(row  [i]-l)]:=l; 

END; 

IF  (BOE  < TOE)  THEN  BEGIN 
FOR  j -leftmost  TO  TOE  DO 
IF  (occ[k][row[i],j]=l)  THEN  BEGIN 
elt— elt+1; 

DL[elt]:= j+N*(row[i]-l); 

Flag  [j+N  * (row  [i]- 1 )]— 1 ; 

END; 

END; 

END  {ofi}; 

{*  last  row  of  Efset  *} 

j-0; 

REPEAT  j-j+1; 

UNTIL  (occ[k][row[inc],j]  = 1); 

Leftmost— j; 
j— N+l; 

REPEAT  j:=j-l; 

UNTIL  (occ[k][row[inc],j]  = 1); 
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rightmost:  =j; 

FOR  j:=leftmost  TO  rightmost  DO 
IF  (occ[k][row[inc],j]=l)  THEN  BEGIN 
elt:=elt+l; 

DL  [elt] : =j+N*  (row  [inc]- 1 ) ; 

Flag  [j+N*  (row[inc]- 1 )]  := 1 ; 

END; 

{*  DL  is  expanded  to  include  the  vertices  of  UL  *} 

FOR  i : =inc- 1 DOWNTO  2 DO  BEGIN  {irow  number  in  Francis's  Efset  alg} 
j:=N+l; 

REPEAT  j:=j-l; 

UNTIL  (occ[k][row[i],j]  = 1); 
rightmost:=j; 

TOE : =trunc  (Effset.  estar  [i]) ; 

BOE:=trunc(Effset.estar[i-l]); 

IF  (TOE  >=  BOE)  THEN 
FOR  j -rightmost  DOWNTO  BOE  DO 

IF  ((occ[k][row[i],j]=l)  and  (Flag[j+N*(row[i]-l)]=0))  THEN  BEGIN 
elt:=elt+l; 

DL[elt]:=j  + N*(row[i]-l); 

Flag[j+N*  (row[i]- 1 )]  :=1 ; 

END; 

IF  (BOE  > TOE)  THEN 
FOR  j:=TOE  TO  rightmost  DO 

IF  ((occ[k][row[i],j]=l)  and  (Flag[j+N*(row[i]-l)]=0))  THEN  BEGIN 
elt:=elt+l; 

DL[elt]:=j+N*(row[i]-l); 

Flag  [j+N*  (row[i]- 1 )]  :=1 ; 

END; 

END  { for  i} ; 

{*  1st  row  of  Effset  on  the  way  UP  *} 

j:=0; 

REPEAT  j:=j+l; 

UNTIL  (occ[k][row[l],j]  = 1); 

Leftmost:=j; 

j:=N+l; 

REPEAT  j:=j-l; 

UNTIL  (occ[k][row[l],j]  = 1); 
rightmost:  =j; 

FOR  j:=rightmost  DOWNTO  leftmost  DO 

IF  ((occ [k] [row [ 1 ],j]=  1 ) and  (Flag[j+N*(row[l]-l)]=0))  THEN  BEGIN 
elt:=elt+l; 

DL[elt]:=j+N*(row[l]-l); 

Flag[j+N*(row[l]-l)]:=l; 

END; 

DL[elt+l]:=DL[l]; 

FOR  j:=l  TO  quantity [parttypel+1  DO  Cycles[k][j]:=DL[j]; 

END;  {**  Procedure  Boundary  **} 

|*********************************.,.*.(.#********!|!*********************#.)c.(.,|cJ 
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PROCEDURE  WAI(k:integer); 

VAR 

maxangle,angle,dl,d2,d3  : real; 

insindex,mdex,inspos,vtx,j,ptl,pt2,i,unflaged  vtx:  integer; 
tempcycle : Circ; 

BEGIN 

unflaged  vtx:=0; 

FOR  vtx:=l  TO  M*N  DO  BEGIN 
j:=(vtx  mod  N); 

IF  0=0)  THEN  BEGIN 

j:=N; 

i:=(vtx  div  N); 

END 

ELSE  i:=(vtx  div  N)+l; 

IF  (Flag[vtx]  = 0)  and  (occ[k][i,j]=l)  THEN  unflaged  vtx:  =unflaged  vtx+1; 

END; 

FORvtx:=l  TO  unflaged  vtx  DO  BEGIN 
maxangle:=2; 

inspos:=0;  { position  where  unflaged  part  is  to  be  inserted  } 

FOR  index:=l  TO  M*N  DO  BEGIN 
j:=(indexmodN); 

IF  0=0)  THEN  BEGIN 

j:=N; 

i:=0ndex  div  N); 

END 

ELSE  i:=(index  div  N)+l; 

IF  (Flagfindex]  = 0)  and  (occ[lc][i,j]=l)  THEN  BEGIN 

j:=i; 

REPEAT 

j:=j+i; 

ptl:=Cycles[k][j-l]; 
pt2:=Cycles[k][j]; 
d 1 : =dist(index,pt  1 ) ; 
d2:=dist(index,pt2); 
d3  :=dist(ptl,pt2); 

angle:=(dl*dl+d2*d2-d3*d3)/(2*dl*d2); 

{*  maximum  angle  corresponds  to  minimum  cosine  *} 

IF  (angle  < maxangle)  THEN  BEGIN 
inspos:=j; 
maxangle:=angle; 
insindex:=index; 

END; 

UNTIL  ((Cycles[k][j+1]=0)  or  0=quantity[parttype]+l)); 

END;  {**  of  if  statement  **} 

END;  {**  for  index  **} 

{*  insert  a PART  into  current  CIRCUIT  *} 

F lag  [ins  index]  := 1 ; 

FOR  j:=l  TO  quantity [parttype]+l  DO  tempcycle[j]:=Cycles[k][j]; 

Cycles  [k][inspos]:=insindex; 

FOR  j:=inspos+l  TO  quantity [parttype]+l  DO  Cycles[k][j]:=tempcycle[j-1]; 

END;  {of  vtx} 

END;  {**  Procedure  WAI  **} 

s***************************************** ****************************** 
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PROCEDURE  OPT2(VAR  S:Circ); 

VAR  i,j, counter, impcounter, best  i,best  J,ptl,pt2,pt3.pt4  : integer; 
curDelta,bestDelta,dl,d2,d3,d4  : real; 

Temp  S : Circ; 


BEGIN 


FORi:=l  TO  sizes  [parttype]  DO  Temp  S[i]:=S[i]; 
impcounter:=0; 

REPEAT 

impcounter:  =impcounter+ 1 ; 

counter:=0;  { shows  whether  improvement  is  made  in  current  iteration  } 
bestDelta:=0;  { this  is  an  upper  bound  for  cost  difference  } 
best  i:=0; 
besfj:=0; 

FOR  i:=l  TO  sizes  [parttype]  DO  BEGIN 
FOR  j:=l  TO  sizes  [parttype]  DO  BEGIN 
IF  (i<>j)  THEN  BEGIN 
{ evaluate  cost  of  i,j  interchange  } 
ptl:=temp_S[i]; 

IF  (i=sizes[parttype])  THEN  pt2:=temp  S[l]  ELSE  pt2:=temp  S[i+1]; 
pt3:=temp_S[j]; 

IF  (j=sizes [parttype])  THEN  pt4:=temp  S[l]  ELSE  pt4:=temp  S[j+1]; 

dl:=dist(ptl,pt3); 

d2:=dist(pt2,pt4); 

d3  :=dist(ptl  ,pt2); 

d4:=dist(pt3,pt4); 

curDelta:=dl+d2-d3-d4; 

{*  if  improvement  is  made,  store  i,j,bestDelta  *} 

IF  (curDelta  <=  bestDelta)  THEN  BEGIN 
counter:=counter+ 1 ; 
bestDelta:=curDelta; 
bestj;=j; 
best  i:=i; 

{*  interchange  best  i with  best  j *} 

{*  and  reverse  the  in  between  segment  *} 

IF  (counter  o 0)  THEN  BEGIN 
IF  (best  i <=  best_j)  THEN 

FORi:=best  i+1  TO  bestj  DO  Temp  S[i]:=S[bestJ+best  i+l-i] 
ELSE 


FOR  i;=best j+ 1 TO  best  i DO  Temp  S[i]:=S[besfj+best  i+l-i]; 

FOR  i:=l  TO  sizes  [parttype]  DO  S[i]:=Temp  S[i]; 

END;  {*  if  counter  *} 

END;  {*  if  CurDelta  *} 

END;  {*  if  i *} 

END;  {*  for  j *} 

END;  {*  for  i *} 

UNTIL  ((counter=0)  or  (impcounter=5)); 

S [sizes  [parttype]+ 1 ] :=S  [ 1 ] ; 

END;  {**  Procedure  Opt2  **} 

|***********************************************************************\ 
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PROCEDURE  Merge(partfrom,ptype,row  : integer); 
VAR  i,j,k  : integer; 


BEGIN 


j-=0; 

FOR  i:=partfrom+l  TO  sizes[ptype]  DO  BEGIN 

j:=j+i; 

Patch  Array  [row,j]  :=Cy cles  [ptype]  [i]; 

END; 


FOR  k:=l  TO  partfrom  DO  BEGIN 

j:=j+l; 

Patch  Array  [row]  [j] : =Cy  cles  [ptype]  [k] ; 


END; 

END;  {**  Procedure  Merge  **} 

^***********************************************************************j 


PROCEDURE  Rmerge(partfrom, ptype, row  : integer); 
VAR  i,j,k  : integer; 


BEGIN 

j:=0; 

FOR  i:=partfrom-l  DOWNTO  1 DO  BEGIN 

j:=j+i; 

Patch  Array  [row, j]:=Cycles[ptype][i]; 

END; 

FOR  k:=sizes[ptype]  DOWNTO  partfrom  DO  BEGIN 

j'=j+i; 

Patch  Array  [row]  [j] : =Cy  cles  [ptype]  [k] ; 

END; 

END;  {**  Procedure  Rmerge  **} 

!***********************************************************************! 

PROCEDURE  Translate; 

VAR 

j,  yl,  xl,  a,  b : integer; 

BEGIN 

FOR  j:=l  TO  sizes [P ARTTYPE]  DO  BEGIN 

{*  translate  sequence  number  to  rect  lattice  coords  *} 
y 1 :=(Cycles[parttype][j]  mod  N); 

IF  (y  1=0)  THEN  BEGIN 
yl:=N; 

xl:=(Cycles[parttype][j]  div  N); 

END 

ELSE  xl:=(Cycles[parttype][j]  div  N)+l; 

{*  translate  rect  coords  to  Cheb  coords  *} 
a:=round((x  1 +y  1 -M_prime+ 1 )/2); 
b :=round((x  1 -y  1 +M_prime+ 1)/2); 

Cy cles  [parity pe]  [j]:=(b- 1 )*N_primef  a; 

END;  {*forj*} 

END;  {**  Procedure  Translate  **} 

|***********************************************************************j 


PROCEDURE  Patch; 

VAR 

ptype,  partffom,  partto,  typeto,  cur  typeto, 
vtxl,  vtx2,  xl,  yl,  a,  b,  i,  j,  k : integer; 

distance,  mindist  : real; 

Flag  : array  [l.MaxP]  of  byte; 

BEGIN 

FOR  i:=l  TO  P DO 

FORj:=0  TO  Size  DO  Patch  Array[i,j]:=0; 

FOR  i:=l  TO  P DO  Flag[i]:=0; 

{*  find  minimum  distance  from  Cycle[l]  to  other  cycles  *} 
mindist:  =maxint; 

FORptype:=2  TO  P DO 
FOR  j:=l  TO  sizes[l]  DO 
FOR  i:=l  TO  sizes(ptype]  DO  BEGIN 
distance:=chdist(Cycles[ptype][i],Cycles[l][j]); 

IF  (distance  < mindist)  THEN  BEGIN 
mindist : =dis  tance; 
typeto  :=ptype; 
partfrom:=j; 
partto  :=i; 

END; 

END; 

Flag[l]:=l; 

Flag[typeto]:=l; 

Patch  Array[l,0]:=l;  { 0-th  entry  of  each  row,  stores  part  type  } 

k:=2; 

IF  (parttype  > 1)  THEN  Patch  Array  [k,0]:=typeto; 

merge(partffom,l,l);  { copy  part  1 in  row  1 of  patch  array  } 

FOR  k:=2  TO  P-1  DO  BEGIN  { k is  the  current  row  of  patch  array  } 
mindist:  =maxint; 

{*  identify  the  two  possible  vertices  vtxl,vtx2  where  the  next  row  could  end 
IF  (sizes[typeto]  <=  2)  THEN  vtx2:=0; 

IF  (sizes[typeto]  = 1)  THEN  vtxl:=partto; 

IF  (sizes [typeto]  = 2)  THEN  vtxl:=(partto mod 2)+l; 

IF  (sizes [typeto]  > 2)  THEN  BEGIN 
vtxl:=partto-l; 
vtx2:=partto+l; 

IF  (vtxl=0)  THEN  vtxl:=sizes[fypeto]; 

IF  (vtx2=sizes[partto]+l)  THEN  vtx2:=l; 

END; 

cur  typeto  :=fypeto; 

FOR  ptype:=l  TO  P DO 
IF  (Flag [ptype]  = 0)  THEN  BEGIN 

{*  find  closest  neighbor  of  vtxl  which  is  not  in  the  tour  yet  *} 

FOR  i:=l  TO  sizes[ptype]  DO  BEGIN 
distance:=chdist(Cycles [ptype][i],Cycles[cur  typeto] [vtxl]); 

IF  (distance  < mindist)  THEN  BEGIN 
mindist:  =distance; 
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typeto:=ptype; 

partfrom:=vtxl; 

partto:=i; 

END; 

END;  {**  for  i **} 

{*  find  closest  neighbor  of  vtx2  which  is  not  in  the  tour  yet  *} 

IF  (vtx2<>0)  THEN 
FOR  i:=l  TO  sizes[ptype]  DO  BEGIN 
distance:=chdist(Cycles[ptype][i],Cycles[cur  typeto][vtx2]); 

IF  (distance  < mindist)  THEN  BEGIN 
mindist:  =distance; 
typeto:=ptype; 
partlfom:=vtx2; 
partto:=i; 

END; 

END;  {**  for  i **} 

END;  {**  of  if  flag[ptype]=0  **} 

Flag[typeto]:=l; 

Patch  Array  [k+ 1 ,0]  :=typeto; 

IF  (partfrom=vtxl)  THEN  merge(partfrom,cur typeto,k); 

IF  ((vtx2<>0)  and  (partfrom=vtx2))  THEN  rmerge(partfrom,cur  typeto,k); 

END;  { of  k } 


{*  copy  the  last  tour,  in  row  P of  patch  array  *} 
merge(partto- 1 ,typeto,P); 

END;  {**  Procedure  Patch  **} 


{***********************************^**********************************1 
PROCEDURE  compute  length;  {used  in  Procedure  Write  Sequences} 

VAR  i,  j,  pt , sum  : integer; 


BEGIN 

tlength:=0; 

FOR  i:=l  TO  P DO 

FOR  j:=l  TO  sizes  [Patch  array  [i,0]]-l  DO 
tlength:=tlength+chdist(Patch  array [i,j], Patch  array [i,j+ 1]); 

FOR  i:=l  TO  P-1  DO 

tlength:=tlength+max(l,chdist(Patch  array [i,sizes [Patch  array [i,0]]], 
Patch  array  [i+ 1 , 1 ] )) ; 

tlength:=tlength+max(l,  chdist(Patch  array [P,sizes[Patch  array [P,0]]], 
Patch  array[l,l])); 

FORpt:=l  TO  P DO  tlength:=(tlength+(quantity[pt]-sizes[pt])); 
tlength:=width*tlength; 


writeln; 

writelnCTotal  Length  = ',tlength:6:0); 


sum:=0; 
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FOR  i:=  1 TO  P DO  sum:=sum+quantity[i]; 
writelnCLower  bound  = width*sum:6); 

writelnCMaximum  error  = ',(tlength-width*sum)/(width*sum)*  100:6:2); 

END; 

PROCEDURE  TSP  output; 

VAR 

Name  : string[8]; 

i,  j,  sum,  tsp  dist : integer; 

BEGIN 

writeln; 

WriteC  Do  you  want  to  create  an  output  file  ? (Y /N)  '); 

Readln(Y  es); 

IF  (( YES  = ’Y')  or  (YES  = 'y'))  THEN  BEGIN 
{*  start  writing  to  file  *} 

Writeln; 

WriteC  Enter  output  file  name  : '); 

Readln(Name); 

Assign(OutFile,Name); 

Rewrite(OutFile) ; 

sum:=0; 

FOR  i:=l  TO  P DO  sum:=sum+quantity[i]; 
writeln(outfile,sum); 

FOR  i:=2  TO  sum  DO  BEGIN 
FOR  j:=l  TO  i-1  DO  BEGIN 

tsp_dist:=trunc(max(abs(location[i,l]-location[j,l]),abs(location[i,2]-location[J,2]))); 
IF  (tspdist  < width)  THEN  tsp_dist:=width; 

IF  (location[i,3]  o location[j,3]) 

THEN  tsp_dist:=tsp_dist+2*round(max(M_prime,N_prime)); 
write(outfile,' tsp  dist); 

END; 

writeln  (outfile); 

END; 

writeln(outfile,2*round(max(M_prime,N_prime))*(P+l)+tlength); 

close(outfile); 

writeln C clustered  TSP  bound:  ',2*round(max(M_prime,Njrime))*(P+l)+tlength); 
END 

ELSE  i:=i; 

END;  {**  Procedure  TSP  Output  **} 

|*******iM!**********^****  + + + + #**^  + + ^ + + + * + <.*  + + + + + + + + + + + ^<.#**  + + ^ + + + ^ + ,.  + + + + j 

PROCEDURE  WriteSequences; 

VAR  i,j,  k : integer; 

BEGIN 

FOR  i:=l  TO  P DO  BEGIN 
FOR  j:=l  TO  Mjrime  DO  BEGIN 
FORk:=l  TO  N_prime  DO  write(rocc[i][k,j]:3); 
writeln; 

END; 

writeln; 
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IF  (P>3)  THEN 

IF  ((i  mod  3)=0)  THEN  BEGIN 
writeln; 

wnte('Press  <Enter>  TO  continue.'); 

writeln; 

ch:=readkey; 

clrscr; 

END;  {*  if  *} 

END;  {*  for  *} 

FOR  i:=l  TO  P DO  BEGIN 


FOR  j:-0  TO  sizes[Patch  Array[i,0]]  DO  write(Patch  Array[i,j]:3); 
writeln; 

END; 

compute  length; 

TSPOutput; 

END;  {**  Procedure  WriteSequences  **} 

z*********************************************^************************ 


/*********************************************************************** 
PROCEDURE  FSET  SEQ; 


} 

} 


BEGIN 

SortData(LatPt[PartType].l,quantity[parttype]); 

F indSet  (LatPt[PartTy  pe] , quantity  [partty  pe]  ,EffSet,NR); 

Results  (EffSet,quantity  [parttype],NR); 

Boundary  (PartType); 

WAI(PartType); 

OPT2(Cycles[PartType]); 

Translate; 

END;  {**  Procedure  FSet  Seq  **} 

I*********************************************************************** 

PROCEDURE  HAMILTON; 

VAR 

j,counter,k,next,previous,current,col,row,inspos,vtx,i, outer  : integer; 
flag,tempcycle : FlagArr; 

PROCEDURE  Expand(Previous, Current : integer;  VAR  Next : integer); 

VAR 

col,  row,  k,  i,  vertex  : integer, 
clock  : array[0..7]  of  integer; 

BEGIN 

clock[0]:=current-N_prime; 
c lo  ck  [ 1 ] : =curr  ent-N_prime+ 1 ; 
clock  [2] : =current+ 1 ; 
clock[3]:=current+N_prime+l ; 
clo  ck  [4] : -currcnt+N  _prime; 
clockt5]:=current+N  jrime- 1 ; 
clo  ck  [6] : =current- 1 ; 
clock[7]  :=current-N_prime- 1 ; 
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IF  (current  <=  N jrime)  THEN  BEGIN 
clock[7]:=maxint; 
clock[0]:=maxint; 
clock[l]:=maxint; 

END; 

IF  (current  > (M  prime-l)*N  prime)  THEN  BEGIN 
clock[3]:=maxint; 
clock  [4] : =maxint; 
clock[5]:=maxint; 

END; 

IF  ((current  mod  N_prime)  = 0)  THEN  BEGIN 
clock[l]:=maxint; 
clock[2]:=maxint; 
clock[3]:=maxint; 

END; 

IF  ((current  mod  N_prime)  = 1)  THEN  BEGIN 
clock[5]:=maxint; 
clock[6]:=maxint; 
clock[7]:=maxint; 

END; 

IF  counter=l  THEN  i:=0 
ELSE  BEGIN 
i:=-l; 

REPEAT 

i:=i+l; 

UNTIL  (previous=clock[i]); 

END; 

{*  the  following  piece  searches  clockwise  for  next  starting  from  *} 

{*  previous  direction  *} 

k:=0; 

REPEAT 

k:=k+l; 

vertex:=clock[(i+k)  mod  8]; 
col:=(vertex  mod  N_prime); 

IF  (col=0)  THEN  BEGIN 
col:=N_prime; 
row:=(vertex  div  N_prime); 

END 

ELSE  row:=(vertex  div  N pnme)+ 1 ; 

UNTIL  ((vertex  < maxint)  and  (flagfvertex]  = 0)  and 
(rocc[parttype][col,row]=l)); 
next:=vertex; 

END; 

{-+-+-+-+-H — K-+-H — h-H — (--+-+-+-H — h-+-H — I--+-+-+-+-H — h-H — I--+-H — I--+-+-H — (--+-} 
BEGIN  {*  Procedure  Hamilton  *} 

{*  find  first  1 on  rocc  matrix;  the  position  of  that  1 is  (row,col)  *} 

j:=0; 

REPEAT 

j:=j+i; 

col:=(j  mod  N_prime); 

IF  (col=0)  THEN  BEGIN 
col:=N_prime; 


row:=(j  div  N_prime); 

END 

ELSE  row:=(j  div  N_prime)+1; 

UNTIL  (rocc[parttype][row,col]=l); 
counter  :=1; 

Cycles  [partty  pe]  [counter]  :=j ; 
previous : =maxint ; 
current:=j; 

{*  expand  Hamilton  cycle  *} 

FOR  k:=l  TO  M_prime*N_prime  DO  flag[k]:=0; 

REPEAT 

expand(previous,current,next); 
previous : =current; 
current:  =next; 
counter: =counter+ 1 ; 

Cycles  [partty  pe]  [counter] : =next; 
flag[next]:=l; 

UNTIL  (next=Cycles[parttype][l]); 

{*  insert  remaining  unllaged  vertices  *} 

FOR  i:=2  TO  Mjrime-1  DO 
FOR  j:=l  TO  Njrime  DO  BEGIN 
vtx:=(i-l)*N_prime+j; 

IF  ((flag[vtx]=0)  and  (rocc[parttype][j,i]=l))  THEN  BEGIN 
flag[vtx]:=l; 

FORk:=l  TO  Sizes  [partty  pe]  DO  tempcycle[k]:=Cycles[parttype][k]; 

IF  (i=2)  THEN  outer:=vtx-N_prime; 

IF  (i=3)  THEN  outer:-vtx+N  prime; 
inspos:=0; 

REPEAT 

inspos:=inspos+l; 

UNTIL  (Cycles  [parttype][inspos]=outer); 

Cycles  [parttype]  [inspos+ 1 ] :=vtx; 

FORk:=inspos+2  TO  Sizes[parttype]  DO  Cycles  [parttype]  [k]:=tempcycle[k-l]; 
END; 

END; 

END;  {**  Procedure  Hamilton  **} 

I*********************************************************************** 

BEGIN  { Main  Program  } 

1: 

N:=N_prime+M  jrime- 1 ; 

Mr  N prime  i M prime- 1; 
clrscr; 

Initialize; 

Instance; 

clrscr; 

FOR  PartType:=l  TO  P DO  BEGIN 
FOR  i:=l  TO  MaxN  DO 

FOR  j:=l  TO  MaxM  DO  dummyarray[i,j]:=occ[PartType][i,j]; 

FOR  i:=l  TO  MaxN  DO 

FOR  j:=l  TO  MaxM  DO  occ[PartType][i,j]:=dummyarray[j,i]; 

Biconnected(PartType, Sizes  [PartType],Status); 
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IF  Status  THEN  writelnCOccupancy  graph  is  Hamiltonian.') 

ELSE  writelnCOccupancy  graph  is  not  Hamiltonian.'); 

IF  Status  THEN  {Hamilton}  FSet  Seq  ELSE  Fset  Seq; 

END; 

IF  (P>=32)  THEN  BEGIN 
writeCPress  <Enter>  to  continue.'); 
writeln; 
ch:=readkey; 
clrscr; 

END; 

Patch; 

writelnCEnter  any  key  to  continue.'); 
readln; 

W riteSequences ; 

{ writeCPress  <Enter>  to  continue.'); 
writeln; 
ch:=readkey; 

} 

{ writeC  Do  you  want  to  run  another  PC  board?  (Y /N) '); 

Readln(Y  es); 

IF  ((YES='Y)  or  (YES-y'))  THEN  GOTO  1 

ELSE  writeCPress  <Enter>  to  terminate  the  program.'); 


writeln;} 

ch:=readkey; 

END.  { of  Main  Program  }• 


var  Adj  : WA; 

MaxV  : integer; 
status  : Boolean; 

PROCEDURE  DFS(A  : WA;  MaxV  : integer;  root : integer;  var  status  : Boolean); 
var  id,  k,  preroot : integer; 
val  : array[l..MaxM*MaxN]  of  integer; 

PROCEDURE  VISIT(k  : integer); 
var  t : integer; 
begin 
id:=id+ 1 ; 
val[k]:=id; 
for  t:=l  to  MaxV  do 
if  (A[k,t]=l)  then 
If  (val[t]=0)  then  visit(t); 
end;  {**  Procedure  Visit  **} 

BEGIN 

id:=0; 

for  k:=l  to  MaxV  do  val{k]:=0; 
visit(root); 
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if  (root=l)  then  preroot:=MaxV  else  preroot:=root-l; 
for  k:=l  to  MaxV  do 

if  ((val[k]=0)  and  (kopreroot))  then  status:=false; 

END;  {**  Procedure  DFS  **} 

PROCEDURE  TWO  CONNECTED  (A  : WA;MaxV  : integer;var  STATUS  : Boolean); 
var  vtx,  i,  j,  root : integer; 

Adj  : WA; 

BEGIN 
status  :=true; 
vtx:=0; 
repeat 
vtx:=vtx+l; 
for  i:=l  to  MaxV  do 
for  j:=l  to  MaxV  do 
begin 

Adj  [ij] : = A[i  j] ; 

if  ((i=vtx)  or  (j=vtx))  then  Adj[i,j]:=0; 
end; 

if  (vtx<MaxV)  then  root:  vtx+l  else  root:=l; 

DFS(Adj, MaxV, root, status); 
until  ((status=false)  or  (vtx=MaxV)); 

END;  {**  Procedure  TwoConnect  **} 


PROCEDURE  BICONNECTED  (PARTTYPE  : BYTE;  MaxV  : INTEGER;  VAR  STATUS  : BOOLEAN); 
{ determines  biconnectivity  of  occupancy  graph  for  parttype 
STATU S=T  means  biconnected 
STATUS=F  means  disconnected  or  simply  connected  } 

var 

AdjMatrix  : WA; 

Vertexlist : array[l..MaxM*MaxN]  of  integer; 
i,j,elt : integer; 


BEGIN 

FOR  i:=l  TO  MaxV  DO 
FOR  j:=l  TO  MaxV  DO 
AdjMatrix[i,j]  :=0; 
elt:=0; 

FOR  i:=l  TO  M_prime  DO 
FOR  j:=l  TO  N_prime  DO 
if  (rocc[parttype][j,i]=l)  then 
begin 
elt:=elt+l; 

v ertexlis t [elt]  :=(i- 1 )*  N prime*  j ; 
end; 

FOR  i:=l  TO  MaxV  DO 
FOR  j:=i+l  TO  MaxV  DO 

IF  ( ((vertexlist[j]=vertexlist[i]+l)  AND  (vertexlist[i]  mod  N prime<>0)) 

OR  ((vertexlist[j]=vertexlist[i]+N_prime-l)  AND  (vertexlistfi]  mod  N jirimeol)) 
OR  ((vertexlist[j]=vertexlist[i]+N  _prime)) 
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OR  ((vertexlist[j]=vertexlist[i]+N_prime+l)  AND  (vertexlist[i]  mod  N prime<>0))) 
THEN  BEGIN 
AdjMatrix[i,j]:=l; 

AdjMatrix[j,i]:=l; 

END; 

TWO  CONNECTED(AdjMatrix,MaxV, STATUS); 

END; 


REFERENCES 


REFERENCES 


Ahmadi,  J.,  S.  Grot zinger  and  D.  Johnson  [1988],  "Component 
Allocation  and  Partitioning  for  a Dual  Delivery  Placement 
Machine,"  Operations  Research,  vol.  36,  no.  2,  pp.  176-191. 

Allison,  D.C.S.  and  M.T.  Noga  [1984],  "The  Li  Traveling 
Salesman  Problem,"  Information  Processing  Letters,  vol.  18, 
pp.  195-199. 

Ball,  M.  0.  and  M.  J.  Magazine  [1988],  "Sequencing  of 
Insertions  in  Printed  Circuit  Board  Assembly,"  Operations 
Research,  vol.  36,  pp.  192-201. 

Barachet,  L.  L.  [1957],  "Graphic  Solution  of  the  Traveling 
Salesman  Problem,"  Operations  Research,  vol.  5,  pp.  841- 
845. 

Bard,  J.  F.,  R.  W.  Clayton  and  T.  A.  Feo  [1989],  "Optimizing 
Machine  Setup  and  Component  Insertion  in  Printed  Circuit 
Board  Assembly,"  working  paper.  Operations  Research  Group, 
College  of  Engineering,  University  of  Texas,  Austin,  TX 
78712. 

Bartholdi,  J.J.  and  L.K.  Platzman  [1982],  "An  0(N  log  N) 
Planar  Travelling  Salesman  Heuristic  Based  on  Spacefilling 
Curves,"  Operations  Research  Letters . vol.  1,  no.  4,  pp. 
121-125. 

Bartholdi,  J.J.  and  L.K.  Platzman  [1988],  "Heuristics  Based 
on  Spacefilling  Curves  for  Combinatorial  Problems  in 
Euclidean  Space,"  Management  Science,  vol.  34,  no.  3,  pp. 
291-305. 

Bartholdi,  J.J.  Ill,  L.K.  Platzman,  R.L . Collins,  and  W.H. 
Warden,  III  [1983],  "A  Minimal  Technology  Routing  System  for 
Meals  on  Wheels,"  Interfaces . vol.  13,  no.  3,  pp.  1-8. 

Bedini,  R. , G.  G.  Lisini  and  P.  Sterpos  [1979],  "Optimal 
Programming  of  Working  Cycles  for  Industrial  Robots," 

Journal  of  Mechanical  Design,  vol.  101,  pp.  250-257. 

Bellmore,  M.  and  G.  L.  Nemhauser  [1968],  The  Traveling 
Salesman  Problem:  A Survey,"  Operations  Research,  vol.  16, 
pp.  538-558. 

Bondy,  J . A.  and  U.S.R.  Murty  [1985],  "Graph  Theory  with 
Applications  f " Elsevier  Science  Publishing,  New  York,  NY. 

Bozer,  Y.A.,  E.C.  Schorn,  and  G.P.  Sharp  [1990],  "Geometric 
Approaches  to  Solve  the  Chebyshev  Traveling  Salesman 
Problem,"  HE  Transactions,  vol.  22,  no.  3,  pp.  238-254. 


156 


157 


Cavalier,  T.  M.  and  N.  A.  Younis  [1986],  "A  Robotic  Assembly 
Layout  Problem,"  ORSA/TIMS  Bulletin,  Miami  Beach  Meeting. 

Chalmet,  L.G.,  R.L . Francis,  and  A.  Kolen  [1981],  "Finding 
Efficient  Solutions  for  Rectilinear  Distance  Location 
Problems  Efficiently,"  European  Journal  of  Operational 
Research,  vol.  6,  pp.  117,  124. 

Chauny,  F.,  A.  Haurie,  E.  Wagneur,  and  R.  Loulou  [1987], 
"Sequencing  Punch  Operations  in  a Flexible  Manufacturing 
Cell  A Three-Dimensional  Space-Filling  Curve,"  INFOR,  vol. 
25,  no.  1,  pp.  26-45. 

Chern,  T.  C.,  W.  L.  Hsu  and  W.  K.  Shih  [1989],  "An 
0(n2  log  n)  Algorithm  for  the  Hamiltonian  Cycle  Problem  on 
Circular-Arc  Graphs,"  CORS/TIMS/ORSA  Bulletin,  Vancouver, 
Canada. 

Chisman,  J.A.  [1975],  "The  Clustered  Traveling  Salesman 
Problem,"  Computers  and  Operations  Research,  vol.  2,  pp. 
115-119. 

Croes,  G. A.  [1958],  "A  Method  for  Solving  Traveling-Salesman 
Problems,"  Operations  Research,  vol.  6,  pp.  791-812. 

Dieudonne,  J [1969],  Foundations  of  Modern  Analysis. 

Academic  Press,  New  York,  NY. 

Drezner,  Z.  and  S.  Y.  Nof  [1984],  "On  Optimizing  Bin  Picking 
and  Insertion  Plans  for  Assembly  Robots,"  HE  Transactions, 
vol.  16,  no.  3,  pp.  262-270. 

Francis,  R.L.,  H.W.  Hamacher,  C.-Y.  Lee,  and  S.  Yerelan 
[1994],  "On  Automating  Robotic  Assembly  Workplace  Planning," 
to  appear  in  IIE  Transactions,  Research  Report  89-7, 
Department  of  Industrial  & Systems  Engineering,  University 
of  Florida,  Gainesville,  FL. 

Francis,  R.L.,  L.F.  McGinnis,  Jr,,  and  J.A.  White  [1992], 
"Facility  Layout  and  Location;  Analytical  Approach. " 

Prentice  Hall,  Englewood  Cliffs,  NJ. 

Garey,  M.R.  and  D.S.  Johnson  [1979],  Computers  and 
Intractability:  A Guide  to  the  Theory  of  NP-Completness . 
Freeman,  San  Francisco,  CA. 

Hu,  X.  and  M.  W.  Carter  [1989],  "Optimal  Component 
Sequencing  to  Minimize  Total  Production  Time,"  presented  at 
CORS/TIMS/ORSA  Meeting  in  Vancouver,  Department  of 
Industrial  Engineering,  University  of  Toronto,  Toronto, 
Ontario,  Canada. 


158 


Itai,  A.,  C.H.  Papadimitriou  and  J.L.  Szwarcfiter  [1982], 
"Hamilton  Paths  in  Grid  Graphs",  SIAM  Journal  on  Computing, 
vol.  11,  no.  4,  pp.  676-686. 


Lawler,  E.  L,  J.  K.  Lenstra,  A.  H.  G.  Rinnooy  Kan  and  D.  B. 
Shmoys  [1985],  The  Traveling  Salesman  Problem;  a Guided  Tour 
of  Combinatorial  Optimization.  John  Wiley  & Sons,  New  York, 
NY. 

Leipala,  T.  and  0.  Nevalainen  [1989],  "Optimization  of  the 
Movements  of  a Component  Placement  Machine,"  European 
Journal  of  Operational  Research,  vol.  38,  pp.  167-177. 

Lofgren,  C.  B.,  and  L.  F.  McGinnis  [1986a],  "Dynamic 
Scheduling  for  Flexible  Printed  Circuit  Card  Assembly," 
Proceedings  of  the  1986  IEEE  Conference  on  Systems.  Man  and 
Cybernetics . Atlanta,  GA. 

Lofgren,  C.  B.,  and  L.  F.  McGinnis  [1986b],  "Soft 
Configuration  in  Automated  Insertion,"  Proceedings  of  the 
1986  IEEE  International  Conference  on  Robotics  and 
Automation . San  Francisco,  CA. 

Love,  R.F.,  J.G.  Morris,  and  G.O.  Weselowsky  [1988], 
Facilities  Locations;  Models  & Methods.  North-Holland,  New 
York,  NY. 

Luccio,  F.  and  C.  Mugnai  [1978],  "Hamiltonian  Paths  on  a 
Rectangular  Chessboard,"  Proceedings  of  the  16th  Annual 
Allerton  Conference,  pp.  73-78. 

Nof,  S.  Y.  and  Z.  Drezner  [1986],  "Part  Flow  in  the  Robotic 
Assembly  Plan  Problem,"  Material  Flow,  vol.  3,  pp.  197-205. 

Papadimitriou,  C.  H.  and  K.  Steiglitz  [1982],  Combinatorial 
Optimization;  Algorithms  and  Complexity.  Prentice  Hall  Inc., 
Englewood  Cliffs,  NJ. 

Sarin,  S.  C.  [1987],  "Mathematical  Analysis  of  a Robotized 
Production  Cell,"  INFOR.  vol.  25,  no.  1,  pp.  46-56. 

Sarin,  S.  C.  and  W.  E.  Wilhelm  [1984],  "Prototype  Models  for 
Two-Dimensional  Layout  Design  for  Robot  Systems,"  HE 
Transactions . vol.  16,  no.  3,  pp.  206-215. 

Simeone,  B.,  P.  Toth,  G.  Gallo,  F.  Mafioli  and  S.  Pallotino 
[1988],  "FORTRAN  Codes  for  Network  Optimization,"  Annals  of 
Operations  Research.  13 . 

Stewart,  W.R.  Jr.  [1977],  "A  Computationally  Efficient 
Heuristic  for  the  Traveling  Salesman  Problem, " Proceedings 
13th  Annual  Meeting  of  Southeast  TIMS,  pp.  75-85. 


159 


Tang,  C.S.  [1986],  "A  Job  Scheduling  Model  for  a Flexible 
Manufacturing  Machine,"  Proceedings  of  the  1986  IEEE 
International  Conference  on  Robotics  and  Automation,  San 
Francisco,  CA. 

Walas,  R.  A.  and  R.  G.  Askin  [1984],  "An  Algorithm  for  NC 
Turret  Punch  Press  Tool  Location  and  Hit  Sequencing,"  HE 
Transactions . vol.  16,  no.  3,  pp.  280-287. 

Wilhelm,  W.  E.  and  S.  C.  Sarin  [1985],  "A  Structure  for 
Sequencing  Robot  Activities  in  Machine  Tending 
Applications,"  International  Journal  of  Production  Research, 
vol.  23,  pp.  47-64. 

Younis,  N.  A.  and  T.  M.  Cavalier  [1989],  "On  Locating  Part 
Bins  in  a Constrained  Layout  Area  for  an  Automated  Assembly 
Process,"  Computers  and  Industrial  Engineering,  preliminary 
copy. 


BIOGRAPHICAL  SKETCH 

Tomas  Horak  was  born  on  March  17,  1960,  in  Brno, 
Czechoslovakia.  He  is  currently  studying  for  the  Doctor  of 
Philosophy  degree  which  he  expects  to  receive  in  August 
1994.  He  received  his  Bachelor  of  Science  in  mathematics  in 
December,  1982,  Master  of  Science  in  mathematics  in  August, 
1985,  and  Master  of  Science  in  industrial  and  systems 
engineering  in  May,  1988.  He  received  all  degrees  from  the 
University  of  Florida. 


160 


I certify  that  I have  read  this  study  and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality, 
as  a dissertation  for  the  degree  of  Doctor  of  Philosophy. 


Richard  L.  Francis 
Professor  of  Industrial 


& Systems  Engineering 


I certify  that  I have  read  this  study  and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality, 
as  a dissertation  for  the  degree  of  Doctor  of  Philosophy. 


Donald  w.  Hearn 
Professor  of  Industrial 


& Systems  Engineering 


I certify  that  I have  read  this  study  and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality, 
as  a dissertation  for  the  degree  of  Doctor  of  Philosophy. 

Chung-Yee  Lee 

Associate  Professor  of  Industrial 
& Systems  Engineering 


I certify  that  I have  read  this  study  and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality, 
as  a dissertation  for  the  degree  of  Doctor  of  ^Philosophy . 


Suleyman  'Fufekci 
Associate  Professor  of  Industrial 
& Systems  Engineering 


I certify  that  I have  read  this  study  and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality, 
as  a dissertation  for  the  degree  of  Doctor  of  Philosophy. 


This  dissertation  was  submitted  to  the  Graduate 
Faculty  of  the  College  of  Engineering  and  to  the  Graduate 
School  and  was  accepted  as  partial  fulfillment  of  the 
requirements  for  the  f ' “ " Uosophy. 


August  1994 


Professor' of  Decision 

& Information  Sciences 


Dean,  College  of  Engineering 


Karen  A.  Holbrook 
Dean,  Graduate  School 


